20 br_tautau(19961, 2, 0.),
22 br_mumu(19961, 2, 0.),
25 GammaHtot_SM(19961, 2, 0.),
26 log_cs_ggH_8(199, 2, 0.),
27 log_cs_VBF_8(199, 2, 0.),
28 log_cs_WH_8(199, 2, 0.),
29 log_cs_ZH_8(199, 2, 0.),
30 log_cs_ggH_13(199, 2, 0.),
31 log_cs_VBF_13(199, 2, 0.),
32 log_cs_WH_13(199, 2, 0.),
33 log_cs_ZH_13(199, 2, 0.),
34 log_cs_ttH_8(199, 2, 0.),
35 log_cs_ttH_13(199, 2, 0.),
36 log_cs_bbH_8(199, 2, 0.),
37 log_cs_bbH_13(199, 2, 0.),
38 log_cs_ggA_8(199, 2, 0.),
39 log_cs_ttA_8(199, 2, 0.),
40 log_cs_bbA_8(199, 2, 0.),
41 log_cs_ggA_13(199, 2, 0.),
42 log_cs_ttA_13(199, 2, 0.),
43 log_cs_bbA_13(199, 2, 0.),
44 log_cs_ggHp_8(744, 3, 0.),
45 log_cs_ggHp_13(1104, 3, 0.),
46 csrH_top_8(199, 2, 0.),
47 csrH_bottom_8(199, 2, 0.),
48 csrA_top_8(199, 2, 0.),
49 csrA_bottom_8(199, 2, 0.),
50 csrH_top_13(199, 2, 0.),
51 csrH_bottom_13(199, 2, 0.),
52 csrA_top_13(199, 2, 0.),
53 csrA_bottom_13(199, 2, 0.),
54 ATLAS8_pp_phi_gaga(108, 2, 0.),
55 ATLAS8_pp_phi_Zga_llga(141, 2, 0.),
56 ATLAS8_gg_phi_tautau(92, 2, 0.),
57 ATLAS8_bb_phi_tautau(92, 2, 0.),
58 ATLAS8_gg_A_hZ_tautauZ(79, 2, 0.),
59 ATLAS8_gg_A_hZ_bbZ(79, 2, 0.),
60 ATLAS8_gg_phi_tt(53, 2, 0.),
61 ATLAS8_gg_H_WW(13,2,0.),
62 ATLAS8_VBF_H_WW(13,2,0.),
63 ATLAS8_gg_H_ZZ(173,2,0.),
64 ATLAS8_VBF_H_ZZ(173,2,0.),
65 ATLAS8_gg_H_hh(75,2,0.),
66 ATLAS8_pp_phi_gaga_e(108, 2, 0.),
67 ATLAS8_pp_phi_Zga_llga_e(141, 2, 0.),
68 ATLAS8_gg_phi_tautau_e(92, 2, 0.),
69 ATLAS8_bb_phi_tautau_e(92, 2, 0.),
70 ATLAS8_gg_A_hZ_tautauZ_e(79, 2, 0.),
71 ATLAS8_gg_A_hZ_bbZ_e(79, 2, 0.),
72 ATLAS8_gg_phi_tt_e(53, 2, 0.),
73 ATLAS8_gg_H_WW_e(13,2,0.),
74 ATLAS8_VBF_H_WW_e(13,2,0.),
75 ATLAS8_gg_H_ZZ_e(173,2,0.),
76 ATLAS8_VBF_H_ZZ_e(173,2,0.),
77 ATLAS8_gg_H_hh_e(75,2,0.),
78 CMS8_mu_pp_H_VV(172, 2, 0.),
79 CMS8_mu_pp_H_VV_e(172, 2, 0.),
80 CMS8_gg_A_hZ_bbll(16, 2, 0.),
81 CMS8_pp_H_hh(71, 2, 0.),
82 CMS8_pp_H_hh_gagabb(85, 2, 0.),
83 CMS8_pp_H_hh_bbbb(167, 2, 0.),
84 CMS8_bb_phi_bb(81, 2, 0.),
85 CMS8_gg_phi_tautau(92,2,0.),
86 CMS8_bb_phi_tautau(92,2,0.),
87 CMS8_gg_phi_gaga(141,2,0.),
88 CMS8_pp_A_Zga_llga(101,2,0.),
89 CMS8_pp_phi_Zga(101,2,0.),
90 CMS8_gg_H_hh_bbtautau(10,2,0.),
91 CMS8_gg_A_hZ_tautaull(14,2,0.),
92 CMS8_pp_A_HZ_bbll(28718, 3, 0.),
93 CMS8_pp_H_AZ_bbll(29050, 3, 0.),
94 CMS8_pp_A_HZ_tautaull(400, 3, 0.),
95 CMS8_pp_H_AZ_tautaull(400, 3, 0.),
96 CMS8_gg_A_hZ_bbll_e(16, 2, 0.),
97 CMS8_pp_H_hh_e(71, 2, 0.),
98 CMS8_pp_H_hh_gagabb_e(85, 2, 0.),
99 CMS8_pp_H_hh_bbbb_e(167, 2, 0.),
100 CMS8_bb_phi_bb_e(81, 2, 0.),
101 CMS8_gg_phi_tautau_e(92,2,0.),
102 CMS8_bb_phi_tautau_e(92,2,0.),
103 CMS8_gg_phi_gaga_e(141,2,0.),
104 CMS8_pp_A_Zga_llga_e(101,2,0.),
105 CMS8_gg_H_hh_bbtautau_e(10,2,0.),
106 CMS8_gg_A_hZ_tautaull_e(14,2,0.),
111 ATLAS13_bb_phi_tt(61,2,0.),
112 ATLAS13_tt_phi_tt(61,2,0.),
113 ATLAS13_gg_phi_tautau(206,2,0.),
114 ATLAS13_bb_phi_tautau(206,2,0.),
115 ATLAS13_pp_phi_gaga(251,2,0.),
116 ATLAS13_pp_phi_Zga(216,2,0.),
117 ATLAS13_gg_phi_Zga_llga(216,2,0.),
118 ATLAS13_gg_H_ZZ_llllnunu(101,2,0.),
119 ATLAS13_VBF_H_ZZ_llllnunu(101,2,0.),
120 ATLAS13_gg_H_ZZ_llnunu(71,2,0.),
121 ATLAS13_gg_H_ZZ_llll(81,2,0.),
122 ATLAS13_VBF_H_ZZ_llll(81,2,0.),
123 ATLAS13_gg_H_ZZ_qqllnunu(271,2,0.),
124 ATLAS13_VBF_H_ZZ_qqllnunu(271,2,0.),
125 ATLAS13_gg_H_ZZ_llqq(271,2,0.),
126 ATLAS13_VBF_H_ZZ_llqq(271,2,0.),
127 ATLAS13_gg_H_ZZ_nunuqq(251,2,0.),
128 ATLAS13_gg_H_WW_enumumu(381,2,0.),
129 ATLAS13_VBF_H_WW_enumumu(281,2,0.),
130 ATLAS13_gg_H_WW_lnuqq(271,2,0.),
131 ATLAS13_VBF_H_WW_lnuqq(271,2,0.),
132 ATLAS13_pp_H_VV_qqqq(181,2,0.),
133 ATLAS13_pp_H_hh_bbbb(271,2,0.),
134 ATLAS13_pp_H_hh_gagabb(26,2,0.),
135 ATLAS13_pp_H_hh_gagaWW(25,2,0.),
136 ATLAS13_gg_A_Zh_Zbb(181,2,0.),
137 ATLAS13_bb_A_Zh_Zbb(181,2,0.),
138 ATLAS13_bb_phi_tt_e(61,2,0.),
139 ATLAS13_tt_phi_tt_e(61,2,0.),
140 ATLAS13_gg_phi_tautau_e(206,2,0.),
141 ATLAS13_bb_phi_tautau_e(206,2,0.),
142 ATLAS13_pp_phi_gaga_e(251,2,0.),
143 ATLAS13_pp_phi_Zga_e(216,2,0.),
144 ATLAS13_gg_phi_Zga_llga_e(216,2,0.),
145 ATLAS13_gg_H_ZZ_llllnunu_e(101,2,0.),
146 ATLAS13_VBF_H_ZZ_llllnunu_e(101,2,0.),
147 ATLAS13_gg_H_ZZ_llnunu_e(71,2,0.),
148 ATLAS13_gg_H_ZZ_llll_e(81,2,0.),
149 ATLAS13_VBF_H_ZZ_llll_e(81,2,0.),
150 ATLAS13_gg_H_ZZ_qqllnunu_e(271,2,0.),
151 ATLAS13_VBF_H_ZZ_qqllnunu_e(271,2,0.),
152 ATLAS13_gg_H_ZZ_llqq_e(271,2,0.),
153 ATLAS13_VBF_H_ZZ_llqq_e(271,2,0.),
154 ATLAS13_gg_H_ZZ_nunuqq_e(251,2,0.),
155 ATLAS13_gg_H_WW_enumumu_e(381,2,0.),
156 ATLAS13_VBF_H_WW_enumumu_e(281,2,0.),
157 ATLAS13_gg_H_WW_lnuqq_e(271,2,0.),
158 ATLAS13_VBF_H_WW_lnuqq_e(271,2,0.),
159 ATLAS13_pp_H_VV_qqqq_e(181,2,0.),
160 ATLAS13_pp_H_hh_bbbb_e(271,2,0.),
161 ATLAS13_pp_H_hh_gagabb_e(26,2,0.),
162 ATLAS13_pp_H_hh_gagaWW_e(25,2,0.),
163 ATLAS13_gg_A_Zh_Zbb_e(181,2,0.),
164 ATLAS13_bb_A_Zh_Zbb_e(181,2,0.),
165 CMS13_pp_phi_bb(66,2,0.),
166 CMS13_gg_phi_tautau(312,2,0.),
167 CMS13_bb_phi_tautau(312,2,0.),
168 CMS13_gg_phi_gaga(351,2,0.),
169 CMS13_pp_phi_Zga_llga(171,2,0.),
170 CMS13_pp_phi_Zga_qqga(236,2,0.),
171 CMS13_ggF_phi_Zga(366,2,0.),
172 CMS13_pp_H_ZZ_llnunu(191,2,0.),
173 CMS13_gg_H_ZZ_llnunu(131,2,0.),
174 CMS13_VBF_H_ZZ_llnunu(131,2,0.),
175 CMS13_pp_H_ZZ_llll(241,2,0.),
176 CMS13_VBFVH_H_ZZ_llll(241,2,0.),
177 CMS13_pp_H_ZZ_llqq(151,2,0.),
178 CMS13_ggFVBF_H_WW_lnulnu(81,2,0.),
179 CMS13_pp_H_hh_bbbb(95,2,0.),
180 CMS13_ggF_H_hh_bbbb(226,2,0.),
181 CMS13_pp_H_hh_gagabb(66,2,0.),
182 CMS13_pp_H_hh_bbtautau(66,2,0.),
183 CMS13_pp_H_hh_bbtautau1(66,2,0.),
184 CMS13_pp_H_hh_bblnulnu(65,2,0.),
185 CMS13_pp_H_hh_bbVV(65,2,0.),
186 CMS13_pp_phi_bb_e(66,2,0.),
187 CMS13_gg_phi_tautau_e(312,2,0.),
188 CMS13_bb_phi_tautau_e(312,2,0.),
189 CMS13_gg_phi_gaga_e(351,2,0.),
190 CMS13_pp_phi_Zga_llga_e(171,2,0.),
191 CMS13_pp_phi_Zga_qqga_e(236,2,0.),
192 CMS13_ggF_phi_Zga_e(366,2,0.),
193 CMS13_pp_H_ZZ_llnunu_e(191,2,0.),
194 CMS13_gg_H_ZZ_llnunu_e(131,2,0.),
195 CMS13_VBF_H_ZZ_llnunu_e(131,2,0.),
196 CMS13_pp_H_ZZ_llll_e(241,2,0.),
197 CMS13_VBFVH_H_ZZ_llll_e(241,2,0.),
198 CMS13_pp_H_ZZ_llqq_e(151,2,0.),
199 CMS13_ggFVBF_H_WW_lnulnu_e(81,2,0.),
200 CMS13_pp_H_hh_bbbb_e(95,2,0.),
201 CMS13_ggF_H_hh_bbbb_e(226,2,0.),
202 CMS13_pp_H_hh_gagabb_e(66,2,0.),
203 CMS13_pp_H_hh_bbtautau_e(66,2,0.),
204 CMS13_pp_H_hh_bbtautau1_e(66,2,0.),
205 CMS13_pp_H_hh_bblnulnu_e(65,2,0.),
206 CMS13_pp_H_hh_bbVV_e(65,2,0.),
207 temp1(1,2,0.), temp2(1,2,0.), temp3(1,2,0.), temp4(1,2,0.), temp5(1,2,0.), temp6(1,2,0.), temp7(1,2,0.), temp8(1,2,0.), temp9(1,2,0.), temp10(1,2,0.),
208 temp11(1,2,0.), temp12(1,2,0.), temp13(1,2,0.), temp14(1,2,0.), temp15(1,2,0.), temp16(1,2,0.), temp17(1,2,0.), temp18(1,2,0.), temp19(1,2,0.), temp20(1,2,0.),
209 temp21(1,2,0.), temp22(1,2,0.), temp23(1,2,0.), temp24(1,2,0.), temp25(1,2,0.), temp26(1,2,0.), temp27(1,2,0.), temp28(1,2,0.), temp29(1,2,0.), temp30(1,2,0.),
210 temp31(1,2,0.), temp32(1,2,0.), temp33(1,2,0.), temp34(1,2,0.), temp35(1,2,0.), temp36(1,2,0.), temp37(1,2,0.), temp38(1,2,0.), temp39(1,2,0.), temp40(1,2,0.),
211 temp1e(1,2,0.), temp2e(1,2,0.), temp3e(1,2,0.), temp4e(1,2,0.), temp5e(1,2,0.), temp6e(1,2,0.), temp7e(1,2,0.), temp8e(1,2,0.), temp9e(1,2,0.), temp10e(1,2,0.),
212 temp11e(1,2,0.), temp12e(1,2,0.), temp13e(1,2,0.), temp14e(1,2,0.), temp15e(1,2,0.), temp16e(1,2,0.), temp17e(1,2,0.), temp18e(1,2,0.), temp19e(1,2,0.), temp20e(1,2,0.),
213 temp21e(1,2,0.), temp22e(1,2,0.), temp23e(1,2,0.), temp24e(1,2,0.), temp25e(1,2,0.), temp26e(1,2,0.), temp27e(1,2,0.), temp28e(1,2,0.), temp29e(1,2,0.), temp30e(1,2,0.),
214 temp31e(1,2,0.), temp32e(1,2,0.), temp33e(1,2,0.), temp34e(1,2,0.), temp35e(1,2,0.), temp36e(1,2,0.), temp37e(1,2,0.), temp38e(1,2,0.), temp39e(1,2,0.), temp40e(1,2,0.),
215 ATLAS8_pp_Hpm_taunu(83,2,0.),
216 ATLAS8_pp_Hpm_tb(41,2,0.),
217 ATLAS8_pp_Hpm_taunu_e(83,2,0.),
218 ATLAS8_pp_Hpm_tb_e(41,2,0.),
219 CMS8_pp_Hp_taunu(43,2,0.),
220 CMS8_pp_Hp_tb(43,2,0.),
221 CMS8_pp_Hp_taunu_e(43,2,0.),
222 CMS8_pp_Hp_tb_e(43,2,0.),
223 ATLAS13_pp_Hpm_taunu(181,2,0.),
224 ATLAS13_pp_Hp_tb1(71,2,0.),
225 ATLAS13_pp_Hp_tb2(181,2,0.),
226 ATLAS13_pp_Hpm_taunu_e(181,2,0.),
227 ATLAS13_pp_Hp_tb1_e(71,2,0.),
228 ATLAS13_pp_Hp_tb2_e(181,2,0.),
229 CMS13_pp_Hpm_taunu(283,2,0.),
230 CMS13_pp_Hpm_taunu_e(283,2,0.),
231 arraybsgamma(1111, 3, 0.),
232 myTHDM(static_cast<const
THDM*> (&SM_i)),
261 const int NumPar,
const double params[])
const {
265 for(
int j=0; j<NumPar; j++)
266 bCache &= (params[j] == cache[j][i].real());
267 if (bCache)
return i;
273 const int NumPar,
const double params[])
const {
277 for(
int j=0; j<NumPar; j++)
278 bCache &= (params[j] == cache[j][i]);
279 if (bCache)
return i;
288 for(
int j=0; j<NumPar+1; j++)
289 cache[j][i] = cache[j][i-1];
292 for(
int j=0; j<NumPar; j++) {
294 cache[NumPar][0] = newResult;
299 const double params[],
const double newResult)
const {
302 for(
int j=0; j<NumPar+1; j++)
303 cache[j][i] = cache[j][i-1];
306 for(
int j=0; j<NumPar; j++) {
307 cache[j][0] = params[j];
308 cache[NumPar][0] = newResult;
318 double params[] = {MZ2, MW2,
mHh2};
332 double params[] = {MZ2, MW2,
mHl2};
346 double params[] = {MZ2,
mHh2};
360 double params[] = {MZ2,
mHl2};
374 double params[] = {MZ2, MW2,
mHh2};
388 double params[] = {MZ2, MW2,
mHl2};
402 double params[] = {MZ2,
mHh2};
416 double params[] = {MZ2,
mHl2};
430 double params[] = {MZ2,
mHl2};
444 double params[] = {MZ2,
mHh2};
458 double params[] = {MZ2,
mHp2,
mHl2};
472 double params[] = {MZ2,
mHp2,
mHh2};
486 double params[] = {MZ2,
mA2,
mHl2};
500 double params[] = {MZ2,
mA2,
mHh2};
514 double params[] = {MZ2,
mHl2};
528 double params[] = {MZ2,
mHl2,
mHp2};
542 double params[] = {MZ2,
mHl2,
mA2};
556 double params[] = {MZ2,
mHl2};
570 double params[] = {MZ2,
mHl2,
mHh2};
584 double params[] = {MZ2,
mHl2,
mHh2};
598 double params[] = {MZ2,
mHl2,
mHp2};
612 double params[] = {MZ2,
mHl2,
mA2};
626 double params[] = {MZ2,
mHh2};
640 double params[] = {MZ2,
mHh2,
mHp2};
654 double params[] = {MZ2,
mHh2,
mA2};
668 double params[] = {MZ2,
mHh2,
mHl2};
682 double params[] = {MZ2,
mHh2,
mHl2};
696 double params[] = {MZ2,
mHh2};
710 double params[] = {MZ2,
mHh2,
mHp2};
724 double params[] = {MZ2,
mHh2,
mA2};
738 double params[] = {MZ2,
mHp2,
mHl2};
752 double params[] = {MZ2,
mHp2,
mHh2};
766 double params[] = {MZ2,
mHp2,
mHl2};
780 double params[] = {MZ2,
mHp2,
mHh2};
794 double params[] = {MZ2,
mA2,
mHl2};
808 double params[] = {MZ2,
mA2,
mHh2};
822 double params[] = {MZ2,
mA2,
mHl2};
836 double params[] = {MZ2,
mA2,
mHh2};
852 double params[] = {MZ2,
mHl2};
866 double params[] = {MZ2,
mHh2};
880 double params[] = {MZ2,
mHp2,
mHl2};
894 double params[] = {MZ2,
mHp2,
mHh2};
908 double params[] = {MZ2,
mHp2,
mA2};
922 double params[] = {MZ2,
mA2,
mHl2};
936 double params[] = {MZ2,
mA2,
mHh2};
950 double params[] = {MZ2,
mHl2};
964 double params[] = {MZ2,
mHl2,
mHp2};
978 double params[] = {MZ2,
mHl2,
mA2};
992 double params[] = {MZ2,
mHl2};
1006 double params[] = {MZ2,
mHl2,
mHh2};
1020 double params[] = {MZ2,
mHl2,
mHh2};
1034 double params[] = {MZ2,
mHl2,
mHp2};
1048 double params[] = {MZ2,
mHl2,
mA2};
1062 double params[] = {MZ2,
mHh2};
1076 double params[] = {MZ2,
mHh2,
mHp2};
1090 double params[] = {MZ2,
mHh2,
mA2};
1104 double params[] = {MZ2,
mHh2,
mHl2};
1118 double params[] = {MZ2,
mHh2,
mHl2};
1132 double params[] = {MZ2,
mHh2};
1146 double params[] = {MZ2,
mHh2,
mHp2};
1160 double params[] = {MZ2,
mHh2,
mA2};
1174 double params[] = {MZ2,
mHp2,
mHl2};
1188 double params[] = {MZ2,
mHp2,
mHh2};
1202 double params[] = {MZ2,
mHp2,
mA2};
1216 double params[] = {MZ2,
mHp2,
mHl2};
1230 double params[] = {MZ2,
mHp2,
mHh2};
1244 double params[] = {MZ2,
mA2,
mHl2};
1258 double params[] = {MZ2,
mA2,
mHh2};
1272 double params[] = {MZ2,
mA2,
mHp2};
1286 double params[] = {MZ2,
mA2,
mHl2};
1300 double params[] = {MZ2,
mA2,
mHh2};
1316 double params[] = {MZ2,
mA2,
mHp2};
1330 double params[] = {MZ2,
mHh2,
mA2};
1344 double params[] = {MZ2,
mHh2,
mHp2};
1358 double params[] = {MZ2,
mHl2,
mA2};
1372 double params[] = {MZ2,
mHl2,
mHp2};
1386 double params[] = {MZ2,
mHp2};
1400 double params[] = {MZ2, MW2,
mHh2};
1414 double params[] = {MZ2, MW2,
mHl2};
1428 double params[] = {MZ2,
mHh2};
1442 double params[] = {MZ2,
mHl2};
1456 double params[] = {MZ2, MW2,
mA2,
mHp2};
1470 double params[] = {MZ2, MW2,
mHh2,
mHp2};
1484 double params[] = {MZ2, MW2,
mHl2,
mHp2};
1498 double params[] = {MZ2, MW2,
mHp2};
1512 double params[] = {MZ2, MW2,
mHh2};
1526 double params[] = {MZ2, MW2,
mHl2};
1540 double params[] = {MZ2,
mHh2,
mA2};
1554 double params[] = {MZ2,
mHl2,
mA2};
1568 double params[] = {MZ2,
mHp2};
1582 double params[] = {MZ2,
mHh2};
1596 double params[] = {MZ2,
mHl2};
1610 std::stringstream br1,br2,br3,br4,br5,br6,br7;
1611 std::stringstream dw1;
1612 std::stringstream cs1,cs2,cs3,cs4,cs5,cs6,cs7,cs8,cs9;
1613 std::stringstream cs11,cs12,cs13,cs14,cs15,cs16,cs17,cs18,cs19;
1614 std::stringstream cs20,cs21;
1615 std::stringstream csr1,csr2,csr3,csr4;
1616 std::stringstream csr11,csr12,csr13,csr14;
1617 std::stringstream ex1,ex2,ex3,ex4,ex5,ex6,ex7,ex8,ex9,ex10,ex11,ex12,ex13,ex14,ex15,ex16,ex17,ex18,ex19,ex20,ex21,ex22,ex23;
1618 std::stringstream ex1e,ex2e,ex3e,ex4e,ex5e,ex6e,ex7e,ex8e,ex9e,ex10e,ex11e,ex12e,ex13e,ex14e,ex15e,ex16e,ex17e,ex18e,ex19e,ex20e,ex21e,ex22e,ex23e;
1620 std::stringstream ex24,ex25,ex26,ex27,ex28,ex29,ex30,ex31,ex32,ex33,ex34,ex35,ex36,ex37,ex38,ex39,ex40,ex41,ex42,ex43,ex44,\
1621 ex45,ex46,ex47,ex48,ex49,ex50,ex51,ex52,ex53,ex54,ex55,ex56;
1622 std::stringstream ex24e,ex25e,ex26e,ex27e,ex28e,ex29e,ex30e,ex31e,ex32e,ex33e,ex34e,ex35e,ex36e,ex37e,ex38e,ex39e,ex40e,ex41e,ex42e,ex43e,ex44e,\
1623 ex45e,ex46e,ex47e,ex48e,ex49e,ex50e,ex51e,ex52e,ex53e,ex54e,ex55e,ex56e;
1624 std::stringstream ex57,ex58,ex59,ex60,ex61,ex62,ex63,ex64,ex65,ex66,ex67,ex68,ex69,ex70,ex71,ex72,ex73,ex74,ex75,ex76,ex77,\
1625 ex78,ex79,ex80,ex81,ex82,ex83,ex84,ex85,ex86,ex87,ex88,ex89,ex90,ex91,ex92,ex93,ex94,ex95,ex96,ex97,ex98;
1626 std::stringstream ex57e,ex58e,ex59e,ex60e,ex61e,ex62e,ex63e,ex64e,ex65e,ex66e,ex67e,ex68e,ex69e,ex70e,ex71e,ex72e,ex73e,ex74e,ex75e,ex76e,ex77e,\
1627 ex78e,ex79e,ex80e,ex81e,ex82e,ex83e,ex84e,ex85e,ex86e,ex87e,ex88e,ex89e,ex90e,ex91e,ex92e,ex93e,ex94e,ex95e,ex96e,ex97e,ex98e;
1628 std::stringstream ex99,ex100,ex101,ex102,ex103,ex104,ex105,ex106,ex107,ex108,ex109,ex110,ex111,ex112,ex113,ex114,ex115,ex116,ex117,ex118,ex119,ex120;
1629 std::stringstream ex99e,ex100e,ex101e,ex102e,ex103e,ex104e,ex105e,ex106e,ex107e,ex108e,ex109e,ex110e,ex111e,ex112e,ex113e,ex114e,ex115e,ex116e,ex117e,ex118e,ex119e,ex120e;
1630 std::stringstream ex121,ex122,ex123,ex124;
1631 std::stringstream bsg1;
1633 std::cout<<
"reading tables"<<std::endl;
1636 std::stringstream path;
1637 path << getenv(
"HEPFITTABS") <<
"/THDM/tabs/";
1638 std::string tablepath=path.str();
1640 br1 << tablepath <<
"br1.dat";
1655 br2 << tablepath <<
"br2.dat";
1657 br3 << tablepath <<
"br3.dat";
1659 br4 << tablepath <<
"br4.dat";
1661 br5 << tablepath <<
"br5.dat";
1663 br6 << tablepath <<
"br6.dat";
1665 br7 << tablepath <<
"br7.dat";
1667 dw1 << tablepath <<
"dw1.dat";
1669 cs1 << tablepath <<
"log_cs_ggH_8.dat";
1671 cs11 << tablepath <<
"log_cs_ggH_13.dat";
1673 cs2 << tablepath <<
"log_cs_VBF_8.dat";
1675 cs12 << tablepath <<
"log_cs_VBF_13.dat";
1677 cs3 << tablepath <<
"log_cs_WH_8.dat";
1679 cs13 << tablepath <<
"log_cs_WH_13.dat";
1681 cs4 << tablepath <<
"log_cs_ZH_8.dat";
1683 cs14 << tablepath <<
"log_cs_ZH_13.dat";
1685 cs5 << tablepath <<
"log_cs_ttH_8.dat";
1687 cs15 << tablepath <<
"log_cs_ttH_13.dat";
1689 cs6 << tablepath <<
"log_cs_bbH_8.dat";
1691 cs16 << tablepath <<
"log_cs_bbH_13.dat";
1693 cs7 << tablepath <<
"log_cs_ggA_8.dat";
1695 cs17 << tablepath <<
"log_cs_ggA_13.dat";
1697 cs8 << tablepath <<
"log_cs_ttA_8.dat";
1699 cs18 << tablepath <<
"log_cs_ttA_13.dat";
1701 cs9 << tablepath <<
"log_cs_bbA_8.dat";
1703 cs19 << tablepath <<
"log_cs_bbA_13.dat";
1705 cs20 << tablepath <<
"log_cs_ggHp_8.dat";
1707 cs21 << tablepath <<
"log_cs_ggHp_13.dat";
1709 csr1 << tablepath <<
"csrH_top_8.dat";
1711 csr11 << tablepath <<
"csrH_top_13.dat";
1713 csr2 << tablepath <<
"csrH_bottom_8.dat";
1715 csr12 << tablepath <<
"csrH_bottom_13.dat";
1717 csr3 << tablepath <<
"csrA_top_8.dat";
1719 csr13 << tablepath <<
"csrA_top_13.dat";
1721 csr4 << tablepath <<
"csrA_bottom_8.dat";
1723 csr14 << tablepath <<
"csrA_bottom_13.dat";
1725 ex1 << tablepath <<
"150400936.dat";
1727 ex1e << tablepath <<
"150400936_e.dat";
1729 ex2 << tablepath <<
"150404710.dat";
1731 ex2e << tablepath <<
"150404710_e.dat";
1733 ex3 << tablepath <<
"160306896.dat";
1735 ex3e << tablepath <<
"160306896_e.dat";
1737 ex4 << tablepath <<
"150304114.dat";
1739 ex4e << tablepath <<
"150304114_e.dat";
1741 ex5 << tablepath <<
"14076583.dat";
1743 ex5e << tablepath <<
"14076583_e.dat";
1745 ex6 << tablepath <<
"14096064_a.dat";
1747 ex6e << tablepath <<
"14096064_a_e.dat";
1749 ex7 << tablepath <<
"14096064_b.dat";
1751 ex7e << tablepath <<
"14096064_b_e.dat";
1753 ex8 << tablepath <<
"150204478_a.dat";
1755 ex8e << tablepath <<
"150204478_a_e.dat";
1757 ex9 << tablepath <<
"150204478_b.dat";
1759 ex9e << tablepath <<
"150204478_b_e.dat";
1761 ex10 << tablepath <<
"150608329.dat";
1763 ex10e << tablepath <<
"150608329_e.dat";
1765 ex11 << tablepath <<
"150507018.dat";
1767 ex11e << tablepath <<
"150507018_e.dat";
1769 ex12 << tablepath <<
"CMS-PAS-HIG-14-029_a.dat";
1771 ex12e << tablepath <<
"CMS-PAS-HIG-14-029_a_e.dat";
1773 ex13 << tablepath <<
"CMS-PAS-HIG-14-029_b.dat";
1775 ex13e << tablepath <<
"CMS-PAS-HIG-14-029_b_e.dat";
1777 ex14 << tablepath <<
"150602301.dat";
1779 ex14e << tablepath <<
"150602301_e.dat";
1793 ex15 << tablepath <<
"150900389_a.dat";
1795 ex15e << tablepath <<
"150900389_a_e.dat";
1797 ex16 << tablepath <<
"150900389_b.dat";
1799 ex16e << tablepath <<
"150900389_b_e.dat";
1801 ex17 << tablepath <<
"150904670.dat";
1803 ex17e << tablepath <<
"150904670_e.dat";
1805 ex18 << tablepath <<
"151001181_a.dat";
1807 ex18e << tablepath <<
"151001181_a_e.dat";
1809 ex19 << tablepath <<
"151001181_b.dat";
1811 ex19e << tablepath <<
"151001181_b_e.dat";
1813 ex20 << tablepath <<
"150705930_a.dat";
1815 ex20e << tablepath <<
"150705930_a_e.dat";
1817 ex21 << tablepath <<
"150705930_b.dat";
1819 ex21e << tablepath <<
"150705930_b_e.dat";
1821 ex22 << tablepath <<
"CMS-PAS-HIG-15-013.dat";
1823 ex22e << tablepath <<
"CMS-PAS-HIG-15-013_e.dat";
1825 ex23 << tablepath <<
"CMS-PAS-HIG-16-014.dat";
1827 ex23e << tablepath <<
"CMS-PAS-HIG-16-014_e.dat";
1830 ex24 << tablepath <<
"ATLAS-CONF-2016-104_b.dat";
1832 ex24e << tablepath <<
"ATLAS-CONF-2016-104_b_e.dat";
1834 ex25 << tablepath <<
"ATLAS-CONF-2016-104_a.dat";
1836 ex25e << tablepath <<
"ATLAS-CONF-2016-104_a_e.dat";
1838 ex26 << tablepath <<
"170907242_a.dat";
1840 ex26e << tablepath <<
"170907242_a_e.dat";
1842 ex27 << tablepath <<
"170907242_b.dat";
1844 ex27e << tablepath <<
"170907242_b_e.dat";
1846 ex28 << tablepath <<
"170704147.dat";
1848 ex28e << tablepath <<
"170704147_e.dat";
1850 ex29 << tablepath <<
"ATLAS-CONF-2016-044.dat";
1852 ex29e << tablepath <<
"ATLAS-CONF-2016-044_e.dat";
1854 ex30 << tablepath <<
"ATLAS-CONF-2016-056.dat";
1856 ex30e << tablepath <<
"ATLAS-CONF-2016-056_e.dat";
1858 ex31 << tablepath <<
"ATLAS-CONF-2016-079_a.dat";
1860 ex31e << tablepath <<
"ATLAS-CONF-2016-079_a_e.dat";
1862 ex32 << tablepath <<
"ATLAS-CONF-2016-079_b.dat";
1864 ex32e << tablepath <<
"ATLAS-CONF-2016-079_b_e.dat";
1866 ex33 << tablepath <<
"ATLAS-CONF-2016-082_a.dat";
1868 ex33e << tablepath <<
"ATLAS-CONF-2016-082_a_e.dat";
1870 ex34 << tablepath <<
"ATLAS-CONF-2016-082_b.dat";
1872 ex34e << tablepath <<
"ATLAS-CONF-2016-082_b_e.dat";
1874 ex35 << tablepath <<
"ATLAS-CONF-2016-082_c.dat";
1876 ex35e << tablepath <<
"ATLAS-CONF-2016-082_c_e.dat";
1878 ex36 << tablepath <<
"171001123_a.dat";
1880 ex36e << tablepath <<
"171001123_a_e.dat";
1882 ex37 << tablepath <<
"171001123_b.dat";
1884 ex37e << tablepath <<
"171001123_b_e.dat";
1886 ex38 << tablepath <<
"171007235_a.dat";
1888 ex38e << tablepath <<
"171007235_a_e.dat";
1890 ex39 << tablepath <<
"ATLAS-CONF-2016-049.dat";
1892 ex39e << tablepath <<
"ATLAS-CONF-2016-049_e.dat";
1894 ex40 << tablepath <<
"ATLAS-CONF-2016-004.dat";
1896 ex40e << tablepath <<
"ATLAS-CONF-2016-004_e.dat";
1898 ex41 << tablepath <<
"ATLAS-CONF-2016-071.dat";
1900 ex41e << tablepath <<
"ATLAS-CONF-2016-071_e.dat";
1902 ex42 << tablepath <<
"171206518_a.dat";
1904 ex42e << tablepath <<
"171206518_a_e.dat";
1906 ex43 << tablepath <<
"171206518_b.dat";
1908 ex43e << tablepath <<
"171206518_b_e.dat";
1910 ex44 << tablepath <<
"CMS-PAS-HIG-16-025.dat";
1912 ex44e << tablepath <<
"CMS-PAS-HIG-16-025_e.dat";
1914 ex45 << tablepath <<
"CMS-PAS-HIG-16-037_a.dat";
1916 ex45e << tablepath <<
"CMS-PAS-HIG-16-037_a_e.dat";
1918 ex46 << tablepath <<
"CMS-PAS-HIG-16-037_b.dat";
1920 ex46e << tablepath <<
"CMS-PAS-HIG-16-037_b_e.dat";
1922 ex47 << tablepath <<
"160902507.dat";
1924 ex47e << tablepath <<
"160902507_e.dat";
1926 ex48 << tablepath <<
"CMS-PAS-EXO-16-034.dat";
1928 ex48e << tablepath <<
"CMS-PAS-EXO-16-034_e.dat";
1930 ex49 << tablepath <<
"CMS-PAS-EXO-16-035.dat";
1932 ex49e << tablepath <<
"CMS-PAS-EXO-16-035_e.dat";
1934 ex50 << tablepath <<
"CMS-PAS-HIG-16-033_a.dat";
1936 ex50e << tablepath <<
"CMS-PAS-HIG-16-033_a_e.dat";
1938 ex51 << tablepath <<
"CMS-PAS-HIG-16-033_b.dat";
1940 ex51e << tablepath <<
"CMS-PAS-HIG-16-033_b_e.dat";
1942 ex52 << tablepath <<
"CMS-PAS-HIG-16-023.dat";
1944 ex52e << tablepath <<
"CMS-PAS-HIG-16-023_e.dat";
1946 ex53 << tablepath <<
"180603548.dat";
1948 ex53e << tablepath <<
"180603548_e.dat";
1950 ex54 << tablepath <<
"180600408.dat";
1952 ex54e << tablepath <<
"180600408_e.dat";
1954 ex55 << tablepath <<
"CMS-PAS-HIG-16-029.dat";
1956 ex55e << tablepath <<
"CMS-PAS-HIG-16-029_e.dat";
1958 ex56 << tablepath <<
"CMS-PAS-HIG-16-011.dat";
1960 ex56e << tablepath <<
"CMS-PAS-HIG-16-011_e.dat";
1963 ex57 << tablepath <<
"t1.dat";
1965 ex57e << tablepath <<
"t1_e.dat";
1967 ex58 << tablepath <<
"t2.dat";
1969 ex58e << tablepath <<
"t2_e.dat";
1971 ex59 << tablepath <<
"t3.dat";
1973 ex59e << tablepath <<
"t3_e.dat";
1975 ex60 << tablepath <<
"t4.dat";
1977 ex60e << tablepath <<
"t4_e.dat";
1979 ex61 << tablepath <<
"t5.dat";
1981 ex61e << tablepath <<
"t5_e.dat";
1983 ex62 << tablepath <<
"t6.dat";
1985 ex62e << tablepath <<
"t6_e.dat";
1987 ex63 << tablepath <<
"t7.dat";
1989 ex63e << tablepath <<
"t7_e.dat";
1991 ex64 << tablepath <<
"t8.dat";
1993 ex64e << tablepath <<
"t8_e.dat";
1995 ex65 << tablepath <<
"t9.dat";
1997 ex65e << tablepath <<
"t9_e.dat";
1999 ex66 << tablepath <<
"t10.dat";
2001 ex66e << tablepath <<
"t10_e.dat";
2003 ex67 << tablepath <<
"t11.dat";
2005 ex67e << tablepath <<
"t11_e.dat";
2007 ex68 << tablepath <<
"t12.dat";
2009 ex68e << tablepath <<
"t12_e.dat";
2011 ex69 << tablepath <<
"t13.dat";
2013 ex69e << tablepath <<
"t13_e.dat";
2015 ex70 << tablepath <<
"t14.dat";
2017 ex70e << tablepath <<
"t14_e.dat";
2019 ex71 << tablepath <<
"t15.dat";
2021 ex71e << tablepath <<
"t15_e.dat";
2023 ex72 << tablepath <<
"t16.dat";
2025 ex72e << tablepath <<
"t16_e.dat";
2027 ex73 << tablepath <<
"t17.dat";
2029 ex73e << tablepath <<
"t17_e.dat";
2031 ex74 << tablepath <<
"t18.dat";
2033 ex74e << tablepath <<
"t18_e.dat";
2035 ex75 << tablepath <<
"t19.dat";
2037 ex75e << tablepath <<
"t19_e.dat";
2039 ex76 << tablepath <<
"t20.dat";
2041 ex76e << tablepath <<
"t20_e.dat";
2043 ex77 << tablepath <<
"t21.dat";
2045 ex77e << tablepath <<
"t21_e.dat";
2047 ex78 << tablepath <<
"t22.dat";
2049 ex78e << tablepath <<
"t22_e.dat";
2051 ex79 << tablepath <<
"t23.dat";
2053 ex79e << tablepath <<
"t23_e.dat";
2055 ex80 << tablepath <<
"t24.dat";
2057 ex80e << tablepath <<
"t24_e.dat";
2059 ex81 << tablepath <<
"t25.dat";
2061 ex81e << tablepath <<
"t25_e.dat";
2063 ex82 << tablepath <<
"t26.dat";
2065 ex82e << tablepath <<
"t26_e.dat";
2067 ex83 << tablepath <<
"t27.dat";
2069 ex83e << tablepath <<
"t27_e.dat";
2071 ex84 << tablepath <<
"t28.dat";
2073 ex84e << tablepath <<
"t28_e.dat";
2075 ex85 << tablepath <<
"t29.dat";
2077 ex85e << tablepath <<
"t29_e.dat";
2079 ex86 << tablepath <<
"t30.dat";
2081 ex86e << tablepath <<
"t30_e.dat";
2083 ex87 << tablepath <<
"t31.dat";
2085 ex87e << tablepath <<
"t31_e.dat";
2087 ex88 << tablepath <<
"t32.dat";
2089 ex88e << tablepath <<
"t32_e.dat";
2091 ex89 << tablepath <<
"t33.dat";
2093 ex89e << tablepath <<
"t33_e.dat";
2095 ex90 << tablepath <<
"t34.dat";
2097 ex90e << tablepath <<
"t34_e.dat";
2099 ex91 << tablepath <<
"t35.dat";
2101 ex91e << tablepath <<
"t35_e.dat";
2103 ex92 << tablepath <<
"t36.dat";
2105 ex92e << tablepath <<
"t36_e.dat";
2107 ex93 << tablepath <<
"t37.dat";
2109 ex93e << tablepath <<
"t37_e.dat";
2111 ex94 << tablepath <<
"t38.dat";
2113 ex94e << tablepath <<
"t38_e.dat";
2115 ex95 << tablepath <<
"t39.dat";
2117 ex95e << tablepath <<
"t39_e.dat";
2119 ex96 << tablepath <<
"t40.dat";
2121 ex96e << tablepath <<
"t40_e.dat";
2124 ex97 << tablepath <<
"CMS-PAS-HIG-16-034.dat";
2126 ex97e << tablepath <<
"CMS-PAS-HIG-16-034_e.dat";
2128 ex98 << tablepath <<
"14078150.dat";
2130 ex98e << tablepath <<
"14078150.dat";
2133 ex99 << tablepath <<
"14126663.dat";
2135 ex99e << tablepath <<
"14126663_e.dat";
2137 ex100 << tablepath <<
"151203704.dat";
2139 ex100e << tablepath <<
"151203704_e.dat";
2141 ex101 << tablepath <<
"150807774_a.dat";
2143 ex101e << tablepath <<
"150807774_a_e.dat";
2145 ex102 << tablepath <<
"150807774_b.dat";
2147 ex102e << tablepath <<
"150807774_b_e.dat";
2149 ex103 << tablepath <<
"ATLAS-CONF-2016-088.dat";
2151 ex103e << tablepath <<
"ATLAS-CONF-2016-088_e.dat";
2153 ex104 << tablepath <<
"ATLAS-CONF-2016-089.dat";
2155 ex104e << tablepath <<
"ATLAS-CONF-2016-089_e.dat";
2157 ex105 << tablepath <<
"ATLAS-CONF-2016-104_c.dat";
2159 ex105e << tablepath <<
"ATLAS-CONF-2016-104_c_e.dat";
2161 ex106 << tablepath <<
"CMS-PAS-HIG-16-031.dat";
2163 ex106e << tablepath <<
"CMS-PAS-HIG-16-031_e.dat";
2166 ex107 << tablepath <<
"CMS-PAS-HIG-17-002.dat";
2168 ex107e << tablepath <<
"CMS-PAS-HIG-17-002_e.dat";
2170 ex108 << tablepath <<
"170804188.dat";
2172 ex108e << tablepath <<
"170804188_e.dat";
2174 ex109 << tablepath <<
"171203143.dat";
2176 ex109e << tablepath <<
"171203143_e.dat";
2178 ex110 << tablepath <<
"171004960.dat";
2180 ex110e << tablepath <<
"171004960_e.dat";
2183 ex111 << tablepath <<
"171007235_b.dat";
2185 ex111e << tablepath <<
"171007235_b_e.dat";
2187 ex112 << tablepath <<
"170800212.dat";
2189 ex112e << tablepath <<
"170800212_e.dat";
2191 ex113 << tablepath <<
"171206386_a.dat";
2193 ex113e << tablepath <<
"171206386_a_e.dat";
2195 ex114 << tablepath <<
"171206386_b.dat";
2197 ex114e << tablepath <<
"171206386_b_e.dat";
2199 ex115 << tablepath <<
"170809638_a.dat";
2201 ex115e << tablepath <<
"170809638_a_e.dat";
2203 ex116 << tablepath <<
"170809638_b.dat";
2205 ex116e << tablepath <<
"170809638_b_e.dat";
2207 ex117 << tablepath <<
"170804445.dat";
2209 ex117e << tablepath <<
"170804445_e.dat";
2211 ex118 << tablepath <<
"CMS-PAS-B2G-16-023.dat";
2213 ex118e << tablepath <<
"CMS-PAS-B2G-16-023_e.dat";
2215 ex119 << tablepath <<
"CMS-PAS-HIG-16-001_a.dat";
2217 ex119e << tablepath <<
"CMS-PAS-HIG-16-001_a_e.dat";
2219 ex120 << tablepath <<
"CMS-PAS-HIG-16-001_b.dat";
2221 ex120e << tablepath <<
"CMS-PAS-HIG-16-001_b_e.dat";
2224 ex121 << tablepath <<
"160302991_a.dat";
2226 ex122 << tablepath <<
"160302991_b.dat";
2228 ex123 << tablepath <<
"160302991_c.dat";
2230 ex124 << tablepath <<
"160302991_d.dat";
2233 bsg1 << tablepath <<
"bsgammatable.dat";
2241 double params[] = {mass};
2257 double params[] = {mass};
2273 double params[] = {mass};
2289 double params[] = {mass};
2305 double params[] = {mass};
2321 double params[] = {mass};
2337 double params[] = {mass};
2353 double params[] = {mass};
2369 double params[] = {mass};
2375 double newResult = 0.0;
2376 if (mass>=20. && mass <=2000.) {
2388 double params[] = {mass};
2394 double newResult = 0.0;
2395 if (mass>=20. && mass <=2000.) {
2407 double params[] = {mass};
2413 double newResult = 0.0;
2414 if (mass>=20. && mass <=2000.) {
2426 double params[] = {mass};
2432 double newResult = 0.0;
2433 if (mass>=20. && mass <=2000.) {
2445 double params[] = {mass};
2451 double newResult = 0.0;
2452 if (mass>=20. && mass <=2000.) {
2464 double params[] = {mass};
2470 double newResult = 0.0;
2471 if (mass>=20. && mass <=2000.) {
2483 double params[] = {mass};
2489 double newResult = 0.0;
2490 if (mass>=20. && mass <=2000.) {
2502 double params[] = {mass};
2508 double newResult = 0.0;
2509 if (mass>=20. && mass <=2000.) {
2521 double params[] = {mass};
2527 double newResult = 0.0;
2528 if (mass>=20. && mass <=2000.) {
2540 double params[] = {mass};
2546 double newResult = 0.0;
2547 if (mass>=20. && mass <=2000.) {
2559 double params[] = {mass};
2565 double newResult = 0.0;
2566 if (mass>=20. && mass <=2000.) {
2578 double params[] = {mass};
2584 double newResult = 0.0;
2585 if (mass>=20. && mass <=2000.) {
2597 double params[] = {mass};
2603 double newResult = 0.0;
2604 if (mass>=20. && mass <=2000.) {
2616 double params[] = {mass};
2622 double newResult = 0.0;
2623 if (mass>=20. && mass <=2000.) {
2635 double params[] = {mass};
2641 double newResult = 0.0;
2642 if (mass>=20. && mass <=2000.) {
2654 double params[] = {mass};
2660 double newResult = 0.0;
2661 if (mass>=20. && mass <=2000.) {
2673 double params[] = {mass};
2679 double newResult = 0.0;
2680 if (mass>=20. && mass <=2000.) {
2692 double params[] = {mass};
2698 double newResult = 0.0;
2699 if (mass>=20. && mass <=2000.) {
2711 double params[] = {mHp,
logtb};
2717 double newResult = 0.0;
2718 if (mHp>=180. && mHp <=1400. && logtb>=-1. &&
logtb<=1.75) {
2730 double params[] = {mHp,
logtb};
2736 double newResult = 0.0;
2737 if (mHp>=180. && mHp <=2000. && logtb>=-1. &&
logtb<=1.75) {
2749 double params[] = {mass};
2765 double params[] = {mass};
2781 double params[] = {mass};
2797 double params[] = {mass};
2813 double params[] = {mass};
2829 double params[] = {mass};
2845 double params[] = {mass};
2861 double params[] = {mass};
2877 double params[] = {mass};
2893 double params[] = {mass};
2909 double params[] = {mass};
2925 double params[] = {mass};
2941 double params[] = {mass};
2957 double params[] = {mass};
2973 double params[] = {mass};
2989 double params[] = {mass};
3005 double params[] = {mass};
3021 double params[] = {mass};
3037 double params[] = {mass};
3055 double params[] = {mass};
3071 double params[] = {mass};
3087 double params[] = {mass};
3103 double params[] = {mass};
3119 double params[] = {mass};
3135 double params[] = {mass};
3151 double params[] = {mass};
3167 double params[] = {mass};
3183 double params[] = {mass};
3199 double params[] = {mass};
3215 double params[] = {mass};
3231 double params[] = {mass};
3247 double params[] = {mass};
3263 double params[] = {mass};
3279 double params[] = {mass};
3295 double params[] = {mass};
3311 double params[] = {mass};
3327 double params[] = {mass};
3343 double params[] = {mass};
3359 double params[] = {mass};
3375 double params[] = {mass};
3391 double params[] = {mass};
3407 double params[] = {mass};
3423 double params[] = {mass};
3439 double params[] = {mass};
3455 double params[] = {mass};
3471 double params[] = {mass};
3487 double params[] = {mass};
3503 double params[] = {mass};
3519 double params[] = {mass};
3535 double params[] = {mass};
3610 double params[] = {mass};
3626 double params[] = {mass};
3642 double params[] = {mass};
3658 double params[] = {mass};
3674 double params[] = {mass};
3690 double params[] = {mass};
3706 double params[] = {mA, mH};
3722 double params[] = {mA, mH};
3738 double params[] = {mA, mH};
3754 double params[] = {mA, mH};
3770 double params[] = {mass};
3786 double params[] = {mass};
3802 double params[] = {mass};
3818 double params[] = {mass};
3834 double params[] = {mass};
3850 double params[] = {mass};
3866 double params[] = {mass};
3882 double params[] = {mass};
3898 double params[] = {mass};
3914 double params[] = {mass};
3930 double params[] = {mass};
3946 double params[] = {mass};
3962 double params[] = {mass};
3978 double params[] = {mass};
3994 double params[] = {mass};
4010 double params[] = {mass};
4026 double params[] = {mass};
4042 double params[] = {mass};
4058 double params[] = {mass};
4074 double params[] = {mass};
4090 double params[] = {mass};
4106 double params[] = {mass};
4122 double params[] = {mass};
4138 double params[] = {mass};
4154 double params[] = {mass};
4170 double params[] = {mass};
4186 double params[] = {mass};
4202 double params[] = {mass};
4218 double params[] = {mass};
4234 double params[] = {mass};
4250 double params[] = {mass};
4266 double params[] = {mass};
4282 double params[] = {mass};
4298 double params[] = {mass};
4314 double params[] = {mass};
4330 double params[] = {mass};
4346 double params[] = {mass};
4362 double params[] = {mass};
4378 double params[] = {mass};
4394 double params[] = {mass};
4410 double params[] = {mass};
4426 double params[] = {mass};
4442 double params[] = {mass};
4458 double params[] = {mass};
4474 double params[] = {mass};
4490 double params[] = {mass};
4506 double params[] = {mass};
4522 double params[] = {mass};
4538 double params[] = {mass};
4554 double params[] = {mass};
4570 double params[] = {mass};
4586 double params[] = {mass};
4602 double params[] = {mass};
4618 double params[] = {mass};
4634 double params[] = {mass};
4650 double params[] = {mass};
4666 double params[] = {mass};
4682 double params[] = {mass};
4698 double params[] = {mass};
4714 double params[] = {mass};
4730 double params[] = {mass};
4746 double params[] = {mass};
4762 double params[] = {mass};
4778 double params[] = {mass};
4794 double params[] = {mass};
4810 double params[] = {mass};
4826 double params[] = {mass};
4842 double params[] = {mass};
4858 double params[] = {mass};
4874 double params[] = {mass};
4890 double params[] = {mass};
4906 double params[] = {mass};
4922 double params[] = {mass};
4938 double params[] = {mass};
4954 double params[] = {mass};
4970 double params[] = {mass};
4986 double params[] = {mass};
5002 double params[] = {mass};
5018 double params[] = {mass};
5034 double params[] = {mass};
5050 double params[] = {mass};
5066 double params[] = {mass};
5082 double params[] = {mass};
5098 double params[] = {mass};
5114 double params[] = {mass};
5130 double params[] = {mass};
5146 double params[] = {mass};
5162 double params[] = {mass};
5178 double params[] = {mass};
5194 double params[] = {mass};
5210 double params[] = {mass};
5226 double params[] = {mass};
5242 double params[] = {mass};
5258 double params[] = {mass};
5274 double params[] = {mass};
5290 double params[] = {mass};
5306 double params[] = {mass};
5310 return(
ie1c[NumPar][i] );
5322 double params[] = {mass};
5326 return(
ie1ec[NumPar][i] );
5338 double params[] = {mass};
5342 return(
ie2c[NumPar][i] );
5354 double params[] = {mass};
5358 return(
ie2ec[NumPar][i] );
5370 double params[] = {mass};
5374 return(
ie3c[NumPar][i] );
5386 double params[] = {mass};
5390 return(
ie3ec[NumPar][i] );
5402 double params[] = {mass};
5406 return(
ie4c[NumPar][i] );
5418 double params[] = {mass};
5422 return(
ie4ec[NumPar][i] );
5434 double params[] = {mass};
5438 return(
ie5c[NumPar][i] );
5450 double params[] = {mass};
5454 return(
ie5ec[NumPar][i] );
5466 double params[] = {mass};
5470 return(
ie6c[NumPar][i] );
5482 double params[] = {mass};
5486 return(
ie6ec[NumPar][i] );
5498 double params[] = {mass};
5502 return(
ie7c[NumPar][i] );
5514 double params[] = {mass};
5518 return(
ie7ec[NumPar][i] );
5530 double params[] = {mass};
5534 return(
ie8c[NumPar][i] );
5546 double params[] = {mass};
5550 return(
ie8ec[NumPar][i] );
5562 double params[] = {mass};
5566 return(
ie9c[NumPar][i] );
5578 double params[] = {mass};
5582 return(
ie9ec[NumPar][i] );
5594 double params[] = {mass};
5598 return(
ie10c[NumPar][i] );
5610 double params[] = {mass};
5614 return(
ie10ec[NumPar][i] );
5626 double params[] = {mass};
5630 return(
ie11c[NumPar][i] );
5642 double params[] = {mass};
5646 return(
ie11ec[NumPar][i] );
5658 double params[] = {mass};
5662 return(
ie12c[NumPar][i] );
5674 double params[] = {mass};
5678 return(
ie12ec[NumPar][i] );
5690 double params[] = {mass};
5694 return(
ie13c[NumPar][i] );
5706 double params[] = {mass};
5710 return(
ie13ec[NumPar][i] );
5722 double params[] = {mass};
5726 return(
ie14c[NumPar][i] );
5738 double params[] = {mass};
5742 return(
ie14ec[NumPar][i] );
5754 double params[] = {mass};
5758 return(
ie15c[NumPar][i] );
5770 double params[] = {mass};
5774 return(
ie15ec[NumPar][i] );
5786 double params[] = {mass};
5790 return(
ie16c[NumPar][i] );
5802 double params[] = {mass};
5806 return(
ie16ec[NumPar][i] );
5818 double params[] = {mass};
5822 return(
ie17c[NumPar][i] );
5834 double params[] = {mass};
5838 return(
ie17ec[NumPar][i] );
5850 double params[] = {mass};
5854 return(
ie18c[NumPar][i] );
5866 double params[] = {mass};
5870 return(
ie18ec[NumPar][i] );
5882 double params[] = {mass};
5886 return(
ie19c[NumPar][i] );
5898 double params[] = {mass};
5902 return(
ie19ec[NumPar][i] );
5914 double params[] = {mass};
5918 return(
ie20c[NumPar][i] );
5930 double params[] = {mass};
5934 return(
ie20ec[NumPar][i] );
5946 double params[] = {mass};
5950 return(
ie21c[NumPar][i] );
5962 double params[] = {mass};
5966 return(
ie21ec[NumPar][i] );
5978 double params[] = {mass};
5982 return(
ie22c[NumPar][i] );
5994 double params[] = {mass};
5998 return(
ie22ec[NumPar][i] );
6010 double params[] = {mass};
6014 return(
ie23c[NumPar][i] );
6026 double params[] = {mass};
6030 return(
ie23ec[NumPar][i] );
6042 double params[] = {mass};
6046 return(
ie24c[NumPar][i] );
6058 double params[] = {mass};
6062 return(
ie24ec[NumPar][i] );
6074 double params[] = {mass};
6078 return(
ie25c[NumPar][i] );
6090 double params[] = {mass};
6094 return(
ie25ec[NumPar][i] );
6106 double params[] = {mass};
6110 return(
ie26c[NumPar][i] );
6122 double params[] = {mass};
6126 return(
ie26ec[NumPar][i] );
6138 double params[] = {mass};
6142 return(
ie27c[NumPar][i] );
6154 double params[] = {mass};
6158 return(
ie27ec[NumPar][i] );
6170 double params[] = {mass};
6174 return(
ie28c[NumPar][i] );
6186 double params[] = {mass};
6190 return(
ie28ec[NumPar][i] );
6202 double params[] = {mass};
6206 return(
ie29c[NumPar][i] );
6218 double params[] = {mass};
6222 return(
ie29ec[NumPar][i] );
6234 double params[] = {mass};
6238 return(
ie30c[NumPar][i] );
6250 double params[] = {mass};
6254 return(
ie30ec[NumPar][i] );
6266 double params[] = {mass};
6270 return(
ie31c[NumPar][i] );
6282 double params[] = {mass};
6286 return(
ie31ec[NumPar][i] );
6298 double params[] = {mass};
6302 return(
ie32c[NumPar][i] );
6314 double params[] = {mass};
6318 return(
ie32ec[NumPar][i] );
6330 double params[] = {mass};
6334 return(
ie33c[NumPar][i] );
6346 double params[] = {mass};
6350 return(
ie33ec[NumPar][i] );
6362 double params[] = {mass};
6366 return(
ie34c[NumPar][i] );
6378 double params[] = {mass};
6382 return(
ie34ec[NumPar][i] );
6394 double params[] = {mass};
6398 return(
ie35c[NumPar][i] );
6410 double params[] = {mass};
6414 return(
ie35ec[NumPar][i] );
6426 double params[] = {mass};
6430 return(
ie36c[NumPar][i] );
6442 double params[] = {mass};
6446 return(
ie36ec[NumPar][i] );
6458 double params[] = {mass};
6462 return(
ie37c[NumPar][i] );
6474 double params[] = {mass};
6478 return(
ie37ec[NumPar][i] );
6490 double params[] = {mass};
6494 return(
ie38c[NumPar][i] );
6506 double params[] = {mass};
6510 return(
ie38ec[NumPar][i] );
6522 double params[] = {mass};
6526 return(
ie39c[NumPar][i] );
6538 double params[] = {mass};
6542 return(
ie39ec[NumPar][i] );
6554 double params[] = {mass};
6558 return(
ie40c[NumPar][i] );
6570 double params[] = {mass};
6574 return(
ie40ec[NumPar][i] );
6586 double params[] = {mass};
6602 double params[] = {mass};
6618 double params[] = {mass};
6634 double params[] = {mass};
6650 double params[] = {mass};
6666 double params[] = {mass};
6682 double params[] = {mass};
6698 double params[] = {mass};
6714 double params[] = {mass};
6730 double params[] = {mass};
6746 double params[] = {mass};
6762 double params[] = {mass};
6778 double params[] = {mass};
6794 double params[] = {mass};
6810 double params[] = {mass};
6826 double params[] = {mass};
6842 double params[] = {
logtb, logmHp};
6858 std::ifstream INfile;
6859 std::string lineTab;
6860 INfile.open( filename.c_str() );
6862 std::cout<<
"error: in THDMcache, table doesn't exist!"<<std::endl;
6870 while(INfile.good()){
6871 while(getline(INfile, lineTab)){
6872 if( lineTab[0]==
'#' )
continue;
6874 std::istringstream streamTab(lineTab);
6876 while(streamTab >>v){
6894 int rowN=arrayTab.
size_i();
6896 double xmin = arrayTab(0,0);
6897 double xmax = arrayTab(rowN-1,0);
6898 double interval = arrayTab(1,0)-arrayTab(0,0);
6899 int Nintervals = (x-xmin)/interval;
6911 y =(arrayTab(Nintervals+1,1)-arrayTab(Nintervals,1))/(arrayTab(Nintervals+1,0)
6912 -arrayTab(Nintervals,0))*(x-arrayTab(Nintervals,0))+arrayTab(Nintervals,1);
6921 int rowN=arrayTab.
size_i();
6923 double xmin = arrayTab(0,0);
6924 double xmax = arrayTab(rowN-1,0);
6925 double ymin = arrayTab(0,1);
6926 double ymax = arrayTab(rowN-1,1);
6927 double intervalx = arrayTab(1,0)-arrayTab(0,0);
6930 while(arrayTab(i,1)-arrayTab(i-1,1)==0&&i<30000);
6931 double intervaly = arrayTab(i,1)-arrayTab(i-1,1);
6932 int Nintervalsx = (x-xmin)/intervalx;
6933 int Nintervalsy = (y-ymin)/intervaly;
6934 if(x<xmin||x>xmax||y<ymin||y>ymax){
6939 double x1=arrayTab(i*Nintervalsy+Nintervalsx,0);
6940 double x2=arrayTab(i*Nintervalsy+Nintervalsx+1,0);
6941 double y1=arrayTab(i*Nintervalsy+Nintervalsx,1);
6942 double y2=arrayTab(i*(Nintervalsy+1)+Nintervalsx,1);
6943 return (arrayTab(i*Nintervalsy+Nintervalsx,2) * (x2-x) * (y2-y)
6944 +arrayTab(i*Nintervalsy+Nintervalsx+1,2) * (x-x1) * (y2-y)
6945 +arrayTab(i*(Nintervalsy+1)+Nintervalsx,2) * (x2-x) * (y-y1)
6946 +arrayTab(i*(Nintervalsy+1)+Nintervalsx+1,2) * (x-x1) * (y-y1))
6951 double THDMcache::ghHpHm(
const double mHp2,
const double tanb,
const double m12_2,
const double bma,
const double mHl2,
const double vev)
const {
6967 double THDMcache::g_HH_HpHm(
const double mHp2,
const double mHh2,
const double tanb,
const double m12_2,
const double bma,
const double vev)
const {
6995 +TAUc*(1.0+(1.0-TAUc)*
f_func(TAUc))+TAUt*(1.0+(1.0-TAUt)*
f_func(TAUt)));
7012 +TAUt*(1.0+(1.0-TAUt)*
f_func(TAUt)));
7020 double params[] = {
mA2,
Mc,
Mt};
7046 +TAUs*(1.0+(1.0-TAUs)*
f_func(TAUs))+TAUb*(1.0+(1.0-TAUb)*
f_func(TAUb)));
7063 +TAUb*(1.0+(1.0-TAUb)*
f_func(TAUb)));
7071 double params[] = {
mA2,
Ms,
Mb};
7097 +TAUmu*(1.0+(1.0-TAUmu)*
f_func(TAUmu))
7098 +TAUtau*(1.0+(1.0-TAUtau)*
f_func(TAUtau)));
7115 TAUtau*(1.0+(1.0-TAUtau)*
f_func(TAUtau)));
7139 double params[] = {mH,
MW};
7145 double TAUw=4.0*
MW*
MW/(mH*mH);
7154 double params[] = {
mHp2, mH};
7160 double TAUhp=4.0*
mHp2/(mH*mH);
7306 -
Int2(TAUtau,LAMtau));
7326 +
Int1(TAUtau,LAMtau)-
Int2(TAUtau,LAMtau));
7353 double params[] = {mH,
cW2,
MW,
MZ};
7359 double TAUw=4.0*
MW*
MW/(mH*mH);
7363 +((1.0+2.0/TAUw)*sW2/
cW2-(5.0+2.0/TAUw))*
Int1(TAUw,LAMw));
7377 double TAUhp=4.0*
mHp2/(mH*mH);
7392 return pow(asin(
sqrt(1.0/x)),2);
7413 return tau*lambda/(tau-lambda)/2.0+tau*tau*lambda*lambda/((tau-lambda)
7414 *(tau-lambda))/2.0*(
f_func(tau)-
f_func(lambda))+tau*tau*lambda/((tau-lambda)
7421 return -tau*lambda/(tau-lambda)/2.0*(
f_func(tau)-
f_func(lambda));
7433 double params[] = {a2, b2, c2};
7440 double newResult = 0.0;
7441 double x = (a2-b2-c2)*(a2-b2-c2)-4.0*b2*c2;
7442 if(x>0) newResult =
sqrt(std::fabs(x/a2))/2.0;
7475 double BrSM_htoWW = 2.15e-1;
7476 double BrSM_htoZZ = 2.64e-2;
7477 double BrSM_htogg = 8.57e-2;
7478 double BrSM_htoZga = 1.54e-3;
7479 double BrSM_htocc = 2.91e-2;
7504 rh_QuQu=cosa*cosa/(sinb*sinb);
7505 rh_VV=sin_ba*sin_ba;
7522 double ABSgagaTHDM=0.0;
7523 double ABSgagaSM=0.0;
7535 double ABSZgaTHDM=0.0;
7536 double ABSZgaSM=0.0;
7539 rh_gg=cosa/sinb*cosa/sinb;
7541 rh_ll=cosa/sinb*cosa/sinb;
7542 I_h_F=cosa/sinb*(fermU+fermD+fermL);
7543 A_h_F = cosa/sinb*(A_h_Ux+A_h_Dx+A_h_Lx)/
sqrt(sW2*
cW2);
7546 rh_gg=-cosa/sinb*sina/cosb+(cosa/sinb+sina/cosb)
7547 *(Sigmaggh_tt8*cosa/sinb+Sigmaggh_bb8*sina/cosb)/SigmaggF8;
7549 rh_ll=sina/cosb*sina/cosb;
7550 I_h_F=cosa/sinb*fermU -sina/cosb*(fermD+fermL);
7551 A_h_F = (cosa/sinb*A_h_Ux-sina/cosb*(A_h_Dx+A_h_Lx))/
sqrt(sW2*
cW2);
7554 rh_gg=cosa/sinb*cosa/sinb;
7556 rh_ll=sina/cosb*sina/cosb;
7557 I_h_F = cosa/sinb*(fermU+fermD) -sina/cosb*fermL;
7558 A_h_F = (cosa/sinb*(A_h_Ux+A_h_Dx)-sina/cosb*A_h_Lx)/
sqrt(sW2*
cW2);
7561 rh_gg=-cosa/sinb*sina/cosb+(cosa/sinb+sina/cosb)
7562 *(Sigmaggh_tt8*cosa/sinb+Sigmaggh_bb8*sina/cosb)/SigmaggF8;
7564 rh_ll=cosa/sinb*cosa/sinb;
7565 I_h_F = cosa/sinb*(fermU+fermL) -sina/cosb*fermD;
7566 A_h_F = (cosa/sinb*(A_h_Ux+A_h_Lx)-sina/cosb*A_h_Dx)/
sqrt(sW2*
cW2);
7569 throw std::runtime_error(
"modelflag can be only any of \"type1\", \"type2\", \"typeX\" or \"typeY\"");
7572 ABSgagaTHDM=(I_h_F+I_h_W+I_h_Hp).abs2();
7573 ABSgagaSM=(fermU+fermL+fermD+I_hSM_W).abs2();
7574 rh_gaga=ABSgagaTHDM/ABSgagaSM;
7576 ABSZgaTHDM=(A_h_F+A_h_W+A_h_Hp).abs2();
7577 ABSZgaSM=((A_h_Ux+A_h_Lx+A_h_Dx)/
sqrt(sW2*
cW2)+A_hSM_W).abs2();
7578 rh_Zga=ABSZgaTHDM/ABSZgaSM;
7586 pph13 = (SigmaggF13*
rh_gg + SigmaVBFVh13*
rh_VV + Sigmatth13*
rh_QuQu + Sigmabbh13*
rh_QdQd)/(SigmaggF13 + SigmaVBFVh13 + Sigmatth13 + Sigmabbh13);
7614 double cos_2b=cosb*cosb-sinb*sinb;
7629 double rHH_QuQu=sina/sinb*sina/sinb;
7630 double rHH_QdQd=0.0;
7665 rHH_gg=sina/sinb*sina/sinb;
7666 rHH_QdQd=sina/sinb*sina/sinb;
7667 rHH_ll=sina/sinb*sina/sinb;
7668 I_HH_F=sina/sinb*(I_HH_Ux+I_HH_Dx+I_HH_Lx);
7669 A_HH_F = sina/sinb*(A_HH_Ux+A_HH_Dx+A_HH_Lx)/
sqrt(sW2*
cW2);
7672 rHH_gg=sina/sinb*cosa/cosb+rSigmaggH_t8*sina/sinb*(sina/sinb-cosa/cosb)
7673 +rSigmaggH_b8*cosa/cosb*(cosa/cosb-sina/sinb);
7674 rHH_QdQd=cosa/cosb*cosa/cosb;
7675 rHH_ll=cosa/cosb*cosa/cosb;
7676 I_HH_F=sina/sinb*I_HH_Ux+cosa/cosb*(I_HH_Dx+I_HH_Lx);
7677 A_HH_F = (sina/sinb*A_HH_Ux+cosa/cosb*(A_HH_Dx+A_HH_Lx))/
sqrt(sW2*
cW2);
7680 rHH_gg=sina/sinb*sina/sinb;
7681 rHH_QdQd=sina/sinb*sina/sinb;
7682 rHH_ll=cosa/cosb*cosa/cosb;
7683 I_HH_F=sina/sinb*(I_HH_Ux+I_HH_Dx)+cosa/cosb*I_HH_Lx;
7684 A_HH_F = (sina/sinb*(A_HH_Ux+A_HH_Dx)+cosa/cosb*A_HH_Lx)/
sqrt(sW2*
cW2);
7687 rHH_gg=sina/sinb*cosa/cosb+rSigmaggH_t8*sina/sinb*(sina/sinb-cosa/cosb)
7688 +rSigmaggH_b8*cosa/cosb*(cosa/cosb-sina/sinb);
7689 rHH_QdQd=cosa/cosb*cosa/cosb;
7690 rHH_ll=sina/sinb*sina/sinb;
7691 I_HH_F=sina/sinb*(I_HH_Ux+I_HH_Lx)+cosa/cosb*I_HH_Dx;
7692 A_HH_F = (sina/sinb*(A_HH_Ux+A_HH_Lx)+cosa/cosb*A_HH_Dx)/
sqrt(sW2*
cW2);
7695 throw std::runtime_error(
"modelflag can be only any of \"type1\", \"type2\", \"typeX\" or \"typeY\"");
7699 double Gamma_Hgaga=GF*
Ale*
Ale*mHh*mHh*mHh/(
sqrt(2.0)*128.0*M_PI*M_PI*M_PI)
7700 *(I_HH_F+I_HH_W+I_HH_Hp).abs2();
7704 *(1.0-
MZ*
MZ/(mHh*mHh))*(1.0-
MZ*
MZ/(mHh*mHh))*(1.0-
MZ*
MZ/(mHh*mHh))
7705 *(A_HH_F+A_HH_W+A_HH_Hp).abs2();
7708 double Gamma_Hgg=
rHH_gg*GF*
Als*
Als*mHh*mHh*mHh/(
sqrt(2.0)*16.0*M_PI*M_PI*M_PI)
7709 *(9.0/4.0)*(I_HH_Ux/4.0+I_HH_Dx).abs2();
7719 if (mHh>=20. && mHh <=2000.) {
7746 *std::fabs((cos_ba*cos_ba/(4.0*sinb*cosb*sinb*cosb)
7748 *sina*cosa,2))/(
vev*
vev))/(8.0*mHh*M_PI);
7761 *sin_ba*sin_ba/(2.0*M_PI*
vev*
vev);
7765 GammaHtot= ((BrSM_Htott+BrSM_Htocc)*rHH_QuQu
7766 +BrSM_Htobb*rHH_QdQd
7767 +(BrSM_Htotautau+BrSM_Htomumu)*rHH_ll
7769 +Gamma_Hgaga+Gamma_HZga+Gamma_Hgg
7770 +GammaHhh+GammaHHpHm+GammaHAA+GammaHAZ+GammaHHpW;
7809 double Br_Ztoee=0.03363;
7810 double Br_Ztomumu=0.03366;
7811 double Br_Ztotautau=0.0337;
7812 double Br_Ztoinv=0.2;
7813 double Br_Wtoenu=0.1071;
7814 double Br_Wtomunu=0.1063;
7815 double Br_Wtotaunu=0.1138;
7917 +
THDM_BR_h_ZZ*2.0*Br_Ztoinv*(Br_Ztoee+Br_Ztomumu+Br_Ztotautau*0.124));
8018 double nftos=1.95996398454;
8020 if(mHh>=65.0 && mHh<90.0)
8024 if(mA>=50.0 && mA<=1000.0)
8030 else if(mHh>=90.0 && mHh<100.0)
8042 if(mA>=50.0 && mA<=1000.0)
8053 else if(mHh>=100.0 && mHh<130.0)
8065 if(mA>=50.0 && mA<=1000.0)
8076 else if(mHh>=130.0 && mHh<140.0)
8090 if(mA>=50.0 && mA<=1000.0)
8105 else if(mHh>=140.0 && mHh<145.0)
8123 if(mA>=50.0 && mA<=1000.0)
8138 else if(mHh>=145.0 && mHh<150.0)
8158 if(mA>=50.0 && mA<=1000.0)
8173 else if(mHh>=150.0 && mHh<175.0)
8193 if(mA>=50.0 && mA<=1000.0)
8208 else if(mHh>=175.0 && mHh<200.0)
8228 if(mA>=40.0 && mA<=910.0)
8233 if(mA>=50.0 && mA<=1000.0)
8248 else if(mHh>=200.0 && mHh<220.0)
8272 if(mA>=40.0 && mA<=910.0)
8277 if(mA>=50.0 && mA<=1000.0)
8312 else if(mHh>=220.0 && mHh<250.0)
8336 if(mA>=40.0 && mA<=910.0)
8341 if(mA>=50.0 && mA<=1000.0)
8376 else if(mHh>=250.0 && mHh<260.0)
8406 if(mA>=40.0 && mA<=910.0)
8411 if(mA>=50.0 && mA<=1000.0)
8456 else if(mHh>=260.0 && mHh<270.0)
8492 if(mA>=40.0 && mA<=910.0)
8497 if(mA>=50.0 && mA<=1000.0)
8550 else if(mHh>=270.0 && mHh<275.0)
8586 if(mA>=40.0 && mA<=910.0)
8591 if(mA>=50.0 && mA<=1000.0)
8644 else if(mHh>=275.0 && mHh<300.0)
8682 if(mA>=40.0 && mA<=910.0)
8687 if(mA>=50.0 && mA<=1000.0)
8742 else if(mHh>=300.0 && mHh<350.0)
8786 if(mA>=40.0 && mA<=910.0)
8791 if(mA>=50.0 && mA<=1000.0)
8868 else if(mHh>=350.0 && mHh<400.0)
8910 if(mA>=40.0 && mA<=910.0)
8915 if(mA>=50.0 && mA<=1000.0)
8994 else if(mHh>=400.0 && mHh<500.0)
9036 if(mA>=40.0 && mA<=910.0)
9041 if(mA>=50.0 && mA<=1000.0)
9122 else if(mHh>=500.0 && mHh<550.0)
9164 if(mA>=40.0 && mA<=910.0)
9169 if(mA>=50.0 && mA<=1000.0)
9254 else if(mHh>=550.0 && mHh<600.0)
9298 if(mA>=40.0 && mA<=910.0)
9303 if(mA>=50.0 && mA<=1000.0)
9390 else if(mHh>=600.0 && mHh<650.0)
9430 if(mA>=40.0 && mA<=910.0)
9435 if(mA>=50.0 && mA<=1000.0)
9520 else if(mHh>=650.0 && mHh<760.0)
9562 if(mA>=40.0 && mA<=910.0)
9567 if(mA>=50.0 && mA<=1000.0)
9654 else if(mHh>=760.0 && mHh<850.0)
9696 if(mA>=40.0 && mA<=910.0)
9701 if(mA>=50.0 && mA<=1000.0)
9788 else if(mHh>=850.0 && mHh<900.0)
9824 if(mA>=40.0 && mA<=910.0)
9829 if(mA>=50.0 && mA<=1000.0)
9916 else if(mHh>=900.0 && mHh<1000.0)
9950 if(mA>=40.0 && mA<=910.0)
9955 if(mA>=50.0 && mA<=1000.0)
10032 else if(mHh>=1000.0 && mHh<1100.0)
10108 else if(mHh>=1100.0 && mHh<1200.0)
10180 else if(mHh>=1200.0 && mHh<1500.0)
10246 else if(mHh>=1500.0 && mHh<1600.0)
10308 else if(mHh>=1600.0 && mHh<2000.0)
10368 else if(mHh>=2000.0 && mHh<2250.0)
10424 else if(mHh>=2250.0 && mHh<2400.0)
10476 else if(mHh>=2400.0 && mHh<2500.0)
10526 else if(mHh>=2500.0 && mHh<2530.0)
10572 else if(mHh>=2530.0 && mHh<2700.0)
10614 else if(mHh>=2700.0 && mHh<3000.0)
10654 else if(mHh>=3000.0 && mHh<3200.0)
10667 else if(mHh>=3200.0 && mHh<4000.0)
10697 double rA_QuQu=(cosb*cosb)/(sinb*sinb);
10698 double rA_QdQd=0.0;
10724 rA_gg=-cosb/sinb*cosb/sinb+2.0*cosb/sinb*cosb/sinb*(rSigmaggA_t8+rSigmaggA_b8);
10725 rA_QdQd=cosb/sinb*cosb/sinb;
10726 rA_ll=cosb/sinb*cosb/sinb;
10727 I_A_F=cosb/sinb*(I_A_Ux-I_A_Dx-I_A_Lx);
10728 A_A_F=cosb/sinb*(A_A_Ux-A_A_Dx-A_A_Lx);
10731 rA_gg= 1.0+(cosb/sinb-sinb/cosb)*(rSigmaggA_t8*cosb/sinb-rSigmaggA_b8*sinb/cosb);
10732 rA_QdQd=sinb/cosb*sinb/cosb;
10733 rA_ll=sinb/cosb*sinb/cosb;
10734 I_A_F=cosb/sinb*I_A_Ux+sinb/cosb*(I_A_Dx+I_A_Lx);
10735 A_A_F=cosb/sinb*A_A_Ux+sinb/cosb*(A_A_Dx+A_A_Lx);
10738 rA_gg=-cosb/sinb*cosb/sinb+2.0*cosb/sinb*cosb/sinb*(rSigmaggA_t8+rSigmaggA_b8);
10739 rA_QdQd=cosb/sinb*cosb/sinb;
10740 rA_ll=sinb/cosb*sinb/cosb;
10741 I_A_F=cosb/sinb*(I_A_Ux-I_A_Dx)+sinb/cosb*I_A_Lx;
10742 A_A_F=cosb/sinb*(A_A_Ux-A_A_Dx)+sinb/cosb*A_A_Lx;
10745 rA_gg=1.0+(cosb/sinb-sinb/cosb)*(rSigmaggA_t8*cosb/sinb-rSigmaggA_b8*sinb/cosb);
10746 rA_QdQd=sinb/cosb*sinb/cosb;
10747 rA_ll=cosb/sinb*cosb/sinb;
10748 I_A_F=cosb/sinb*(I_A_Ux-I_A_Lx)+sinb/cosb*I_A_Dx;
10749 A_A_F=cosb/sinb*(A_A_Ux-A_A_Lx)+sinb/cosb*A_A_Dx;
10752 throw std::runtime_error(
"modelflag can be only any of \"type1\", \"type2\", \"typeX\" or \"typeY\"");
10756 double Gamma_Agaga=GF*
Ale*
Ale*mA*mA*mA/(
sqrt(2.0)*128.0*M_PI*M_PI*M_PI)
10760 *(1.0-
MZ*
MZ/(mA*mA))*(1.0-
MZ*
MZ/(mA*mA))*(1.0-
MZ*
MZ/(mA*mA))
10763 double Gamma_Agg=
rA_gg*GF*
Als*
Als*mA*mA*mA/(
sqrt(2.0)*16.0*M_PI*M_PI*M_PI)
10764 *(9.0/4.0)*(I_A_Ux/4.0+I_A_Dx).abs2();
10787 *sin_ba*sin_ba/(2.0*M_PI*
vev*
vev);
10790 *cos_ba*cos_ba/(2.0*M_PI*
vev*
vev);
10795 GammaAtot= ((BrSM_Atott+BrSM_Atocc)*rA_QuQu
10796 +BrSM_Atobb*rA_QdQd
10797 +(BrSM_Atotautau+BrSM_Atomumu)*rA_ll)*GammaAtotSM
10798 +Gamma_Agaga+
Gamma_AZga+Gamma_Agg+GammaAHZ+GammaAhZ+GammaAHpW;
10814 double Br_Ztoee=0.03363;
10815 double Br_Ztomumu=0.03366;
10916 double nftos=1.95996398454;
10918 if(mA>=65.0 && mA<90.0)
10922 if(mHh>=50.0 && mHh<=1000.0)
10928 else if(mA>=90.0 && mA<100.0)
10940 if(mHh>=50.0 && mHh<=1000.0)
10951 else if(mA>=100.0 && mA<150.0)
10965 if(mHh>=50.0 && mHh<=1000.0)
10976 else if(mA>=150.0 && mA<175.0)
10992 if(mHh>=50.0 && mHh<=1000.0)
11003 else if(mA>=175.0 && mA<200.0)
11019 if(mHh>=50.0 && mHh<=910.0)
11024 if(mHh>=50.0 && mHh<=1000.0)
11035 else if(mA>=200.0 && mA<220.0)
11055 if(mHh>=50.0 && mHh<=910.0)
11060 if(mHh>=50.0 && mHh<=1000.0)
11081 else if(mA>=220.0 && mA<225.0)
11107 if(mHh>=50.0 && mHh<=910.0)
11112 if(mHh>=50.0 && mHh<=1000.0)
11133 else if(mA>=225.0 && mA<250.0)
11161 if(mHh>=50.0 && mHh<=910.0)
11166 if(mHh>=50.0 && mHh<=1000.0)
11187 else if(mA>=250.0 && mA<300.0)
11215 if(mHh>=50.0 && mHh<=910.0)
11220 if(mHh>=50.0 && mHh<=1000.0)
11245 else if(mA>=300.0 && mA<350.0)
11273 if(mHh>=50.0 && mHh<=910.0)
11278 if(mHh>=50.0 && mHh<=1000.0)
11305 else if(mA>=350.0 && mA<400.0)
11331 if(mHh>=50.0 && mHh<=910.0)
11336 if(mHh>=50.0 && mHh<=1000.0)
11365 else if(mA>=400.0 && mA<500.0)
11393 if(mHh>=50.0 && mHh<=910.0)
11398 if(mHh>=50.0 && mHh<=1000.0)
11431 else if(mA>=500.0 && mA<550.0)
11459 if(mHh>=50.0 && mHh<=910.0)
11464 if(mHh>=50.0 && mHh<=1000.0)
11499 else if(mA>=550.0 && mA<600.0)
11527 if(mHh>=50.0 && mHh<=910.0)
11532 if(mHh>=50.0 && mHh<=1000.0)
11569 else if(mA>=600.0 && mA<650.0)
11593 if(mHh>=50.0 && mHh<=910.0)
11598 if(mHh>=50.0 && mHh<=1000.0)
11635 else if(mA>=650.0 && mA<850.0)
11659 if(mHh>=50.0 && mHh<=910.0)
11664 if(mHh>=50.0 && mHh<=1000.0)
11703 else if(mA>=850.0 && mA<900.0)
11725 if(mHh>=50.0 && mHh<=910.0)
11730 if(mHh>=50.0 && mHh<=1000.0)
11769 else if(mA>=900.0 && mA<1000.0)
11789 if(mHh>=50.0 && mHh<=910.0)
11794 if(mHh>=50.0 && mHh<=1000.0)
11833 else if(mA>=1000.0 && mA<1200.0)
11871 else if(mA>=1200.0 && mA<1600.0)
11905 else if(mA>=1600.0 && mA<2000.0)
11937 else if(mA>=2000.0 && mA<2250.0)
11963 else if(mA>=2250.0 && mA<2400.0)
11985 else if(mA>=2400.0 && mA<2500.0)
12005 else if(mA>=2500.0 && mA<2700.0)
12023 else if(mA>=2700.0 && mA<3000.0)
12039 else if(mA>=3000.0 && mA<3200.0)
12050 else if(mA>=3200.0 && mA<4000.0)
12069 double gHpt=-1.0/
tanb;
12074 double SigmaHp8=0.0;
12075 double SigmaHp13=0.0;
12108 throw std::runtime_error(
"modelflag can be only any of \"type1\", \"type2\", \"typeX\" or \"typeY\"");
12111 double GammaHptaunu=
HSTheta(mHp-
Mtau)*(Mtau2*(
mHp2-Mtau2)*(
mHp2-Mtau2)*gHptau*gHptau)/(8.0*mHp*
mHp2*M_PI*vev2);
12112 double GammaHptb=
HSTheta(mHp-
Mt-
Mb)*(Vtb*Vtb/(8.0*mHp*M_PI*vev2))*3.0*(-4.0*gHpb*gHpt*Mb2*Mt2
12113 -gHpb*gHpb*Mb2*(Mb2-
mHp2+Mt2)-gHpt*gHpt*Mt2*(Mb2-
mHp2+Mt2))
12120 *MW2*MW2/mHp/(2.0*M_PI*vev2);
12121 GammaHptot= GammaHptaunu + GammaHptb + GammaHpHlW + GammaHpHhW + GammaHpAW;
12123 double Br_Hptotaunu=GammaHptaunu/
GammaHptot;
12157 double nftos=1.95996398454;
12159 if(mHp>=180.0 && mHp<200.0)
12171 else if(mHp>=200.0 && mHp<300.0)
12191 else if(mHp>=300.0 && mHp<600.0)
12213 else if(mHp>=600.0 && mHp<1000.0)
12229 else if(mHp>=1000.0 && mHp<1400.0)
12240 else if(mHp>=1400.0 && mHp<2000.0)
12263 if( RGEorder ==
"LO" ) flag=0;
12264 else if( RGEorder ==
"approxNLO" ) flag=1;
12265 else if( RGEorder ==
"NLO" ) flag=2;
12267 throw std::runtime_error(
"RGEorder can be only any of \"LO\", \"approxNLO\" or \"NLO\"");
12272 double g3_at_MZ=
sqrt(4.0*M_PI*
Als);
12275 double Ybottom1_at_MZ=0.0;
12276 double Ybottom2_at_MZ=0.0;
12277 double Ytau1_at_MZ=0.0;
12278 double Ytau2_at_MZ=0.0;
12281 double Mb_at_MZ=2.96;
12282 double Mtau_at_MZ=1.75;
12285 Ybottom2_at_MZ=(
sqrt(2.0)*Mb_at_MZ)/(
vev*sinb);
12286 Ytau2_at_MZ=(
sqrt(2.0)*Mtau_at_MZ)/(
vev*sinb);
12289 Ybottom1_at_MZ=(
sqrt(2.0)*Mb_at_MZ)/(
vev*cosb);
12290 Ytau1_at_MZ=(
sqrt(2.0)*Mtau_at_MZ)/(
vev*cosb);
12294 Ybottom2_at_MZ=(
sqrt(2.0)*Mb_at_MZ)/(
vev*sinb);
12295 Ytau1_at_MZ=(
sqrt(2.0)*Mtau_at_MZ)/(
vev*cosb);
12299 Ybottom1_at_MZ=(
sqrt(2.0)*Mb_at_MZ)/(
vev*cosb);
12300 Ytau2_at_MZ=(
sqrt(2.0)*Mtau_at_MZ)/(
vev*sinb);
12304 throw std::runtime_error(
"modelflag can be only any of \"type1\", \"type2\", \"typeX\" or \"typeY\"");
12342 double InitVals[14];
12343 InitVals[0]=g1_at_MZ;
12344 InitVals[1]=g2_at_MZ;
12345 InitVals[2]=g3_at_MZ;
12346 InitVals[3]=Ytop_at_MZ;
12347 InitVals[4]=Ybottom1_at_MZ+Ybottom2_at_MZ;
12348 InitVals[5]=Ytau1_at_MZ+Ytau2_at_MZ;
12349 InitVals[6]=m11_2_at_MZ;
12350 InitVals[7]=m22_2_at_MZ;
12351 InitVals[8]=m12_2_at_MZ;
12352 InitVals[9]=lambda1_at_MZ;
12353 InitVals[10]=lambda2_at_MZ;
12354 InitVals[11]=lambda3_at_MZ;
12355 InitVals[12]=lambda4_at_MZ;
12356 InitVals[13]=lambda5_at_MZ;
12385 throw std::runtime_error(
"modelflag can be only any of \"type1\", \"type2\", \"typeX\" or \"typeY\"");
12399 std::cout <<
"\n \t Warning: UnitarityTestflag is switched on.\n \t All lambda's are hard coded in this mode!\n" << std::endl;
12422 double WFRcomb1a = 0.0;
12423 double WFRcomb1b = 0.0;
12424 double WFRcomb2a = 0.0;
12425 double WFRcomb3a = 0.0;
12426 double WFRcomb3b = 0.0;
12427 double WFRcomb4a = 0.0;
12433 double beta=atan(
tanb);
12434 double alpha=beta-
bma;
12517 +
cos(alpha+beta)*(-8.0*
m12_2+(3.0*
mHl2+2.0*
mHp2)*2.0*sinb*cosb),2))/(64.0*
pow(cosb*sinb,4)) * ddpB0mhmHpmHp
12518 + (
sin(alpha)*
sin(alpha)*
pow((2.0*
mA2-
mHl2)*
cos(alpha-3.0*beta)*2.0*sinb*cosb+
cos(alpha+beta)*(8.0*
m12_2-(2.0*
mA2+3.0*
mHl2)*2.0*sinb*cosb),2))/(128.0*
pow(cosb*sinb,4)) * ddpB0mhmAmA
12538 +
pow((2.0*
mA2-
mHl2)*
cos(alpha-3.0*beta)*2.0*sinb*cosb+
cos(alpha+beta)*(8.0*
m12_2-(2.0*
mA2+3.0*
mHl2)*2.0*sinb*cosb),2)/(64.0*
pow(cosb,4)*sinb*sinb) * ddpB0mAmAmh
12540 + 3.0*
mHh2*
cos(alpha+3.0*beta)+16.0*
m12_2*
sin(alpha+beta),2)/(256.0*
pow(cosb,4)*sinb*sinb) * ddpB0mAmAmH;
12545 +
cos(alpha+beta)*(-8.0*
m12_2+(3.0*
mHl2+2.0*
mHp2)*2.0*sinb*cosb)))/(2.0*
mHp2*sinb*cosb) * B00mHpmh
12547 -2.0*
mHp2*
cos(alpha+3.0*beta)-16.0*
m12_2*
sin(alpha+beta)))/(4.0*
mHp2*sinb*cosb) * B00mHpmH
12549 +
cos(alpha+beta)*(-8.0*
m12_2+(2.0*
mA2+3.0*
mHl2)*2.0*sinb*cosb)))/(4.0*
mA2*sinb*cosb) * B00mAmh
12551 +3.0*
mHh2*
cos(alpha+3.0*beta)+16.0*
m12_2*
sin(alpha+beta)))/(8.0*
mA2*sinb*cosb) * B00mAmH
12561 +(9.0*
m12_2*
m12_2)/(sinb*sinb*cosb*cosb))*
sin(2.0*alpha)*
sin(2.0*alpha)))/(8.0*(
mHh2-
mHl2)*sinb*sinb*cosb*cosb) * B0mhmHmh
12581 /(8.0*(
mHh2-
mHl2)*sinb*sinb*cosb*cosb) * B0mHmHmh
12594 +
cos(alpha+beta)*(-8.0*
m12_2+(3.0*
mHl2+2.0*
mHp2)*2.0*sinb*cosb)))/(2.0*
mHp2*sinb*cosb) * B0mHpmHpmh
12596 -2.0*
mHp2*
cos(alpha+3.0*beta)-16.0*
m12_2*
sin(alpha+beta)))/(4.0*
mHp2*sinb*cosb) * B0mHpmHpmH
12600 +
cos(alpha+beta)*(-8.0*
m12_2+(2.0*
mA2+3.0*
mHl2)*2.0*sinb*cosb)))/(4.0*
mA2*sinb*cosb) * B0mAmAmh
12602 +3.0*
mHh2*
cos(alpha+3.0*beta)+16.0*
m12_2*
sin(alpha+beta)))/(8.0*
mA2*sinb*cosb) * B0mAmAmH;
12615 -
mHl2*(3.0*
sin(
bma)+
sin(3.0*
bma)+
sin(3.0*alpha+beta)+3.0*
sin(alpha+3.0*beta)),2)/(1024.0*
pow(cosb*sinb,4)) * ddpB0mhmhmh
12616 + 0.5*
pow(
cos(alpha)/sinb +
sin(alpha)/cosb,2)
12619 *
sin(
bma)*
sin(
bma))/(16.0*cosb*cosb*sinb*sinb) * ddpB0mhmHmH
12621 +
cos(alpha+beta)*(-4.0*
m12_2+(3.0*
mHl2+2.0*
mHp2)*cosb*sinb),2)/(32.0*
pow(cosb*sinb,4)) * ddpB0mhmHpmHp
12623 +
cos(alpha+beta)*(4.0*
m12_2-(2.0*
mA2+3.0*
mHl2)*cosb*sinb),2)/(64.0*
pow(cosb*sinb,4)) * ddpB0mhmAmA
12627 + 0.25*
pow(
cos(alpha)/sinb +
sin(alpha)/cosb,2)
12630 *
sin(
bma)*
sin(
bma))/(8.0*cosb*cosb*sinb*sinb) * ddpB0mHmHmh
12634 -16.0*
m12_2*
sin(alpha+beta),2)/(512*
pow(cosb*sinb,4)) * ddpB0mHmHpmHp
12636 +16.0*
m12_2*
sin(alpha+beta),2)/(1024*
pow(cosb*sinb,4)) * ddpB0mHmAmA
12643 -(
m12_2*
sin(alpha+beta))/(cosb*cosb*sinb*sinb),2) * ddpB0mHpmHpmH
12648 +
cos(alpha+beta)*(8.0*
m12_2-(2.0*
mA2+3.0*
mHl2)*2.0*cosb*sinb),2)/(128.0*
pow(cosb*sinb,4)) * ddpB0mAmAmh
12650 +16.0*
m12_2*
sin(alpha+beta),2)/(512.0*
pow(cosb*sinb,4)) * ddpB0mAmAmH;
12667 +
cos(alpha+beta)*(-8.0*
m12_2+(3.0*
mHl2+2.0*
mHp2)*2.0*sinb*cosb),2))/(64.0*
pow(cosb*sinb,4)) * ddpB0mhmHpmHp
12668 + (
sin(alpha)*
sin(alpha)*
pow((2.0*
mA2-
mHl2)*
cos(alpha-3.0*beta)*2.0*sinb*cosb+
cos(alpha+beta)*(8.0*
m12_2-(2.0*
mA2+3.0*
mHl2)*2.0*sinb*cosb),2))/(128.0*
pow(cosb*sinb,4)) * ddpB0mhmAmA
12688 +
pow((2.0*
mA2-
mHl2)*
cos(alpha-3.0*beta)*2.0*sinb*cosb+
cos(alpha+beta)*(8.0*
m12_2-(2.0*
mA2+3.0*
mHl2)*2.0*sinb*cosb),2)/(64.0*
pow(cosb,4)*sinb*sinb) * ddpB0mAmAmh
12690 + 3.0*
mHh2*
cos(alpha+3.0*beta)+16.0*
m12_2*
sin(alpha+beta),2)/(256.0*
pow(cosb,4)*sinb*sinb) * ddpB0mAmAmH;
12695 +
cos(alpha+beta)*(-8.0*
m12_2+(3.0*
mHl2+2.0*
mHp2)*2.0*sinb*cosb)))/(2.0*
mHp2*sinb*cosb) * B00mHpmh
12697 -2.0*
mHp2*
cos(alpha+3.0*beta)-16.0*
m12_2*
sin(alpha+beta)))/(4.0*
mHp2*sinb*cosb) * B00mHpmH
12699 +
cos(alpha+beta)*(-8.0*
m12_2+(2.0*
mA2+3.0*
mHl2)*2.0*sinb*cosb)))/(4.0*
mA2*sinb*cosb) * B00mAmh
12701 +3.0*
mHh2*
cos(alpha+3.0*beta)+16.0*
m12_2*
sin(alpha+beta)))/(8.0*
mA2*sinb*cosb) * B00mAmH
12711 +(9.0*
m12_2*
m12_2)/(sinb*sinb*cosb*cosb))*
sin(2.0*alpha)*
sin(2.0*alpha)))/(8.0*(
mHh2-
mHl2)*sinb*sinb*cosb*cosb) * B0mhmHmh
12731 /(8.0*(
mHh2-
mHl2)*sinb*sinb*cosb*cosb) * B0mHmHmh
12744 +
cos(alpha+beta)*(-8.0*
m12_2+(3.0*
mHl2+2.0*
mHp2)*2.0*sinb*cosb))/(2.0*
mHp2*sinb*cosb) * B0mHpmHpmh
12750 +
cos(alpha+beta)*(-8.0*
m12_2+(2.0*
mA2+3.0*
mHl2)*2.0*sinb*cosb)))/(4.0*
mA2*sinb*cosb) * B0mAmAmh
12752 +3.0*
mHh2*
cos(alpha+3.0*beta)+16.0*
m12_2*
sin(alpha+beta)))/(8.0*
mA2*sinb*cosb) * B0mAmAmH;
12762 -
mHl2*(3.0*
sin(
bma)+
sin(3.0*
bma)+
sin(3.0*alpha+beta)+3.0*
sin(alpha+3.0*beta)),2)/(1024.0*
pow(cosb*sinb,4)) * ddpB0mhmhmh
12763 + 0.5*
pow(
cos(alpha)/sinb +
sin(alpha)/cosb,2)
12766 *
sin(
bma)*
sin(
bma))/(16.0*cosb*cosb*sinb*sinb) * ddpB0mhmHmH
12768 +
cos(alpha+beta)*(-4.0*
m12_2+(3.0*
mHl2+2.0*
mHp2)*cosb*sinb),2)/(32.0*
pow(cosb*sinb,4)) * ddpB0mhmHpmHp
12770 +
cos(alpha+beta)*(4.0*
m12_2-(2.0*
mA2+3.0*
mHl2)*cosb*sinb),2)/(64.0*
pow(cosb*sinb,4)) * ddpB0mhmAmA
12774 + 0.25*
pow(
cos(alpha)/sinb +
sin(alpha)/cosb,2)
12777 *
sin(
bma)*
sin(
bma))/(8.0*cosb*cosb*sinb*sinb) * ddpB0mHmHmh
12781 -16.0*
m12_2*
sin(alpha+beta),2)/(512*
pow(cosb*sinb,4)) * ddpB0mHmHpmHp
12783 +16.0*
m12_2*
sin(alpha+beta),2)/(1024*
pow(cosb*sinb,4)) * ddpB0mHmAmA
12788 +
cos(alpha+beta)*(8.0*
m12_2-(2.0*
mA2+3.0*
mHl2)*2.0*cosb*sinb),2)/(128.0*
pow(cosb*sinb,4)) * ddpB0mAmAmh
12790 +16.0*
m12_2*
sin(alpha+beta),2)/(512.0*
pow(cosb*sinb,4)) * ddpB0mAmAmH;
gslpp::complex B0p_MZ2_mHl2_0_0(const double MZ2, const double mHl2) const
double ip_ex_pp_Hp_taunu_CMS8_cache[2][CacheSize]
double THoEX_ttF_A_tt_ATLAS13
double R_pp_H_hh_bbVV_CMS13
double ip_csr_ggA_b_13_cache[2][CacheSize]
gslpp::matrix< double > CMS13_pp_H_hh_bbbb_e
complex cos(const complex &z)
gslpp::matrix< double > CMS13_bb_phi_tautau_e
double R_bbF_H_tt_ATLAS13
gslpp::complex I_A_U_cache[4][CacheSize]
double THDM_BR_h_cc
branching ratio to two quarks in the THDM.
double R_pp_H_AZ_bbll_CMS8
double ie10ec[2][CacheSize]
gslpp::matrix< double > ATLAS13_pp_Hpm_taunu_e
ATLAS expected upper cross section limits at 13 TeV, depending on the charged Higgs mass...
static const int CacheSize
Cache size.
double ie13ec[2][CacheSize]
double THoEX_ggF_A_hZ_bbZ_ATLAS13
double ie16e(double mass)
double ip_ex_pp_phi_Zga_llga_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a scalar resonance decaying to a boson ...
double ip_cs_VBFtoH_13(double mass)
Interpolating function for the H production cross section via vector boson fusion at 13 TeV...
gslpp::matrix< double > CMS13_ggF_H_hh_bbbb
double ip_ex_pp_Hpm_taunu_CMS13_cache[2][CacheSize]
double ip_ex_pp_H_hh_gagaWW_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
gslpp::complex I_H_W_cache[3][CacheSize]
gslpp::matrix< double > ATLAS13_VBF_H_WW_enumumu
double ip_ex_gg_A_hZ_bbZ_ATLAS8_cache_e[2][CacheSize]
double R_ggF_H_ZZ_nunuqq_ATLAS13
gslpp::matrix< double > temp6e
gslpp::complex B0_MZ2_mHp2_0_mHh2(const double MZ2, const double mHp2, const double mHh2) const
gslpp::complex I_H_Hp_cache[3][CacheSize]
gslpp::complex B0_MZ2_mHl2_mHh2_mHh2(const double MZ2, const double mHl2, const double mHh2) const
gslpp::matrix< double > ATLAS13_tt_phi_tt
gslpp::complex B00_MZ2_MW2_mA2_mHp2(const double MZ2, const double MW2, const double mA2, const double mHp2) const
gslpp::matrix< double > log_cs_ttA_8
double pp_A_Zga_llga_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ggF_A_hZ_tautaull_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::matrix< double > ATLAS13_pp_H_hh_gagabb_e
double R_pp_H_Zga_qqga_CMS13
double ip_ex_mu_pp_H_VV_CMS8_cache_e[2][CacheSize]
double GammaHtot
Total decay width of the heavy CP-even Higgs .
gslpp::complex B0_MZ2_0_0_mHh2(const double MZ2, const double mHh2) const
Particle getQuarks(const QCD::quark q) const
A get method to access a quark as an object of the type Particle.
double ip_ex_pp_phi_Zga_qqga_CMS13_cache[2][CacheSize]
double ie13e(double mass)
gslpp::complex B0p_MZ2_0_mA2_mHl2_cache[4][CacheSize]
double ip_ex_VBF_H_ZZ_llqq_ATLAS13_cache[2][CacheSize]
double THoEX_ggF_H_ZZ_ATLAS8
double ttF_H_tt_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double VBF_H_WW_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_ex_gg_A_hZ_bbZ_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced pseudoscalar reso...
double THoEX_ggF_A_hZ_tautauZ_ATLAS8
double ie37c[2][CacheSize]
gslpp::complex B00_MZ2_0_mHh2_mA2(const double MZ2, const double mHh2, const double mA2) const
gslpp::matrix< double > temp8e
double pp_H_Zga_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_ex_ggVV_H_WW_lnulnu_CMS13_cache_e[2][CacheSize]
double ip_ex_pp_phi_Zga_llga_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a pseudoscalar resonance decaying to a Z...
gslpp::matrix< double > temp40
double R_ggF_H_ZZ_qqllnunu_ATLAS13
double THoEX_bbF_A_tautau_CMS8
gslpp::complex B00_MZ2_0_mHp2_mHp2(const double MZ2, const double mHp2) const
double ghHpHm_cache[7][CacheSize]
gslpp::complex B00_MZ2_0_MW2_mHh2_cache[4][CacheSize]
double ip_ex_gg_H_ZZ_llqq_ATLAS13_cache[2][CacheSize]
gslpp::matrix< double > ATLAS13_gg_phi_Zga_llga
double ip_ex_gg_A_hZ_bbll_CMS8_cache_e[2][CacheSize]
double mu_pp_H_VV_TH8
Signal strength for the process with $VV=WW,ZZ$ at the LHC with 8 TeV.
double MWTHDM(const double MW) const
double ip_ex_gg_H_ZZ_qqllnunu_ATLAS13_cache[2][CacheSize]
double ip_ex_gg_H_ZZ_llnunu_ATLAS13_cache_e[2][CacheSize]
gslpp::complex B0p_MZ2_mA2_0_mHl2_cache[4][CacheSize]
double ip_ex_pp_H_VV_qqqq_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a scalar resonance decaying to two or ...
gslpp::matrix< double > ATLAS13_VBF_H_WW_lnuqq
gslpp::matrix< double > ATLAS13_gg_A_Zh_Zbb_e
const double & real() const
double ip_csr_ggA_b_8_cache[2][CacheSize]
double ggF_A_tautau_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double pp_H_Zga_llga_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ie9ec[2][CacheSize]
double ip_ex_pp_H_hh_bbtautau1_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
double R_ggF_H_ZZ_llnunu_CMS13
double THoEX_pp_H_Zga_llga_CMS13
double ip_ex_mu_pp_H_VV_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two massive...
double R_ggF_A_hZ_tautaull_CMS8
gslpp::complex B00_MZ2_0_MW2_mHl2_cache[4][CacheSize]
double ip_ex_pp_H_ZZ_llll_CMS13_cache[2][CacheSize]
double THoEX_pp_A_Zga_llga_ATLAS8
gslpp::matrix< double > ATLAS13_VBF_H_ZZ_llllnunu_e
gslpp::matrix< double > temp36
double THoEX_pp_Hp_tb_ATLAS13_2
double Br_HtoAA
branching ratio to two CP-odd Higgs bosons.
gslpp::complex B0p_MZ2_mHp2_0_mA2(const double MZ2, const double mHp2, const double mA2) const
double ip_ex_VBF_H_ZZ_llll_ATLAS13_cache[2][CacheSize]
double ip_ex_gg_H_hh_bbtautau_CMS8_cache_e[2][CacheSize]
double computeSigmaggH(const double sqrt_s) const
The ggH cross section in the Standard Model.
double ggF_A_hZ_bbZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double THoEX_VBF_VH_H_ZZ_llll_CMS13
double ie12e(double mass)
gslpp::matrix< double > ATLAS8_pp_phi_Zga_llga
double ip_ex_pp_phi_gaga_ATLAS8_cache[2][CacheSize]
double ip_ex_VBF_H_ZZ_llqq_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a vector boson fusion produced scalar re...
gslpp::matrix< double > CMS8_pp_Hp_taunu_e
CMS expected upper cross section limits at 8 TeV, depending on the charged Higgs mass...
gslpp::complex B00_MZ2_0_mHl2_mA2_cache[4][CacheSize]
double ip_ex_VBF_H_ZZ_ATLAS8_cache_e[2][CacheSize]
double ip_ex_tt_phi_tt_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a tt associated scalar resonance decayin...
gslpp::complex A_h_D(const double mHl2, const double cW2, const double Md, const double Ms, const double Mb, const double MZ) const
double ie22ec[2][CacheSize]
double ip_ex_pp_Hpm_tb_ATLAS8_cache_e[2][CacheSize]
double ie28c[2][CacheSize]
gslpp::complex B0_MZ2_mHh2_mA2_mA2_cache[4][CacheSize]
gslpp::complex B0p_MZ2_0_mHp2_mHh2_cache[4][CacheSize]
double R_pp_Hpm_taunu_ATLAS8
An observable class for the quartic Higgs potential coupling .
double ip_cs_ZtoZH_13_cache[2][CacheSize]
double ip_ex_pp_phi_hh_gagabb_CMS8_cache_e[2][CacheSize]
gslpp::matrix< double > ATLAS13_pp_H_hh_gagaWW_e
double ip_ex_gg_H_ZZ_llllnunu_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double ip_ex_bb_phi_tautau_CMS8_cache[2][CacheSize]
double ie23e(double mass)
double ie32c[2][CacheSize]
double ip_ex_pp_H_hh_gagabb_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
double ggF_H_tt_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double R_pp_H_hh_bbbb_CMS13
double ip_ex_gg_H_WW_lnuqq_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
gslpp::matrix< double > CMS8_gg_H_hh_bbtautau
double THoEX_pp_A_bb_CMS13
double THoEX_pp_Hp_taunu_CMS8
gslpp::matrix< double > CMS13_ggF_phi_Zga_e
double ip_ex_gg_phi_tautau_ATLAS13_cache[2][CacheSize]
double ie40e(double mass)
gslpp::complex B00_MZ2_MZ2_MZ2_mHh2_cache[3][CacheSize]
double ip_ex_gg_phi_gaga_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a gluon-gluon produced scalar resonance de...
gslpp::matrix< double > temp30e
double R_pp_H_hh_bblnulnu_CMS13
gslpp::complex B0p_MZ2_0_mHp2_mA2_cache[4][CacheSize]
double THoEX_ggF_H_tautau_ATLAS8
gslpp::matrix< double > temp36e
double R_ttF_A_tt_ATLAS13
double R_ttF_H_tt_ATLAS13
double THoEX_VBF_H_ZZ_llqq_ATLAS13
double ip_ex_VBF_H_WW_enumunu_ATLAS13_cache[2][CacheSize]
double THoEX_ggF_A_hZ_tautaull_CMS8
double ggF_tth13
Ratio of THDM and SM cross sections for ggF and tth production of h at 13 TeV.
double ip_ex_VBF_H_WW_ATLAS8_cache[2][CacheSize]
double ip_ex_gg_H_ZZ_llqq_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double ip_ex_pp_H_hh_bbtautau_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
gslpp::complex B0_MZ2_mHl2_mHh2_mHl2_cache[4][CacheSize]
gslpp::complex A_H_W_cache[5][CacheSize]
gslpp::matrix< double > ATLAS13_bb_phi_tautau_e
gslpp::matrix< double > CMS8_pp_A_Zga_llga_e
gslpp::matrix< double > CMS8_pp_H_AZ_bbll
gslpp::matrix< double > CMS13_gg_phi_gaga
double R_pp_A_gaga_ATLAS13
double ggF_H_hh_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::complex B0_MZ2_mHl2_mHh2_mHl2(const double MZ2, const double mHl2, const double mHh2) const
gslpp::matrix< double > temp33e
Particle getLeptons(const QCD::lepton p) const
A get method to retrieve the member object of a lepton.
double ie36ec[2][CacheSize]
double ip_Br_HPtocc(double mass)
Interpolating function for the SM branching ratio to two charm quarks.
double ip_ex_pp_H_hh_CMS8_cache[2][CacheSize]
double ip_ex_VBF_H_ZZ_llll_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a vector boson fusion produced scalar re...
double THoEX_VBF_H_ZZ_llll_ATLAS13
double ie9c[2][CacheSize]
gslpp::matrix< double > log_cs_bbH_8
SM Higgs production cross section table at 8 TeV obtained with SusHi 1.5, depending on the Higgs mass...
gslpp::complex B00_MZ2_MZ2_MZ2_mHh2(const double MZ2, const double mHh2) const
gslpp::complex B0_MZ2_0_MW2_mHh2(const double MZ2, const double MW2, const double mHh2) const
int CacheCheckReal(const double cache[][CacheSize], const int NumPar, const double params[]) const
Check whether for the latest set of parameters a value is in the cache.
double ip_ex_pp_A_HZ_tautaull_CMS8(double mA, double mH)
Interpolating function for the observed CMS upper limit on a pseudoscalar resonance decaying to whic...
double ip_ex_pp_phi_Zga_llga_ATLAS13_cache[2][CacheSize]
gslpp::complex B0_MZ2_MZ2_MZ2_mHh2(const double MZ2, const double mHh2) const
double bbF_H_tt_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::complex B0p_MZ2_mHh2_mHh2_mHh2(const double MZ2, const double mHh2) const
gslpp::matrix< double > arraybsgamma
table, depending on logtb and the logarithm of the charged Higgs mass.
double rHH_VV
Squared relative coupling of to two massive vector bosons.
double ie30ec[2][CacheSize]
gslpp::complex B00_MZ2_MW2_MW2_mHh2_cache[4][CacheSize]
gslpp::complex B0_MZ2_0_mA2_mHl2_cache[4][CacheSize]
double ie20c[2][CacheSize]
double ip_ex_pp_Hpm_taunu_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a singly charged scalar resonance decayi...
double ip_ex_bb_phi_bb_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a bottom quark produced scalar resonance d...
double ip_ex_gg_phi_tautau_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
gslpp::complex B0_MZ2_mA2_0_mHl2_cache[4][CacheSize]
void CacheShift(gslpp::complex cache[][CacheSize], const int NumPar, const double params[], const gslpp::complex newResult) const
Adds a new result and its parameters into the cache.
gslpp::complex B0_MZ2_mHl2_0_mHp2(const double MZ2, const double mHl2, const double mHp2) const
double ip_Br_HPtott(double mass)
Interpolating function for the SM branching ratio to two top quarks.
double ip_ex_gg_A_hZ_tautauZ_ATLAS8_cache_e[2][CacheSize]
gslpp::matrix< double > temp10
double ip_ex_pp_Hp_tb_ATLAS13_2(double mass)
Interpolating function for the observed ATLAS upper limit on a singly charged scalar resonance decayi...
double ip_cs_WtoWH_13(double mass)
Interpolating function for the W associated H production cross section at 13 TeV. ...
double pp_H_AZ_bbll_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::matrix< double > CMS13_pp_H_ZZ_llqq_e
gslpp::complex B0_MZ2_mHh2_0_mHp2(const double MZ2, const double mHh2, const double mHp2) const
gslpp::matrix< double > temp27e
double ip_ex_pp_Hp_tb_ATLAS13_1_cache[2][CacheSize]
gslpp::matrix< double > ATLAS8_pp_Hpm_taunu
ATLAS observed upper cross section limits at 8 TeV, depending on the charged Higgs mass...
double g_HH_HpHm_cache[7][CacheSize]
double ip_ex_mu_pp_H_VV_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two massive...
double Br_HtoHpW
branching ratio to two light Higgs bosons.
double THDM_BR_h_gg
branching ratio to two gluons in the THDM.
gslpp::matrix< double > CMS8_gg_A_hZ_tautaull_e
double THoEX_pp_H_ZZ_llll_CMS13
gslpp::complex B0_MZ2_mHh2_mHp2_mHp2(const double MZ2, const double mHh2, const double mHp2) const
gslpp::complex A_HH_L_cache[6][CacheSize]
double ip_ex_gg_H_ZZ_llqq_ATLAS13_cache_e[2][CacheSize]
double ip_ex_gg_H_ZZ_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double ip_cs_ggtoH_13_cache[2][CacheSize]
gslpp::matrix< double > ATLAS13_pp_H_hh_bbbb_e
double ip_ex_tt_phi_tt_ATLAS13_cache_e[2][CacheSize]
gslpp::matrix< double > temp31
double getNLOuniscale() const
double ip_ex_ggF_phi_Zga_CMS13_cache[2][CacheSize]
A class for constructing and defining operations on real matrices.
double ip_ex_pp_Hpm_taunu_ATLAS13_cache[2][CacheSize]
double R_ggF_H_tautau_ATLAS13
gslpp::matrix< double > temp23e
gslpp::complex B0_MZ2_mHp2_0_mHh2_cache[4][CacheSize]
gslpp::complex B0p_MZ2_mHh2_mA2_mA2_cache[4][CacheSize]
gslpp::matrix< double > CMS13_pp_phi_bb
CMS observed upper cross section limits at 13 TeV, depending on the Higgs mass.
double ie7c[2][CacheSize]
double ip_ex_gg_H_WW_ATLAS8_cache_e[2][CacheSize]
double ggF_tth8
Ratio of THDM and SM cross sections for ggF and tth production of h at 8 TeV.
gslpp::matrix< double > CMS8_pp_H_hh
double THoEX_pp_H_hh_bbtautau_CMS13
gslpp::complex B0p_MZ2_mA2_mA2_mHl2_cache[4][CacheSize]
double GammaHptot
Total decay width of the charged Higgs .
double THoEX_ggF_A_Zga_CMS13
double pp_H_AZ_tautaull_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::matrix< double > temp3
std::string getModelTypeflag() const
gslpp::matrix< double > CMS8_gg_phi_tautau_e
double ip_ex_gg_phi_Zga_llga_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double R_pp_H_AZ_tautaull_CMS8
double sumModBRs
Sum of the modified branching ratios.
double THoEX_pp_H_Zga_ATLAS13
gslpp::complex B0_MZ2_mHl2_mHp2_mHp2(const double MZ2, const double mHl2, const double mHp2) const
double ip_cs_pptobbA_8(double mass)
Interpolating function for the bottom associated A production cross section at 8 TeV.
double ip_ex_pp_Hp_tb_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a singly charged scalar resonance decaying...
gslpp::matrix< double > CMS8_pp_A_HZ_tautaull
double THoEX_ggF_H_Zga_llga_ATLAS13
gslpp::matrix< double > temp33
double THoEX_pp_H_hh_bbtautau1_CMS13
double ip_ex_pp_H_hh_bbVV_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
double R_VBF_H_ZZ_qqllnunu_ATLAS13
double ie17e(double mass)
double computeThValue()
The quartic coupling .
gslpp::complex A_h_U_cache[7][CacheSize]
double R_ggF_H_WW_enumunu_ATLAS13
gslpp::complex B00_MZ2_MZ2_MZ2_mHl2(const double MZ2, const double mHl2) const
gslpp::complex B0_MZ2_mA2_0_mHl2(const double MZ2, const double mA2, const double mHl2) const
double ip_ex_pp_H_hh_gagabb_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a scalar resonance decaying to two boso...
double pp_A_HZ_tautaull_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double THoEX_pp_H_ZZ_llqq_CMS13
gslpp::matrix< double > ATLAS13_gg_H_ZZ_nunuqq_e
gslpp::matrix< double > log_cs_VBF_8
double computeThValue()
The quartic coupling .
double ip_ex_pp_H_hh_bbtautau1_CMS13_cache[2][CacheSize]
gslpp::matrix< double > ATLAS13_VBF_H_ZZ_llll_e
double ip_ex_pp_H_ZZ_llll_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
gslpp::complex B0(const double mu2, const double p2, const double m02, const double m12) const
.
double pp_Hpm_taunu_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::complex A_h_U(const double mHl2, const double cW2, const double Mu, const double Mc, const double Mt, const double MZ) const
double ggF_H_gaga_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::matrix< double > ATLAS13_pp_Hp_tb1_e
double ip_ex_VBF_H_ZZ_qqllnunu_ATLAS13_cache[2][CacheSize]
double R_ggF_H_hh_bbtautau_CMS8
double R_pp_H_hh_bbbb_ATLAS13
double ip_ex_VBF_H_ZZ_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a vector boson fusion produced scalar re...
double ip_ex_gg_A_hZ_bbll_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a gluon-gluon produced pseudoscalar resona...
gslpp::matrix< double > ATLAS13_tt_phi_tt_e
double R_ggF_H_tautau_ATLAS8
double ip_ex_bb_phi_tautau_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a bb associated scalar resonance decayin...
double ip_ex_pp_H_hh_bbtautau_CMS13_cache[2][CacheSize]
gslpp::matrix< double > csrH_top_8
Production cross section ratio tables at 8 TeV obtained with HIGLU 4.34, depending on the Higgs mass...
gslpp::matrix< double > temp23
double ie31e(double mass)
double ip_csr_ggA_b_13(double mass)
Interpolating function for the gluon-gluon fusion A cross section ratio of the bottom-loop and the to...
gslpp::matrix< double > ATLAS13_gg_H_ZZ_llll
gslpp::complex B0_MZ2_mHl2_0_mHp2_cache[4][CacheSize]
gslpp::complex I_HH_L(const double mHh2, const double Mmu, const double Mtau) const
double R_pp_H_Zga_ATLAS13
double ggF_H_hh_bbtautau_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ie35e(double mass)
double ip_ex_pp_phi_Zga_qqga_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to a boson an...
gslpp::matrix< double > ATLAS8_gg_H_ZZ
double THoEX_pp_H_Zga_qqga_CMS13
double ip_ex_gg_H_hh_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double ggF_A_gaga_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double BrSM_htobb
SM branching ratio of .
double ip_ex_gg_phi_tautau_ATLAS13_cache_e[2][CacheSize]
gslpp::matrix< double > CMS13_pp_H_hh_bbVV
double ip_Br_HPtotautau_cache[2][CacheSize]
gslpp::matrix< double > CMS8_pp_H_hh_bbbb
double ip_ex_gg_A_hZ_tautauZ_ATLAS8_cache[2][CacheSize]
gslpp::matrix< double > ATLAS8_VBF_H_WW
double ip_ex_mu_pp_H_VV_CMS8_cache[2][CacheSize]
double pp_H_hh_bbtautau_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::complex A_A_U(const double mA2, const double cW2, const double Mc, const double Mt, const double MZ) const
gslpp::complex A_HH_D_cache[6][CacheSize]
gslpp::matrix< double > ATLAS8_gg_phi_tt_e
double R_pp_H_VV_qqqq_ATLAS13
double ie31c[2][CacheSize]
double ie15e(double mass)
gslpp::complex B0p_MZ2_mHp2_0_mHl2_cache[4][CacheSize]
double R_ggF_A_hZ_bbZ_ATLAS13
double ip_ex_VBF_H_WW_lnuqq_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a VBF produced scalar resonance decaying...
double THoEX_ggF_A_tautau_CMS13
double ip_ex_pp_A_Zga_llga_CMS8_cache_e[2][CacheSize]
double ip_ex_pp_H_hh_bbVV_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
double ggF_H_ZZ_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_ex_VBF_H_ZZ_llnunu_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a VBF produced scalar resonance decaying t...
gslpp::complex A_HH_U_cache[6][CacheSize]
gslpp::matrix< double > temp20e
double ip_ex_pp_phi_gaga_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a scalar resonance decaying to two photo...
gslpp::matrix< double > CMS13_pp_H_hh_bbtautau1
double pp_H_VV_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double computeSigmaZH(const double sqrt_s) const
The ZH production cross section in the Standard Model.
double ie38ec[2][CacheSize]
double ip_cs_WtoWH_13_cache[2][CacheSize]
gslpp::matrix< double > csrA_top_13
gslpp::matrix< double > temp11e
double ip_ex_gg_H_ZZ_ATLAS8_cache[2][CacheSize]
double THoEX_pp_H_VV_qqqq_ATLAS13
complex pow(const complex &z1, const complex &z2)
double THoEX_pp_Hp_tb_ATLAS13
double ip_ex_gg_phi_gaga_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a gluon-gluon produced scalar resonance de...
double ip_ex_pp_H_hh_gagabb_ATLAS13_cache[2][CacheSize]
double ip_csr_ggA_t_8_cache[2][CacheSize]
double ggF_A_hZ_bbZ_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_ex_pp_H_hh_bbbb_ATLAS13_cache[2][CacheSize]
double ie27e(double mass)
double ip_cs_ggtoA_8_cache[2][CacheSize]
gslpp::matrix< double > ATLAS13_bb_phi_tautau
double ip_ex_pp_H_ZZ_llnunu_CMS13_cache[2][CacheSize]
double ip_ex_VBF_H_ZZ_llllnunu_ATLAS13_cache[2][CacheSize]
double THoEX_ggF_H_gaga_CMS13
gslpp::matrix< double > ATLAS13_bb_A_Zh_Zbb
double ip_ex_gg_phi_tautau_CMS8_cache_e[2][CacheSize]
double THoEX_pp_H_gaga_ATLAS13
double ie12c[2][CacheSize]
gslpp::complex B0p_MZ2_0_mHp2_mHl2(const double MZ2, const double mHp2, const double mHl2) const
gslpp::complex B0p_MZ2_mHl2_mHh2_mHh2_cache[4][CacheSize]
double ip_cs_pptottA_8_cache[2][CacheSize]
gslpp::matrix< double > log_cs_ttA_13
gslpp::matrix< double > readTable(std::string filename, int rowN, int colN)
This function reads values from a table and returns them as an array.
double THoEX_bbF_A_tautau_ATLAS8
double ggF_H_tautau_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_ex_pp_H_VV_qqqq_ATLAS13_cache_e[2][CacheSize]
double ip_cs_ZtoZH_13(double mass)
Interpolating function for the Z associated H production cross section at 13 TeV. ...
double ggF_A_tautau_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_GammaHPtotSM(double mass)
Interpolating function for the total SM Higgs decay width.
double ie6ec[2][CacheSize]
double ip_ex_pp_H_hh_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
double ie31ec[2][CacheSize]
double THoEX_ggF_H_hh_bbbb_CMS13
gslpp::complex B0p_MZ2_mHh2_0_mHp2(const double MZ2, const double mHh2, const double mHp2) const
gslpp::matrix< double > temp13e
double THoEX_ggF_H_gaga_CMS8
double ip_cs_ggtoHp_13(double mHp, double logtb)
Interpolating function for the H+ production cross section from two gluons at 13 TeV.
double ip_ex_bb_phi_tautau_CMS13_cache[2][CacheSize]
double getBDstartaunu_A() const
double ip_ex_pp_H_ZZ_llqq_CMS13_cache_e[2][CacheSize]
gslpp::matrix< double > ATLAS13_gg_H_ZZ_qqllnunu_e
double pp_A_bb_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double THoEX_ggF_H_WW_lnuqq_ATLAS13
double computeGammaHTotal() const
The Higgs total width in the Standard Model.
double ip_ex_VBF_H_ZZ_qqllnunu_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a VBF produced scalar resonance decaying...
double BrSM_htogaga
SM branching ratio of .
gslpp::complex B0_MZ2_mA2_0_mHh2_cache[4][CacheSize]
double ie25e(double mass)
gslpp::complex B00_MZ2_MZ2_mHh2_mA2_cache[4][CacheSize]
double THoEX_pp_H_AZ_tautaull_CMS8
double ggF_H_Zga_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::matrix< double > log_cs_ggHp_13
Charged Higgs production cross section table at 13 TeV from LHCHXSWGMSSMCharged, depending on the cha...
double ip_ex_pp_H_hh_bbVV_CMS13_cache_e[2][CacheSize]
double VBF_VH_H_ZZ_llll_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double Br_Htohh
branching ratio to two light Higgs bosons.
double THoEX_bbF_A_tautau_CMS13
double R_ggF_A_Zga_llga_ATLAS13
gslpp::matrix< double > CMS13_pp_H_ZZ_llqq
double THoEX_ggF_A_tautau_ATLAS8
double R_pp_Hpm_taunu_CMS13
double ip_ex_pp_H_hh_bbbb_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a spin-2 resonance decaying to two boso...
double ggF_VBF_H_WW_lnulnu_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_cs_WtoWH_8_cache[2][CacheSize]
gslpp::complex A_h_L_cache[7][CacheSize]
double R_pp_H_hh_bbtautau_CMS13
double ip_ex_bb_phi_tt_ATLAS13_cache_e[2][CacheSize]
double ip_ex_pp_phi_Zga_llga_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to a boson an...
double ip_ex_pp_H_hh_bblnulnu_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
gslpp::matrix< double > CMS8_mu_pp_H_VV
CMS observed upper signal strength limits at 8 TeV, depending on the Higgs mass. ...
double ip_ex_pp_H_hh_bblnulnu_CMS13_cache_e[2][CacheSize]
gslpp::matrix< double > csrH_bottom_13
gslpp::complex B0_MZ2_mA2_mA2_mHl2_cache[4][CacheSize]
double BrSM_htotautau
SM branching ratio of .
double THoEX_pp_A_Zga_llga_CMS8
double ip_ex_pp_Hpm_taunu_ATLAS13_cache_e[2][CacheSize]
gslpp::matrix< double > temp9
double ie38c[2][CacheSize]
double bbF_H_tautau_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ggF_H_WW_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::matrix< double > ATLAS13_VBF_H_ZZ_qqllnunu
gslpp::complex B0p_MZ2_mHl2_mHp2_mHp2(const double MZ2, const double mHl2, const double mHp2) const
double THoEX_ggF_H_tautau_ATLAS13
gslpp::matrix< double > temp19
gslpp::complex g_func(const double x) const
g function for the Int1 function
double ip_ex_bb_phi_bb_CMS8_cache_e[2][CacheSize]
double ie17c[2][CacheSize]
gslpp::matrix< double > ATLAS13_pp_phi_gaga
gslpp::matrix< double > temp12e
double ip_ex_VBF_H_ZZ_qqllnunu_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a VBF produced scalar resonance decaying...
double R_pp_H_hh_bbtautau1_CMS13
gslpp::complex getV_tb() const
A member for returning the value of the CKM element .
gslpp::complex B00_MZ2_MW2_mA2_mHp2_cache[5][CacheSize]
gslpp::matrix< double > temp34e
gslpp::complex B0p_MZ2_mHh2_mHh2_mHh2_cache[3][CacheSize]
double ie21c[2][CacheSize]
double ie18ec[2][CacheSize]
double ip_ex_gg_phi_tautau_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a gluon-gluon produced scalar resonance de...
double getAle() const
A get method to retrieve the fine-structure constant .
double R_bbF_A_hZ_bbZ_ATLAS13
gslpp::matrix< double > ATLAS8_bb_phi_tautau_e
gslpp::matrix< double > temp22
double ip_ex_gg_A_hZ_tautauZ_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced pseudoscalar reso...
gslpp::complex B0p_MZ2_mA2_0_mHl2(const double MZ2, const double mA2, const double mHl2) const
gslpp::complex Int1(const double tau, const double lambda) const
function for Z gamma coupling to h, H and A
double ip_ex_gg_H_ZZ_nunuqq_ATLAS13_cache_e[2][CacheSize]
double ip_ex_VBF_H_ZZ_llnunu_CMS13_cache[2][CacheSize]
static const complex & i()
gslpp::matrix< double > ATLAS8_VBF_H_WW_e
double ip_ex_gg_phi_gaga_CMS13_cache[2][CacheSize]
double ip_ex_pp_phi_gaga_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a scalar resonance decaying to two photo...
double ip_ex_VBF_H_WW_enumunu_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a vector boson fusion produced scalar re...
double ip_ex_pp_A_Zga_llga_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a pseudoscalar resonance decaying to a Z b...
double R_pp_H_gaga_ATLAS8
double ip_ex_bb_phi_tautau_ATLAS8_cache[2][CacheSize]
double R_bbF_H_tautau_CMS8
double THoEX_pp_Hpm_taunu_CMS13
double ip_ex_pp_H_hh_gagaWW_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double ip_ex_pp_Hpm_taunu_ATLAS8_cache_e[2][CacheSize]
gslpp::matrix< double > temp2
double ie29c[2][CacheSize]
gslpp::matrix< double > log_cs_ggH_8
SM Higgs production cross section tables at 8 TeV from the LHC Higgs Cross Section Working Group...
gslpp::complex B0p_MZ2_mHp2_0_mA2_cache[4][CacheSize]
double getBDtaunu_SM() const
double ip_ex_gg_phi_tautau_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a gluon-gluon produced scalar resonance de...
gslpp::complex I_A_L_cache[4][CacheSize]
double THoEX_bbF_A_hZ_bbZ_ATLAS13
double R_ggF_A_hZ_tautauZ_ATLAS8
gslpp::matrix< double > temp34
double bbF_H_tautau_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::matrix< double > ATLAS13_gg_H_WW_lnuqq_e
gslpp::complex I_A_D(const double mA2, const double Ms, const double Mb) const
gslpp::matrix< double > ATLAS13_gg_H_ZZ_llll_e
double ip_Br_HPtott_cache[2][CacheSize]
double THoEX_ggF_A_hZ_bbll_CMS8
double R_ggF_H_ZZ_llll_ATLAS13
double ie11c[2][CacheSize]
gslpp::matrix< double > temp13
THDM branching ratio of .
gslpp::matrix< double > ATLAS8_gg_H_hh_e
CKM getCKM() const
A get method to retrieve the member object of type CKM.
double ip_cs_pptottA_13_cache[2][CacheSize]
gslpp::matrix< double > ATLAS13_bb_phi_tt
ATLAS observed upper cross section limits at 13 TeV, depending on the Higgs mass.
gslpp::complex I_h_L_cache[5][CacheSize]
gslpp::matrix< double > CMS13_pp_H_hh_bbbb
double R_ggF_H_hh_gagaWW_ATLAS13
double R_pp_A_HZ_tautaull_CMS8
double ip_cs_pptobbA_13_cache[2][CacheSize]
gslpp::complex B00_MZ2_MZ2_mHl2_mA2_cache[4][CacheSize]
gslpp::complex B0p_MZ2_mHh2_mHh2_mHl2_cache[4][CacheSize]
gslpp::matrix< double > CMS8_bb_phi_bb_e
~THDMcache()
THDMcache destructor.
double ip_ex_gg_A_Zh_Zbb_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced pseudoscalar reso...
double THoEX_pp_Hpm_taunu_ATLAS8
gslpp::matrix< double > CMS8_gg_phi_tautau
gslpp::complex B0p_MZ2_mA2_mA2_mHh2(const double MZ2, const double mA2, const double mHh2) const
gslpp::complex B00_MZ2_0_mHl2_mA2(const double MZ2, const double mHl2, const double mA2) const
gslpp::matrix< double > temp29e
gslpp::matrix< double > ATLAS13_gg_H_ZZ_llnunu_e
double ie5ec[2][CacheSize]
double ie37ec[2][CacheSize]
gslpp::matrix< double > CMS8_bb_phi_tautau
gslpp::matrix< double > CMS13_pp_phi_Zga_llga_e
double ip_ex_pp_phi_hh_gagabb_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
double ip_cs_pptottH_8(double mass)
Interpolating function for the top associated H production cross section at 8 TeV.
double R_ggF_A_gaga_CMS13
gslpp::complex B00_MZ2_0_MW2_mHl2(const double MZ2, const double MW2, const double mHl2) const
double ggF_H_tautau_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_ex_gg_H_ZZ_qqllnunu_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
gslpp::matrix< double > csrA_bottom_8
gslpp::complex B0p_MZ2_mHl2_mHp2_mHp2_cache[4][CacheSize]
int HSTheta(const double x) const
Heaviside function.
double ip_csr_ggA_t_8(double mass)
Interpolating function for the gluon-gluon fusion A cross section ratio of the top-loop and the total...
double ip_ex_VBF_H_ZZ_llqq_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a vector boson fusion produced scalar re...
double THoEX_pp_H_hh_bblnulnu_CMS13
gslpp::matrix< double > log_cs_VBF_13
gslpp::complex I_A_D_cache[4][CacheSize]
gslpp::complex B0p_MZ2_mHp2_mHp2_mHh2(const double MZ2, const double mHp2, const double mHh2) const
gslpp::matrix< double > CMS8_pp_Hp_taunu
CMS observed upper cross section limits at 8 TeV, depending on the charged Higgs mass...
double ip_ex_tt_phi_tt_ATLAS13_cache[2][CacheSize]
An observable class for the quadratic Higgs potential coupling .
A model class for the Standard Model.
double ie28ec[2][CacheSize]
gslpp::matrix< double > CMS13_pp_H_hh_gagabb
double ip_ex_pp_H_ZZ_llqq_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
double ip_ex_ggVV_H_WW_lnulnu_CMS13_cache[2][CacheSize]
double ip_ex_VBF_H_ZZ_llllnunu_ATLAS13_cache_e[2][CacheSize]
double ie36e(double mass)
double ip_ex_pp_Hp_tb_ATLAS13_2_cache_e[2][CacheSize]
double ie26ec[2][CacheSize]
gslpp::matrix< double > temp28
double THoEX_ggF_H_ZZ_nunuqq_ATLAS13
gslpp::matrix< double > ATLAS13_gg_phi_tautau_e
double bbF_A_bb_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ghHpHm(const double mHp2, const double tanb, const double m12_2, const double bma, const double mHl2, const double vev) const
gslpp::complex B0p_MZ2_mHp2_mHp2_mHl2(const double MZ2, const double mHp2, const double mHl2) const
double THoEX_ggF_H_ZZ_llllnunu_ATLAS13
double THoEX_VBF_H_WW_lnuqq_ATLAS13
double ip_ex_pp_Hp_tb_ATLAS13_2_cache[2][CacheSize]
double ie17ec[2][CacheSize]
double ip_ex_bb_phi_tautau_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a bottom quark produced scalar resonance d...
double ip_ex_pp_H_hh_bbtautau1_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
gslpp::complex B0_MZ2_mHh2_mHh2_mHl2(const double MZ2, const double mHh2, const double mHl2) const
double THoEX_ggF_H_WW_enumunu_ATLAS13
double ip_Br_HPtoZZ_cache[2][CacheSize]
double ip_ex_pp_H_AZ_tautaull_CMS8(double mA, double mH)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to which furt...
gslpp::complex B00_MZ2_0_MZ2_mHh2(const double MZ2, const double mHh2) const
gslpp::complex I_h_U(const double mHl2, const double Mu, const double Mc, const double Mt) const
double THoEX_bbF_A_tautau_ATLAS13
double ie25ec[2][CacheSize]
gslpp::matrix< double > ATLAS13_gg_H_ZZ_qqllnunu
gslpp::matrix< double > temp21
gslpp::matrix< double > ATLAS13_VBF_H_ZZ_llqq_e
gslpp::complex B0p_MZ2_mHh2_mHl2_mHl2(const double MZ2, const double mHh2, const double mHl2) const
void assign(const size_t &i, const size_t &j, const double &a)
double ip_ex_pp_A_HZ_tautaull_CMS8_cache[3][CacheSize]
double pp_H_gaga_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
void computeAquantities()
double ggF_H_hh_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double computeSigmaVBF(const double sqrt_s) const
The VBF cross section in the Standard Model.
double ip_ex_bb_phi_bb_CMS8_cache[2][CacheSize]
double rHH_gg
Squared relative coupling of to two photons.
double computeSigmaWH(const double sqrt_s) const
The WH production cross section in the Standard Model.
gslpp::matrix< double > ATLAS13_bb_phi_tt_e
ATLAS expected upper cross section limits at 13 TeV, depending on the Higgs mass.
double R_bbF_A_tautau_ATLAS8
double pp_A_gaga_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::matrix< double > CMS13_pp_phi_Zga_llga
double ip_ex_pp_Hp_taunu_CMS8_cache_e[2][CacheSize]
double THoEX_ggF_H_Zga_CMS13
double ie15ec[2][CacheSize]
gslpp::matrix< double > temp38
gslpp::complex B0_MZ2_mHh2_mHp2_mHp2_cache[4][CacheSize]
gslpp::complex B0_MZ2_0_mHp2_mHh2(const double MZ2, const double mHp2, const double mHh2) const
double ip_ex_gg_H_hh_bbtautau_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a gluon-gluon produced scalar resonance de...
double R_bbF_A_tautau_ATLAS13
gslpp::complex B00_MZ2_0_mHh2_mHp2(const double MZ2, const double mHh2, const double mHp2) const
double ip_cs_ZtoZH_8(double mass)
Interpolating function for the Z associated H production cross section at 8 TeV.
double ip_ex_pp_H_hh_bbbb_ATLAS13_cache_e[2][CacheSize]
gslpp::complex B0_MZ2_mHh2_mHl2_mHl2(const double MZ2, const double mHh2, const double mHl2) const
double ip_csr_ggA_t_13(double mass)
Interpolating function for the gluon-gluon fusion A cross section ratio of the top-loop and the total...
gslpp::matrix< double > br_tautau
double ip_ex_ggF_H_hh_bbbb_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
double THoEX_pp_Hpm_taunu_ATLAS13
double ip_cs_pptobbA_8_cache[2][CacheSize]
double ip_ex_pp_phi_Zga_llga_ATLAS8_cache[2][CacheSize]
gslpp::complex I_H_W(const double mH, const double MW) const
gslpp::matrix< double > temp37
gslpp::matrix< double > br_ZZ
gslpp::matrix< double > temp25
gslpp::matrix< double > ATLAS8_pp_phi_gaga
ATLAS observed upper cross section limits at 8 TeV, depending on the Higgs mass. ...
double ip_ex_pp_H_AZ_bbll_CMS8_cache[3][CacheSize]
gslpp::matrix< double > ATLAS8_gg_H_ZZ_e
double ip_ex_VBF_H_ZZ_llllnunu_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a VBF produced scalar resonance decaying...
double R_VBF_H_ZZ_llnunu_CMS13
double ip_ex_pp_H_hh_bblnulnu_CMS13_cache[2][CacheSize]
double ip_ex_VBF_H_WW_enumunu_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a vector boson fusion produced scalar re...
gslpp::complex I_h_D_cache[5][CacheSize]
gslpp::matrix< double > ATLAS13_gg_H_WW_enumumu
double ie26c[2][CacheSize]
double ip_ex_pp_phi_gaga_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a scalar resonance decaying to two photo...
double ip_ex_gg_A_hZ_tautaull_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a gluon-gluon produced pseudoscalar resona...
double ip_Br_HPtobb_cache[2][CacheSize]
gslpp::matrix< double > temp18
double ie11ec[2][CacheSize]
void computeWFRcombinations()
double R_pp_A_gaga_ATLAS8
double ie14c[2][CacheSize]
gslpp::matrix< double > ATLAS8_bb_phi_tautau
double THoEX_ggF_A_hZ_bbZ_ATLAS8
double ip_ex_pp_phi_bb_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to a b quark p...
double ip_ex_pp_Hpm_taunu_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a singly charged scalar resonance decayi...
double R_VBF_VH_H_ZZ_llll_CMS13
gslpp::matrix< double > temp15e
gslpp::complex B0_MZ2_mA2_mA2_mHl2(const double MZ2, const double mA2, const double mHl2) const
gslpp::complex B0p_MZ2_mHl2_0_mHp2_cache[4][CacheSize]
double ip_ex_bb_phi_tautau_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a bottom quark produced scalar resonance...
gslpp::complex B0_MZ2_MZ2_MZ2_mHh2_cache[3][CacheSize]
double getAlsMz() const
A get method to access the value of .
gslpp::matrix< double > CMS8_bb_phi_tautau_e
double ie32e(double mass)
gslpp::complex I_HH_U(const double mHh2, const double Mc, const double Mt) const
double ie30c[2][CacheSize]
double ip_cs_VBFtoH_8(double mass)
Interpolating function for the H production cross section via vector boson fusion at 8 TeV...
double ip_ex_pp_phi_Zga_llga_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a pseudoscalar resonance decaying to a Z...
double pp_Hp_tb_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_ex_gg_H_ZZ_llll_ATLAS13_cache[2][CacheSize]
gslpp::matrix< double > temp1
data templates
double ip_cs_pptottH_13_cache[2][CacheSize]
double pp_Hp_taunu_TH8
Cross section for the process at the LHC with 8 TeV.
double ie8ec[2][CacheSize]
gslpp::matrix< double > CMS13_gg_H_ZZ_llnunu
double ip_ex_pp_H_hh_bbbb_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
double THoEX_bbF_A_bb_CMS8
double ie15c[2][CacheSize]
double ip_ex_pp_H_hh_gagabb_CMS13_cache_e[2][CacheSize]
gslpp::complex B00_MZ2_MW2_mHp2_mHp2_cache[4][CacheSize]
gslpp::complex B0p_MZ2_mA2_0_mHh2_cache[4][CacheSize]
gslpp::complex B0_MZ2_mHh2_0_0(const double MZ2, const double mHh2) const
double R_VBF_H_WW_enumunu_ATLAS13
gslpp::matrix< double > CMS13_pp_phi_Zga_qqga_e
gslpp::matrix< double > log_cs_ZH_13
gslpp::complex B0_MZ2_mHh2_mA2_mA2(const double MZ2, const double mHh2, const double mA2) const
double ip_Br_HPtotautau(double mass)
Interpolating function for the SM branching ratio to two tau leptons.
gslpp::complex B0_MZ2_mHh2_0_mA2_cache[4][CacheSize]
gslpp::matrix< double > br_tt
SM Higgs branching ratio tables (obtained with HDECAY 6.10), depending on the Higgs mass...
double THoEX_pp_A_gaga_ATLAS8
gslpp::complex B0_MZ2_mA2_0_mHh2(const double MZ2, const double mA2, const double mHh2) const
double R_pp_A_Zga_llga_CMS13
gslpp::complex A_A_L_cache[6][CacheSize]
double ip_ex_pp_H_AZ_tautaull_CMS8_cache[3][CacheSize]
gslpp::matrix< double > ATLAS13_gg_A_Zh_Zbb
double ip_ex_pp_H_VV_qqqq_ATLAS13_cache[2][CacheSize]
double THoEX_ggF_A_Zga_llga_ATLAS13
gslpp::complex B0_MZ2_0_MW2_mHl2(const double MZ2, const double MW2, const double mHl2) const
double ip_csr_ggA_t_13_cache[2][CacheSize]
double getMz() const
A get method to access the mass of the boson .
gslpp::matrix< double > ATLAS13_gg_phi_tautau
double ip_csr_ggA_b_8(double mass)
Interpolating function for the gluon-gluon fusion A cross section ratio of the bottom-loop and the to...
double R_ggF_H_Zga_llga_ATLAS13
double ip_ex_gg_H_ZZ_llnunu_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a gluon-gluon produced scalar resonance de...
double ip_ex_gg_A_hZ_tautaull_CMS8_cache_e[2][CacheSize]
double ip_ex_gg_H_hh_ATLAS8_cache_e[2][CacheSize]
double THoEX_pp_H_bb_CMS13
gslpp::complex B00_MZ2_0_MW2_mHh2(const double MZ2, const double MW2, const double mHh2) const
double ip_cs_pptobbH_13_cache[2][CacheSize]
gslpp::complex B00_MZ2_MZ2_mHh2_mA2(const double MZ2, const double mHh2, const double mA2) const
gslpp::matrix< double > CMS13_pp_H_ZZ_llnunu
gslpp::complex B0_MZ2_MW2_MW2_mHl2(const double MZ2, const double MW2, const double mHl2) const
double THoEX_bbF_A_tt_ATLAS13
complex log10(const complex &z)
gslpp::matrix< double > br_WW
gslpp::complex B0p_MZ2_0_0_mHh2(const double MZ2, const double mHh2) const
double ip_ex_bb_phi_tautau_CMS13_cache_e[2][CacheSize]
double ip_ex_VBF_H_WW_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a vector boson fusion produced scalar re...
gslpp::complex B0_MZ2_mHl2_mA2_mA2_cache[4][CacheSize]
double ip_Br_HPtoWW(double mass)
Interpolating function for the SM branching ratio to two bosons.
double ip_ex_gg_H_ZZ_llllnunu_ATLAS13_cache_e[2][CacheSize]
double THoEX_ggF_H_tautau_CMS8
gslpp::matrix< double > temp22e
double computeThValue()
The quartic coupling .
double ie39ec[2][CacheSize]
gslpp::matrix< double > CMS8_gg_A_hZ_bbll_e
CMS expected upper cross section limits at 8 TeV, depending on the Higgs mass.
double ip_ex_VBF_H_WW_ATLAS8_cache_e[2][CacheSize]
double ie8c[2][CacheSize]
double ip_ex_pp_phi_hh_bbbb_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
double ip_ex_gg_H_ZZ_ATLAS8_cache_e[2][CacheSize]
gslpp::complex B0_MZ2_mHp2_mHp2_mHl2_cache[4][CacheSize]
gslpp::complex B0_MZ2_0_0_mHl2_cache[3][CacheSize]
double getBDstartaunu_SM() const
double ip_ex_bsgamma_cache[3][CacheSize]
double R_pp_A_Zga_llga_ATLAS8
double ip_ex_gg_H_ZZ_llnunu_CMS13_cache_e[2][CacheSize]
double ip_cs_VBFtoH_8_cache[2][CacheSize]
double ie36c[2][CacheSize]
double ie4ec[2][CacheSize]
double THoEX_ggF_A_gaga_CMS8
gslpp::complex A_A_D_cache[6][CacheSize]
gslpp::matrix< double > temp38e
gslpp::matrix< double > ATLAS13_VBF_H_ZZ_llll
double ip_ex_pp_H_hh_bbtautau_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
void CacheShiftReal(double cache[][CacheSize], const int NumPar, const double params[], const double newResult) const
Adds a new result and its parameters into the cache.
gslpp::matrix< double > log_cs_WH_13
double ip_ex_ggF_phi_Zga_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a gluon-gluon produced scalar resonance de...
gslpp::matrix< double > temp16
double R_pp_A_Zga_llga_CMS8
double ie19ec[2][CacheSize]
An observable class for the quartic Higgs potential coupling .
gslpp::complex I_HH_D_cache[4][CacheSize]
double ip_ex_gg_H_WW_enumunu_ATLAS13_cache_e[2][CacheSize]
double ie23ec[2][CacheSize]
double ie16ec[2][CacheSize]
double ip_ex_VBF_H_ZZ_qqllnunu_ATLAS13_cache_e[2][CacheSize]
double ip_ex_gg_phi_gaga_CMS8_cache_e[2][CacheSize]
gslpp::matrix< double > temp16e
double ip_ex_VBF_H_ZZ_llllnunu_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a VBF produced scalar resonance decaying...
double rh_gg
Squared relative coupling of to two gluons.
gslpp::matrix< double > CMS13_ggF_phi_Zga
double ip_cs_ggtoHp_8(double mHp, double logtb)
Interpolating function for the H+ production cross section from two gluons at 8 TeV.
gslpp::complex B0_MZ2_mHl2_mA2_mA2(const double MZ2, const double mHl2, const double mA2) const
gslpp::matrix< double > ATLAS13_pp_H_hh_gagaWW
double THoEX_ggF_H_ZZ_llll_ATLAS13
double ie25c[2][CacheSize]
double ip_ex_gg_A_hZ_bbll_CMS8_cache[2][CacheSize]
double ip_ex_bb_phi_tautau_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a bottom quark produced scalar resonance d...
double rh_gaga
Squared relative coupling of to two photons.
double ie34e(double mass)
gslpp::matrix< double > CMS8_gg_phi_gaga
double ip_ex_pp_phi_Zga_llga_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to a boson an...
double R_pp_H_Zga_llga_CMS8
double R_ggF_H_ZZ_llqq_ATLAS13
double ip_cs_ggtoA_13_cache[2][CacheSize]
double ip_ex_bb_phi_tt_ATLAS13_cache[2][CacheSize]
gslpp::complex B0_MZ2_0_mA2_mHl2(const double MZ2, const double mA2, const double mHl2) const
gslpp::complex A_H_W(const double mH, const double cW2, const double MW, const double MZ) const
gslpp::complex B0p_MZ2_0_mHp2_mA2(const double MZ2, const double mHp2, const double mA2) const
double ie34ec[2][CacheSize]
double ie1c[2][CacheSize]
double ip_ex_pp_phi_hh_bbbb_CMS8_cache_e[2][CacheSize]
gslpp::complex I_H_Hp(const double mHp2, const double mH) const
double THoEX_pp_A_gaga_ATLAS13
gslpp::matrix< double > CMS13_VBF_H_ZZ_llnunu
double Br_HtoAZ
branching ratio to two light Higgs bosons.
double ip_ex_pp_H_hh_bbbb_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
double pp_H_ZZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_ex_VBF_H_ZZ_llnunu_CMS13_cache_e[2][CacheSize]
double ip_ex_gg_H_ZZ_llnunu_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a gluon-gluon produced scalar resonance de...
gslpp::matrix< double > ATLAS13_pp_Hp_tb2
gslpp::matrix< double > CMS13_bb_phi_tautau
double ip_ex_bsgamma(double logtb, double logmHp)
Interpolating function for the NNLO value for the branching ratio of decays in the THDM...
double ip_ex_gg_phi_tautau_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a gluon-gluon produced scalar resonance de...
gslpp::matrix< double > temp4e
double THoEX_VBF_H_ZZ_ATLAS8
gslpp::complex B0_MZ2_mA2_mA2_mHh2_cache[4][CacheSize]
void computeHHquantities()
gslpp::matrix< double > ATLAS13_pp_H_VV_qqqq_e
gslpp::complex B0_MZ2_0_MW2_mHl2_cache[4][CacheSize]
gslpp::matrix< double > temp30
double ie19c[2][CacheSize]
double computeThValue()
The quartic coupling .
gslpp::complex I_A_L(const double mA2, const double Mmu, const double Mtau) const
double R_pp_Hp_tb_ATLAS13_1
double ie23c[2][CacheSize]
double ip_ex_gg_phi_tautau_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a gluon-gluon produced scalar resonance de...
gslpp::complex B0_MZ2_mHl2_mHh2_mHh2_cache[4][CacheSize]
gslpp::complex B0_MZ2_mHl2_0_0(const double MZ2, const double mHl2) const
double ip_ex_pp_Hpm_tb_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a singly charged scalar resonance decayi...
double ip_ex_gg_A_hZ_tautaull_CMS8_cache[2][CacheSize]
double ie10c[2][CacheSize]
double ip_ex_pp_phi_bb_CMS13_cache_e[2][CacheSize]
double ip_ex_gg_H_ZZ_llnunu_ATLAS13_cache[2][CacheSize]
double ip_ex_pp_phi_Zga_llga_ATLAS13_cache_e[2][CacheSize]
gslpp::complex B0p_MZ2_mHh2_0_mHp2_cache[4][CacheSize]
gslpp::complex A_H_Hp(const double mHp2, const double mH, const double cW2, const double MZ) const
double ip_ex_pp_phi_Zga_llga_CMS13_cache[2][CacheSize]
gslpp::matrix< double > ATLAS13_gg_H_WW_enumumu_e
double ip_ex_pp_phi_bb_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to a b quark p...
double bbF_A_tt_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::matrix< double > log_cs_bbA_8
double R_pp_H_Zga_llga_ATLAS8
double pp_H_hh_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::matrix< double > ATLAS13_VBF_H_ZZ_llllnunu
gslpp::complex B0p_MZ2_0_mA2_mHl2(const double MZ2, const double mA2, const double mHl2) const
gslpp::complex B0_MZ2_MZ2_MZ2_mHl2_cache[3][CacheSize]
double ip_cs_ZtoZH_8_cache[2][CacheSize]
gslpp::matrix< double > CMS13_pp_H_hh_bbtautau1_e
double ip_ex_VBF_H_WW_enumunu_ATLAS13_cache_e[2][CacheSize]
double ip_cs_ggtoH_8(double mass)
Interpolating function for the H production cross section via gluon-gluon fusion at 8 TeV...
double ip_ex_pp_H_hh_bbtautau_CMS13_cache_e[2][CacheSize]
double ie28e(double mass)
double ie16c[2][CacheSize]
gslpp::complex B00_MZ2_MZ2_mHl2_mA2(const double MZ2, const double mHl2, const double mA2) const
gslpp::matrix< double > CMS13_pp_H_hh_gagabb_e
gslpp::matrix< double > CMS13_ggF_H_hh_bbbb_e
double ip_ex_VBF_H_ZZ_llll_ATLAS13_cache_e[2][CacheSize]
double ip_ex_gg_phi_gaga_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a gluon-gluon produced scalar resonance de...
gslpp::matrix< double > ATLAS13_gg_H_ZZ_llllnunu
gslpp::matrix< double > log_cs_ttH_13
SM Higgs production cross section table at 13 TeV obtained with MadGraph 5, depending on the Higgs ma...
gslpp::complex I_h_L(const double mHl2, const double Me, const double Mmu, const double Mtau) const
double ie37e(double mass)
double pp_H_hh_bbbb_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double THoEX_pp_H_hh_bbbb_CMS13
double R_ggF_A_hZ_bbZ_ATLAS8
double ip_ex_pp_H_hh_gagabb_CMS13_cache[2][CacheSize]
gslpp::complex B0_MZ2_mHh2_mHh2_mHh2_cache[3][CacheSize]
double ie20ec[2][CacheSize]
gslpp::complex B0p_MZ2_0_mA2_mHh2_cache[4][CacheSize]
double pp_H_hh_gagabb_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double R_ggF_H_gaga_CMS13
double THoEX_bbF_H_tt_ATLAS13
double pp_Hpm_taunu_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::matrix< double > temp10e
double ip_ex_pp_H_hh_bbtautau1_CMS13_cache_e[2][CacheSize]
double ie35c[2][CacheSize]
double ie33e(double mass)
gslpp::matrix< double > ATLAS13_VBF_H_ZZ_qqllnunu_e
double Gamma_h
Total h decay rate in the THDM.
double KaellenFunction(const double a2, const double b2, const double c2) const
Kaellen function.
double computeSigmaggH_bb(const double sqrt_s) const
The square of the bottom-quark contribution to the ggH cross section in the Standard Model...
gslpp::matrix< double > temp37e
gslpp::complex B0p_MZ2_mA2_mA2_mHh2_cache[4][CacheSize]
double ggF_A_gaga_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::complex B0p_MZ2_mHh2_0_mA2(const double MZ2, const double mHh2, const double mA2) const
gslpp::matrix< double > br_cc
gslpp::complex B0_MZ2_mA2_mA2_mHh2(const double MZ2, const double mA2, const double mHh2) const
double cW2THDM(const double c02) const
gslpp::complex B0p_MZ2_0_mHp2_mHh2(const double MZ2, const double mHp2, const double mHh2) const
double ip_ex_gg_H_hh_bbtautau_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a gluon-gluon produced scalar resonance de...
double ip_ex_gg_H_hh_bbtautau_CMS8_cache[2][CacheSize]
gslpp::complex B0p_MZ2_mHp2_0_mHh2(const double MZ2, const double mHp2, const double mHh2) const
double ip_ex_pp_H_hh_bbVV_CMS13_cache[2][CacheSize]
const double & getMass() const
A get method to access the particle mass.
double ip_ex_bb_phi_tautau_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a b associated scalar resonance decaying t...
gslpp::matrix< double > CMS13_gg_phi_gaga_e
double ip_ex_gg_A_hZ_bbZ_ATLAS8_cache[2][CacheSize]
gslpp::matrix< double > ATLAS13_gg_H_ZZ_llqq_e
gslpp::matrix< double > CMS13_ggFVBF_H_WW_lnulnu_e
gslpp::complex B0_MZ2_0_0_mHh2_cache[3][CacheSize]
gslpp::matrix< double > ATLAS13_gg_H_ZZ_llqq
double ip_ex_pp_Hpm_taunu_ATLAS8_cache[2][CacheSize]
gslpp::matrix< double > CMS13_VBF_H_ZZ_llnunu_e
double pp_A_gaga_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double interpolate2D(gslpp::matrix< double > arrayTab, double x, double y)
Linearly interpolates a table with two parameter dimensions.
gslpp::matrix< double > CMS13_pp_H_hh_bbVV_e
gslpp::matrix< double > CMS8_pp_H_AZ_tautaull
gslpp::matrix< double > temp3e
double THoEX_pp_H_hh_bbgaga_CMS13
double VBF_Vh
Ratio of THDM and SM cross sections for VBF and Vh production of h.
gslpp::complex B00(const double mu2, const double p2, const double m02, const double m12) const
.
double ie3c[2][CacheSize]
double ip_Br_HPtomumu(double mass)
Interpolating function for the SM branching ratio to two muons.
An RGE running algorithm.
gslpp::complex B0_MZ2_0_mA2_mHh2(const double MZ2, const double mA2, const double mHh2) const
gslpp::matrix< double > log_cs_bbA_13
gslpp::matrix< double > CMS8_bb_phi_bb
double ip_ex_gg_A_hZ_tautaull_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a gluon-gluon produced pseudoscalar resona...
gslpp::complex B0p_MZ2_mHh2_0_0(const double MZ2, const double mHh2) const
double ip_ex_ggF_H_hh_bbbb_CMS13_cache_e[2][CacheSize]
gslpp::matrix< double > temp31e
double R_pp_Hp_taunu_CMS8
double THoEX_pp_H_ZZ_llnunu_CMS13
double ip_ex_gg_H_ZZ_llll_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
gslpp::complex B00_MZ2_MZ2_MZ2_mHl2_cache[3][CacheSize]
double THoEX_pp_A_Zga_qqga_CMS13
double ip_ex_bb_phi_tautau_ATLAS13_cache_e[2][CacheSize]
gslpp::matrix< double > CMS13_pp_H_hh_bbtautau
double THoEX_VBF_H_ZZ_qqllnunu_ATLAS13
double ip_ex_pp_H_hh_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
double ip_Br_HPtobb(double mass)
Interpolating function for the SM branching ratio to two bottom quarks.
double ie29ec[2][CacheSize]
gslpp::matrix< double > CMS13_pp_H_hh_bblnulnu_e
double ip_ex_gg_phi_tautau_ATLAS8_cache[2][CacheSize]
double R_pp_Hpm_tb_ATLAS8
double ie27c[2][CacheSize]
double ip_ex_pp_phi_gaga_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a scalar resonance decaying to two photo...
gslpp::complex I_HH_U_cache[4][CacheSize]
double ip_ex_ggVV_H_WW_lnulnu_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a gluon-gluon or vector boson fusion produ...
gslpp::complex B0_MZ2_MZ2_MZ2_mHl2(const double MZ2, const double mHl2) const
gslpp::complex B0_MZ2_mHl2_0_mA2_cache[4][CacheSize]
gslpp::matrix< double > temp5
double ip_ex_gg_H_ZZ_nunuqq_ATLAS13_cache[2][CacheSize]
gslpp::matrix< double > log_cs_ggA_13
CP-odd Higgs production cross section tables at 13 TeV obtained with HIGLU 4.34, depending on the Hig...
double ip_ex_ggF_H_hh_bbbb_CMS13_cache[2][CacheSize]
double ie20e(double mass)
gslpp::matrix< double > CMS13_pp_phi_bb_e
CMS expected upper cross section limits at 13 TeV, depending on the Higgs mass.
double ip_GammaHPtotSM_cache[2][CacheSize]
gslpp::matrix< double > CMS13_ggFVBF_H_WW_lnulnu
double R_pp_H_ZZ_llnunu_CMS13
double R_pp_H_ZZ_llll_CMS13
gslpp::matrix< double > ATLAS13_VBF_H_WW_lnuqq_e
double THoEX_pp_A_Zga_llga_CMS13
double ie35ec[2][CacheSize]
double ip_ex_pp_H_ZZ_llnunu_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
double THoEX_pp_A_HZ_tautaull_CMS8
gslpp::complex B00_MZ2_0_MZ2_mHl2_cache[3][CacheSize]
double R_ggF_H_tautau_CMS13
gslpp::complex B0_MZ2_mHh2_mHh2_mHh2(const double MZ2, const double mHh2) const
gslpp::complex B0p_MZ2_mHl2_mHl2_mHl2(const double MZ2, const double mHl2) const
gslpp::matrix< double > temp14e
gslpp::complex B00_MZ2_0_mHp2_mHp2_cache[3][CacheSize]
gslpp::complex B00_MZ2_0_mA2_mHp2(const double MZ2, const double mA2, const double mHp2) const
double ip_ex_gg_H_WW_lnuqq_ATLAS13_cache_e[2][CacheSize]
double pp_H_gaga_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_ex_pp_Hpm_taunu_CMS13_cache_e[2][CacheSize]
gslpp::matrix< double > ATLAS13_pp_H_VV_qqqq
double THoEX_VBF_H_ZZ_llllnunu_ATLAS13
double R_pp_H_gaga_ATLAS13
double ip_ex_VBF_H_ZZ_llqq_ATLAS13_cache_e[2][CacheSize]
double ip_ex_pp_H_ZZ_llnunu_CMS13_cache_e[2][CacheSize]
gslpp::matrix< double > temp4
gslpp::matrix< double > temp32e
std::string getRGEorderflag() const
THDM branching ratio of .
gslpp::complex B0p_MZ2_mHh2_mA2_mA2(const double MZ2, const double mHh2, const double mA2) const
double ie33c[2][CacheSize]
gslpp::matrix< double > ATLAS13_gg_H_WW_lnuqq
double getBDtaunu_A() const
gslpp::complex B0_MZ2_mHl2_mHp2_mHp2_cache[4][CacheSize]
double R_bbF_H_tautau_ATLAS8
double THoEX_pp_H_hh_CMS8
double R_ggF_A_tautau_ATLAS8
gslpp::matrix< double > ATLAS8_gg_phi_tt
double ip_ex_pp_H_ZZ_llll_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
gslpp::matrix< double > CMS13_pp_phi_Zga_qqga
gslpp::matrix< double > ATLAS13_VBF_H_WW_enumumu_e
double THoEX_ggF_H_ZZ_llnunu_ATLAS13
gslpp::matrix< double > log_cs_ggH_13
SM Higgs production cross section tables at 13 TeV from the LHC Higgs Cross Section Working Group...
double bbF_H_bb_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ip_cs_ggtoH_13(double mass)
Interpolating function for the H production cross section via gluon-gluon fusion at 13 TeV...
gslpp::matrix< double > ATLAS13_gg_H_ZZ_llnunu
double THoEX_ggF_H_ZZ_llqq_ATLAS13
gslpp::matrix< double > csrA_bottom_13
double ip_ex_pp_Hpm_tb_ATLAS8_cache[2][CacheSize]
gslpp::complex B0p_MZ2_mA2_0_mHp2_cache[4][CacheSize]
gslpp::complex B0p_MZ2_mHh2_mHl2_mHl2_cache[4][CacheSize]
double KaellenFunction_cache[4][CacheSize]
gslpp::matrix< double > CMS13_pp_H_hh_bblnulnu
double ip_ex_pp_phi_gaga_ATLAS13_cache[2][CacheSize]
gslpp::complex B00_MZ2_0_mHh2_mA2_cache[4][CacheSize]
gslpp::matrix< double > log_cs_ggHp_8
Charged Higgs production cross section table at 8 TeV from LHCHXSWGMSSMCharged, depending on the char...
gslpp::matrix< double > temp24
gslpp::matrix< double > br_bb
void read()
Fills all required arrays with the values read from the tables.
double ip_ex_gg_H_ZZ_llllnunu_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double pp_H_hh_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::complex B0_MZ2_mHh2_0_0_cache[3][CacheSize]
gslpp::complex B0p_MZ2_mHl2_mHh2_mHh2(const double MZ2, const double mHl2, const double mHh2) const
gslpp::complex B00_MZ2_MZ2_mHp2_mHp2(const double MZ2, const double mHp2) const
double ip_ex_gg_phi_Zga_llga_ATLAS13_cache[2][CacheSize]
double ip_ex_pp_H_AZ_bbll_CMS8(double mA, double mH)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to which furt...
double R_ggF_H_ZZ_llllnunu_ATLAS13
gslpp::complex B0_MZ2_mHl2_0_mA2(const double MZ2, const double mHl2, const double mA2) const
double ip_ex_ggVV_H_WW_lnulnu_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a gluon-gluon or vector boson fusion produ...
double ip_ex_ggF_phi_Zga_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a gluon-gluon produced scalar resonance de...
double ip_ex_gg_H_ZZ_llll_ATLAS13_cache_e[2][CacheSize]
double R_pp_H_Zga_llga_CMS13
double rh_QdQd
Squared relative coupling of to two down type quarks.
gslpp::complex B0p_MZ2_mA2_0_mHh2(const double MZ2, const double mA2, const double mHh2) const
double THoEX_ggF_H_hh_ATLAS8
double R_VBF_H_ZZ_llqq_ATLAS13
double ie19e(double mass)
double ie26e(double mass)
gslpp::matrix< double > temp21e
double ie2c[2][CacheSize]
double ip_ex_gg_H_ZZ_qqllnunu_ATLAS13_cache_e[2][CacheSize]
double ip_cs_pptobbH_13(double mass)
Interpolating function for the bottom associated H production cross section at 13 TeV...
gslpp::complex B0p_MZ2_mHl2_mHl2_mHl2_cache[3][CacheSize]
gslpp::matrix< double > temp5e
gslpp::complex B0_MZ2_mHh2_0_mA2(const double MZ2, const double mHh2, const double mA2) const
double computeSigmaggH_tt(const double sqrt_s) const
The square of the top-quark contribution to the ggH cross section in the Standard Model...
double R_ggF_A_tautau_CMS8
double ie22c[2][CacheSize]
gslpp::complex B0p_MZ2_mA2_mA2_mHl2(const double MZ2, const double mA2, const double mHl2) const
double ip_ex_bb_phi_tautau_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a bb associated scalar resonance decayin...
double THoEX_bbF_H_tautau_CMS13
gslpp::complex I_HH_L_cache[4][CacheSize]
gslpp::matrix< double > temp28e
double ip_ex_pp_H_VV_qqqq_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a scalar resonance decaying to two or ...
gslpp::complex B00_MZ2_0_mHl2_mHp2_cache[4][CacheSize]
gslpp::complex f_func(const double x) const
f function for the gamma gamma coupling to h, H and A
double ip_csr_ggH_b_8(double mass)
Interpolating function for the gluon-gluon fusion H cross section ratio of the bottom-loop and the to...
double VBF_H_ZZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_ex_gg_H_WW_enumunu_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double ie18e(double mass)
double ie4c[2][CacheSize]
gslpp::matrix< double > temp27
double ip_ex_gg_H_ZZ_llnunu_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double ie12ec[2][CacheSize]
double ip_ex_gg_A_hZ_tautauZ_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced pseudoscalar reso...
gslpp::matrix< double > log_cs_ttH_8
SM Higgs production cross section table at 8 TeV obtained with MadGraph 5, depending on the Higgs mas...
double ip_ex_pp_phi_hh_gagabb_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
double VBF_H_ZZ_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::matrix< double > CMS13_VBFVH_H_ZZ_llll_e
double ip_ex_VBF_VH_H_ZZ_llll_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a vector boson fusion produced scalar reso...
gslpp::matrix< double > temp29
double R_ggF_H_tautau_CMS8
gslpp::matrix< double > CMS8_pp_Hp_tb
gslpp::matrix< double > ATLAS8_gg_A_hZ_tautauZ_e
double rh_Zga
Squared relative coupling of to a boson and a photon.
double ip_ex_pp_phi_gaga_ATLAS13_cache_e[2][CacheSize]
gslpp::complex B0p_MZ2_0_mHp2_mHl2_cache[4][CacheSize]
gslpp::complex B0_MZ2_mHp2_mHp2_mHl2(const double MZ2, const double mHp2, const double mHl2) const
double ip_csr_ggH_b_13(double mass)
Interpolating function for the gluon-gluon fusion H cross section ratio of the bottom-loop and the to...
double ip_cs_pptottH_13(double mass)
Interpolating function for the top associated H production cross section at 13 TeV.
gslpp::matrix< double > ATLAS8_pp_phi_gaga_e
ATLAS expected upper cross section limits at 8 TeV, depending on the Higgs mass. ...
double THoEX_bbF_H_tautau_CMS8
void computeSignalStrengthQuantities()
gslpp::complex B00_MZ2_0_MZ2_mHh2_cache[3][CacheSize]
double R_bbF_H_tautau_ATLAS13
double R_bbF_H_tautau_CMS13
double ie11e(double mass)
double ie39c[2][CacheSize]
double ie40ec[2][CacheSize]
double THoEX_VBF_H_WW_enumunu_ATLAS13
double pp_H_hh_bblnulnu_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_ex_pp_Hp_tb_CMS8_cache[2][CacheSize]
gslpp::matrix< double > ATLAS8_gg_A_hZ_bbZ
double ip_ex_pp_H_hh_gagaWW_ATLAS13_cache[2][CacheSize]
double ie5c[2][CacheSize]
gslpp::complex B0_MZ2_MW2_MW2_mHh2_cache[4][CacheSize]
gslpp::matrix< double > temp8
gslpp::complex A_A_L(const double mA2, const double cW2, const double Mmu, const double Mtau, const double MZ) const
gslpp::matrix< double > csrA_top_8
gslpp::complex B0_MZ2_mHp2_mHp2_mHh2(const double MZ2, const double mHp2, const double mHh2) const
gslpp::complex B0_MZ2_0_mHp2_mHh2_cache[4][CacheSize]
double ip_ex_gg_phi_tt_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double VBF_H_ZZ_llll_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_ex_gg_H_WW_lnuqq_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double pph13
Ratio of THDM and SM cross sections for the production of h at 13 TeV.
double ip_ex_bb_phi_tautau_CMS8_cache_e[2][CacheSize]
double bbF_A_tautau_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::complex B0_MZ2_0_MW2_mHh2_cache[4][CacheSize]
double THoEX_pp_A_HZ_bbll_CMS8
double THoEX_ggF_A_gaga_CMS13
gslpp::complex B0_MZ2_0_MZ2_mHh2_cache[3][CacheSize]
double ip_ex_gg_H_hh_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double ip_ex_gg_phi_tautau_CMS8_cache[2][CacheSize]
double R_ggF_H_ZZ_llnunu_ATLAS13
gslpp::matrix< double > csrH_top_13
Production cross section ratio tables at 13 TeV obtained with HIGLU 4.34, depending on the Higgs mass...
double ip_ex_VBF_H_WW_lnuqq_ATLAS13_cache[2][CacheSize]
gslpp::matrix< double > temp19e
An observable class for the quartic Higgs potential coupling .
double ip_ex_pp_Hp_tb_ATLAS13_1_cache_e[2][CacheSize]
gslpp::matrix< double > br_mumu
double ip_ex_VBF_H_ZZ_llnunu_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a VBF produced scalar resonance decaying t...
double ip_ex_gg_phi_tautau_ATLAS8_cache_e[2][CacheSize]
double THoEX_bbF_H_tautau_ATLAS8
double ip_ex_pp_H_hh_bbbb_CMS13_cache[2][CacheSize]
double ip_ex_pp_H_ZZ_llll_CMS13_cache_e[2][CacheSize]
double R_pp_H_hh_bbgaga_ATLAS13
double ip_ex_pp_A_Zga_llga_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a pseudoscalar resonance decaying to a Z b...
gslpp::matrix< double > CMS8_pp_H_hh_bbbb_e
double g_HH_HpHm(const double mHp2, const double mHh2, const double tanb, const double m12_2, const double bma, const double vev) const
gslpp::complex B0p(const double muIR2, const double p2, const double m02, const double m12) const
.
double R_ggF_A_tautau_CMS13
double R_VBF_H_ZZ_llll_ATLAS13
A base class for symmetric Two-Higgs-Doublet models.
double GammaAtot
Total decay width of the CP-odd Higgs .
gslpp::matrix< double > CMS8_pp_A_Zga_llga
double ip_ex_gg_phi_tt_ATLAS8_cache[2][CacheSize]
gslpp::complex A_h_D_cache[7][CacheSize]
double ip_ex_gg_phi_tautau_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double ip_ex_gg_A_hZ_bbll_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a gluon-gluon produced pseudoscalar resona...
double ip_ex_pp_Hpm_tb_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a singly charged scalar resonance decayi...
double THoEX_ggF_H_WW_ATLAS8
gslpp::complex B00_MZ2_MW2_mHp2_mHp2(const double MZ2, const double MW2, const double mHp2) const
gslpp::matrix< double > temp24e
double THoEX_VBF_H_ZZ_llnunu_CMS13
gslpp::complex B00_MZ2_MW2_mHh2_mHp2_cache[5][CacheSize]
double ip_ex_pp_Hpm_taunu_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a singly charged scalar resonance decayi...
gslpp::complex I_A_U(const double mA2, const double Mc, const double Mt) const
double pp_Hp_tb_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_ex_pp_phi_Zga_llga_CMS13_cache_e[2][CacheSize]
double pp_A_Zga_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double R_pp_Hp_tb_ATLAS13_2
gslpp::matrix< double > CMS13_pp_Hpm_taunu_e
CMS expected upper cross section limits at 13 TeV, depending on the charged Higgs mass...
double ie14ec[2][CacheSize]
gslpp::matrix< double > ATLAS8_pp_Hpm_taunu_e
ATLAS expected upper cross section limits at 8 TeV, depending on the charged Higgs mass...
double THoEX_ggF_H_tt_ATLAS8
gslpp::matrix< double > temp26e
gslpp::matrix< double > ATLAS13_gg_H_ZZ_llllnunu_e
virtual double v() const
The value of at any scale with the number of flavours and full EW corrections. ...
double bbF_A_tautau_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ie39e(double mass)
gslpp::complex Int2(const double tau, const double lambda) const
function for Z gamma coupling to h, H and A
gslpp::matrix< double > CMS8_gg_A_hZ_tautaull
double ip_ex_gg_H_ZZ_llll_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double ip_ex_pp_Hp_tb_ATLAS13_1(double mass)
Interpolating function for the observed ATLAS upper limit on a singly charged scalar resonance decayi...
double ip_ex_gg_H_ZZ_nunuqq_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double THoEX_pp_Hpm_tb_ATLAS8
THDMcache(const StandardModel &SM_i)
THDMcache constructor.
gslpp::matrix< double > ATLAS13_pp_phi_Zga
gslpp::matrix< double > CMS13_gg_phi_tautau
double ip_ex_pp_phi_hh_gagabb_CMS8_cache[2][CacheSize]
double R_pp_H_hh_bbbb_CMS8
gslpp::matrix< double > CMS13_pp_H_ZZ_llll_e
gslpp::complex B0_MZ2_mHl2_mHl2_mHl2(const double MZ2, const double mHl2) const
double ip_ex_gg_H_WW_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double Br_AtoHpW
branching ratio to a charged Higgs and a boson.
double ie29e(double mass)
double ip_cs_ggtoHp_13_cache[3][CacheSize]
double R_VBF_H_ZZ_llllnunu_ATLAS13
double ip_ex_pp_H_hh_gagabb_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
double THoEX_pp_Hp_tb_CMS8
double ie10e(double mass)
gslpp::matrix< double > ATLAS13_bb_A_Zh_Zbb_e
gslpp::complex A_A_D(const double mA2, const double cW2, const double Ms, const double Mb, const double MZ) const
double pp_H_hh_bbVV_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_ex_pp_H_hh_gagabb_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a scalar resonance decaying to two boso...
double ip_ex_VBF_VH_H_ZZ_llll_CMS13_cache_e[2][CacheSize]
gslpp::complex B00_MZ2_0_mA2_mHp2_cache[4][CacheSize]
gslpp::matrix< double > temp6
double ip_ex_gg_A_Zh_Zbb_ATLAS13_cache_e[2][CacheSize]
double R_pp_A_Zga_llga_ATLAS13
double ip_ex_VBF_VH_H_ZZ_llll_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a vector boson fusion produced scalar reso...
gslpp::complex B0p_MZ2_mHp2_mHp2_mHh2_cache[4][CacheSize]
double R_ggF_A_tautau_ATLAS13
gslpp::complex A_H_Hp_cache[5][CacheSize]
double R_pp_Hp_tb_ATLAS13
double ip_cs_VBFtoH_13_cache[2][CacheSize]
gslpp::complex B0p_MZ2_0_mA2_mHh2(const double MZ2, const double mA2, const double mHh2) const
double Br_AtohZ
branching ratio to an and a boson.
gslpp::complex B0_MZ2_mHp2_mHp2_mHh2_cache[4][CacheSize]
double THoEX_pp_H_gaga_ATLAS8
double ip_ex_gg_phi_tautau_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
gslpp::matrix< double > ATLAS13_pp_Hp_tb1
gslpp::matrix< double > temp26
gslpp::complex B0_MZ2_MW2_MW2_mHl2_cache[4][CacheSize]
gslpp::complex B0p_MZ2_mHl2_mHh2_mHl2(const double MZ2, const double mHl2, const double mHh2) const
double THoEX_ggF_H_ZZ_qqllnunu_ATLAS13
double ggF_A_hZ_tautauZ_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double R_ggF_A_hZ_bbll_CMS8
double ip_cs_ggtoA_13(double mass)
Interpolating function for the A production cross section via gluon-gluon fusion at 13 TeV...
gslpp::complex B00_MZ2_MW2_mHl2_mHp2_cache[5][CacheSize]
double ip_ex_gg_H_ZZ_nunuqq_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double ie14e(double mass)
gslpp::complex B00_MZ2_0_mHl2_mHp2(const double MZ2, const double mHl2, const double mHp2) const
double ip_ex_pp_phi_bb_CMS13_cache[2][CacheSize]
double THoEX_pp_A_Zga_llga_ATLAS13
double computeThValue()
The quartic coupling .
double ie2ec[2][CacheSize]
double ip_csr_ggH_t_13(double mass)
Interpolating function for the gluon-gluon fusion H cross section ratio of the top-loop and the total...
double ip_csr_ggH_b_13_cache[2][CacheSize]
double ip_ex_bb_phi_bb_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a bottom quark produced scalar resonance d...
double THoEX_ggF_A_tautau_CMS8
double THoEX_pp_H_AZ_bbll_CMS8
double ggF_H_hh_bbbb_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double ip_ex_pp_phi_hh_bbbb_CMS8_cache[2][CacheSize]
gslpp::matrix< double > ATLAS13_pp_H_hh_gagabb
gslpp::matrix< double > CMS8_pp_A_HZ_bbll
gslpp::matrix< double > ATLAS8_pp_Hpm_tb
gslpp::complex B0_MZ2_mHh2_0_mHp2_cache[4][CacheSize]
gslpp::matrix< double > ATLAS8_pp_phi_Zga_llga_e
double R_bbF_A_tautau_CMS8
gslpp::matrix< double > ATLAS8_gg_H_WW_e
double ip_ex_gg_phi_Zga_llga_ATLAS13_cache_e[2][CacheSize]
double ie24e(double mass)
double ip_cs_ggtoH_8_cache[2][CacheSize]
gslpp::matrix< double > ATLAS13_pp_phi_Zga_e
gslpp::matrix< double > temp18e
double ip_Br_HPtocc_cache[2][CacheSize]
gslpp::complex B0p_MZ2_0_0_mHh2_cache[3][CacheSize]
double THoEX_pp_H_Zga_llga_ATLAS8
gslpp::matrix< double > temp11
double R_pp_A_Zga_qqga_CMS13
double ip_ex_pp_Hp_tb_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a singly charged scalar resonance decaying...
double getBDtaunu_B() const
gslpp::complex B0p_MZ2_0_0_mHl2_cache[3][CacheSize]
double ie33ec[2][CacheSize]
gslpp::complex B0p_MZ2_mHl2_0_mHp2(const double MZ2, const double mHl2, const double mHp2) const
gslpp::matrix< double > ATLAS13_VBF_H_ZZ_llqq
double ip_ex_pp_Hp_taunu_CMS8_e(double mass)
Interpolating function for the expected CMS upper limit on a singly charged scalar resonance decaying...
double ip_ex_gg_phi_tt_ATLAS8_cache_e[2][CacheSize]
double ip_ex_bb_A_Zh_Zbb_ATLAS13_cache[2][CacheSize]
gslpp::matrix< double > ATLAS13_pp_phi_gaga_e
double pp_A_HZ_bbll_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double THoEX_pp_H_hh_bbbb_CMS8
virtual double RGERunner(double InitialValues[], unsigned long int NumberOfRGEs, double Q1, double Q2, int order, double Rpeps, double tNLOuni)
gslpp::complex A_h_L(const double mHl2, const double cW2, const double Me, const double Mmu, const double Mtau, const double MZ) const
double getBDstartaunu_B() const
gslpp::complex B0_MZ2_mHp2_0_mHl2(const double MZ2, const double mHp2, const double mHl2) const
double ip_ex_tt_phi_tt_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a tt associated scalar resonance decayin...
double ip_ex_gg_phi_tautau_CMS13_cache[2][CacheSize]
gslpp::complex B0p_MZ2_mHl2_mA2_mA2(const double MZ2, const double mHl2, const double mA2) const
double THoEX_pp_H_hh_bbgaga_ATLAS13
double ip_ex_pp_H_ZZ_llqq_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
double ggF_H_ZZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
complex log(const complex &z)
double THoEX_ggF_VBF_H_WW_lnulnu_CMS13
double ip_ex_pp_H_hh_bblnulnu_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
virtual double Mw_tree() const
The tree-level mass of the boson, .
gslpp::complex B0_MZ2_mHh2_mHh2_mHl2_cache[4][CacheSize]
double ip_ex_pp_H_hh_gagabb_ATLAS13_cache_e[2][CacheSize]
int CacheCheck(const gslpp::complex cache[][CacheSize], const int NumPar, const double params[]) const
Check whether for the latest set of parameters a value is in the cache.
double ip_ex_gg_H_ZZ_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double R_pp_H_hh_gagabb_CMS8
double THoEX_ggF_H_hh_bbtautau_CMS8
gslpp::complex A_HH_L(const double mHh2, const double cW2, const double Mmu, const double Mtau, const double MZ) const
double ip_cs_pptobbH_8(double mass)
Interpolating function for the bottom associated H production cross section at 8 TeV.
double R_pp_Hpm_taunu_ATLAS13
double ggF_H_gaga_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::matrix< double > log_cs_bbH_13
SM Higgs production cross section table at 13 TeV obtained with SusHi 1.5, depending on the Higgs mas...
double THDM_BR_h_bb
branching ratio to two quarks in the THDM.
gslpp::matrix< double > temp7
gslpp::matrix< double > csrH_bottom_8
gslpp::matrix< double > temp40e
gslpp::complex B0p_MZ2_mHl2_mHh2_mHl2_cache[4][CacheSize]
double ip_ex_pp_phi_Zga_llga_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a scalar resonance decaying to a boson ...
double ip_ex_gg_phi_tautau_CMS13_cache_e[2][CacheSize]
gslpp::matrix< double > temp39
double THoEX_ggF_H_tautau_CMS13
double pp_H_bb_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double R_VBF_H_WW_lnuqq_ATLAS13
double THoEX_pp_H_hh_gagabb_CMS8
gslpp::matrix< double > CMS13_pp_H_hh_bbtautau_e
gslpp::matrix< double > ATLAS8_gg_A_hZ_tautauZ
double ip_ex_VBF_H_WW_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a vector boson fusion produced scalar re...
double ie22e(double mass)
gslpp::matrix< double > ATLAS8_VBF_H_ZZ
gslpp::matrix< double > ATLAS13_gg_H_ZZ_nunuqq
double ip_ex_bb_phi_tt_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a bb associated scalar resonance decayin...
double ip_cs_ggtoHp_8_cache[3][CacheSize]
double ip_ex_gg_A_Zh_Zbb_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced pseudoscalar reso...
double R_pp_A_HZ_bbll_CMS8
double ip_ex_gg_phi_gaga_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a gluon-gluon produced scalar resonance de...
double R_bbF_A_tt_ATLAS13
double ie27ec[2][CacheSize]
gslpp::complex B0p_MZ2_mHp2_0_mHh2_cache[4][CacheSize]
double ie3ec[2][CacheSize]
gslpp::complex I_h_D(const double mHl2, const double Md, const double Ms, const double Mb) const
gslpp::matrix< double > ATLAS13_pp_Hp_tb2_e
double ip_ex_pp_phi_Zga_llga_ATLAS8_cache_e[2][CacheSize]
double ttF_A_tt_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double R_ggF_VBF_H_WW_lnulnu_CMS13
double ip_Br_HPtomumu_cache[2][CacheSize]
An observable class for the quadratic Higgs potential coupling .
double ip_ex_pp_A_HZ_bbll_CMS8(double mA, double mH)
Interpolating function for the observed CMS upper limit on a pseudoscalar resonance decaying to whic...
gslpp::matrix< double > temp14
gslpp::complex B00_MZ2_MW2_MW2_mHh2(const double MZ2, const double MW2, const double mHh2) const
double ip_ex_ggF_H_hh_bbbb_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
gslpp::matrix< double > CMS8_mu_pp_H_VV_e
CMS expected upper signal strength limits at 8 TeV, depending on the Higgs mass. ...
gslpp::complex B0_MZ2_0_MZ2_mHh2(const double MZ2, const double mHh2) const
double ip_ex_gg_A_hZ_bbZ_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced pseudoscalar reso...
gslpp::complex B0_MZ2_0_mHp2_mHl2_cache[4][CacheSize]
double ip_ex_VBF_H_ZZ_llll_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a vector boson fusion produced scalar re...
double ip_ex_gg_H_ZZ_llnunu_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double ip_ex_gg_phi_tt_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double R_bbF_A_tautau_CMS13
double ip_ex_pp_A_Zga_llga_CMS8_cache[2][CacheSize]
double R_pp_H_ZZ_llqq_CMS13
double ip_Br_HPtoWW_cache[2][CacheSize]
gslpp::complex B00_MZ2_MW2_MW2_mHl2(const double MZ2, const double MW2, const double mHl2) const
An observable class for the quartic Higgs potential coupling .
gslpp::complex B0_MZ2_0_MZ2_mHl2(const double MZ2, const double mHl2) const
double ip_ex_bb_phi_tautau_ATLAS13_cache[2][CacheSize]
gslpp::complex B0p_MZ2_mHl2_0_mA2_cache[4][CacheSize]
double pp_Hpm_tb_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::complex B0p_MZ2_mHl2_0_mA2(const double MZ2, const double mHl2, const double mA2) const
double ip_ex_pp_H_hh_bbbb_CMS13_cache_e[2][CacheSize]
gslpp::complex A_HH_U(const double mHh2, const double cW2, const double Mc, const double Mt, const double MZ) const
double pp_H_ZZ_llll_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::complex B0_MZ2_mHl2_0_0_cache[3][CacheSize]
double ie30e(double mass)
gslpp::matrix< double > temp12
gslpp::matrix< double > CMS13_gg_phi_tautau_e
double ip_ex_pp_Hp_tb_ATLAS13_1_e(double mass)
Interpolating function for the expected ATLAS upper limit on a singly charged scalar resonance decayi...
double ie24ec[2][CacheSize]
gslpp::matrix< double > CMS8_gg_A_hZ_bbll
CMS observed upper cross section limits at 8 TeV, depending on the Higgs mass.
gslpp::complex B0_MZ2_mHl2_mHl2_mHl2_cache[3][CacheSize]
double ie18c[2][CacheSize]
double THoEX_ggF_A_tt_ATLAS8
double ip_ex_bb_phi_tautau_ATLAS8_cache_e[2][CacheSize]
double rh_QuQu
Squared relative coupling of to two up type quarks.
double ie6c[2][CacheSize]
double ip_ex_bb_phi_tautau_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a b associated scalar resonance decaying t...
double ip_ex_pp_A_HZ_bbll_CMS8_cache[3][CacheSize]
double ggF_A_hZ_bbll_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::complex B00_MZ2_0_mHh2_mHp2_cache[4][CacheSize]
gslpp::matrix< double > CMS8_pp_H_hh_e
gslpp::complex B0p_MZ2_mA2_0_mHp2(const double MZ2, const double mA2, const double mHp2) const
double THoEX_ggF_H_hh_gagaWW_ATLAS13
gslpp::matrix< double > CMS13_pp_Hpm_taunu
CMS observed upper cross section limits at 13 TeV, depending on the charged Higgs mass...
double ip_Br_HPtoZZ(double mass)
Interpolating function for the SM branching ratio to two bosons.
gslpp::complex A_A_U_cache[6][CacheSize]
double ip_ex_bb_phi_tt_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a bb associated scalar resonance decayin...
gslpp::complex B00_MZ2_MW2_mHl2_mHp2(const double MZ2, const double MW2, const double mHl2, const double mHp2) const
double ie21ec[2][CacheSize]
gslpp::matrix< double > temp15
gslpp::matrix< double > ATLAS13_pp_H_hh_bbbb
double ie34c[2][CacheSize]
A class for defining operations on and functions of complex numbers.
gslpp::complex B0_MZ2_0_0_mHl2(const double MZ2, const double mHl2) const
double ggF_A_tt_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double THoEX_ggF_H_ZZ_llnunu_CMS13
gslpp::complex B0p_MZ2_mHh2_mHh2_mHl2(const double MZ2, const double mHh2, const double mHl2) const
gslpp::matrix< double > CMS13_pp_H_ZZ_llll
double ip_ex_pp_phi_hh_bbbb_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a scalar resonance decaying to two bosons...
gslpp::complex B0_MZ2_mHp2_0_mHl2_cache[4][CacheSize]
double ip_ex_bb_A_Zh_Zbb_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a bb associated pseudoscalar resonance d...
gslpp::complex B0p_MZ2_mHh2_mHp2_mHp2_cache[4][CacheSize]
double THoEX_VBF_H_WW_ATLAS8
double THDM_BR_h_WW
branching ratio to two bosons in the THDM.
double ip_csr_ggH_b_8_cache[2][CacheSize]
complex sin(const complex &z)
double pp_H_hh_gagabb_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double Br_HtoHpHm
branching ratio to two charged Higgs bosons.
double ie38e(double mass)
double ip_ex_gg_phi_gaga_CMS8_cache[2][CacheSize]
double ggF_A_Zga_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::matrix< double > GammaHtot_SM
Total SM decay width (obtained with HDECAY 6.10), depending on the Higgs mass.
double ie7ec[2][CacheSize]
double ip_csr_ggH_t_13_cache[2][CacheSize]
double ip_ex_gg_H_WW_enumunu_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double pp_H_hh_bbbb_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
gslpp::complex B0_MZ2_MW2_MW2_mHh2(const double MZ2, const double MW2, const double mHh2) const
gslpp::matrix< double > temp39e
gslpp::complex B0p_MZ2_mHh2_mHp2_mHp2(const double MZ2, const double mHh2, const double mHp2) const
gslpp::complex B00_MZ2_MZ2_mHp2_mHp2_cache[3][CacheSize]
double ip_cs_pptottA_8(double mass)
Interpolating function for the top associated A production cross section at 8 TeV.
double computeThValue()
The quartic coupling .
gslpp::complex B0p_MZ2_mHl2_mA2_mA2_cache[4][CacheSize]
double computeSigmattH(const double sqrt_s) const
The ttH production cross section in the Standard Model.
double Br_AtoHZ
branching ratio to an and a boson.
gslpp::matrix< double > temp32
double ip_ex_gg_H_WW_enumunu_ATLAS13_cache[2][CacheSize]
gslpp::matrix< double > ATLAS8_gg_H_hh
gslpp::complex B00_MZ2_MW2_MW2_mHl2_cache[4][CacheSize]
gslpp::matrix< double > CMS8_pp_H_hh_gagabb_e
double ip_ex_pp_Hp_taunu_CMS8(double mass)
Interpolating function for the observed CMS upper limit on a singly charged scalar resonance decaying...
double ggF_H_ZZ_llll_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double R_ggF_H_WW_lnuqq_ATLAS13
double ip_ex_VBF_H_WW_lnuqq_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a VBF produced scalar resonance decaying...
gslpp::matrix< double > temp9e
double c02() const
The square of the cosine of the weak mixing angle defined without weak radiative corrections...
void computeHpquantities()
double THoEX_pp_H_hh_bbbb_ATLAS13
gslpp::matrix< double > log_cs_WH_8
double ip_ex_pp_H_hh_CMS8_cache_e[2][CacheSize]
double ie32ec[2][CacheSize]
An observable class for the quartic Higgs potential coupling .
double THoEX_mu_pp_H_VV_CMS8
gslpp::complex B00_MZ2_MW2_mHh2_mHp2(const double MZ2, const double MW2, const double mHh2, const double mHp2) const
double ip_ex_gg_H_ZZ_llllnunu_ATLAS13_cache[2][CacheSize]
double THoEX_pp_Hp_tb_ATLAS13_1
gslpp::matrix< double > temp25e
gslpp::complex B0_MZ2_0_mHp2_mHl2(const double MZ2, const double mHp2, const double mHl2) const
double ip_ex_pp_Hp_tb_CMS8_cache_e[2][CacheSize]
double ip_ex_gg_phi_gaga_CMS13_cache_e[2][CacheSize]
gslpp::matrix< double > temp1e
double ip_ex_gg_H_WW_ATLAS8_cache[2][CacheSize]
gslpp::matrix< double > ATLAS8_gg_A_hZ_bbZ_e
double THDM_BR_h_ZZ
branching ratio to two bosons in the THDM.
gslpp::complex A_HH_D(const double mHh2, const double cW2, const double Ms, const double Mb, const double MZ) const
double ip_cs_ggtoA_8(double mass)
Interpolating function for the A production cross section via gluon-gluon fusion at 8 TeV...
virtual double getMHl() const
A get method to retrieve the Higgs mass .
double ip_cs_pptobbA_13(double mass)
Interpolating function for the bottom associated A production cross section at 13 TeV...
double rh_ll
Squared relative coupling of to two charged leptons.
double THoEX_pp_H_Zga_llga_CMS8
double rA_gg
Squared relative coupling of to two photons.
double ip_ex_pp_Hpm_taunu_CMS13(double mass)
Interpolating function for the observed CMS upper limit on a singly charged scalar resonance decaying...
gslpp::matrix< double > CMS13_gg_H_ZZ_llnunu_e
THDM branching ratio of .
double ip_cs_pptottH_8_cache[2][CacheSize]
double ip_ex_gg_phi_tautau_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double ip_ex_gg_phi_Zga_llga_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double ip_ex_pp_Hpm_taunu_ATLAS13_e(double mass)
Interpolating function for the expected ATLAS upper limit on a singly charged scalar resonance decayi...
double ip_ex_pp_H_hh_bbbb_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a spin-2 resonance decaying to two boso...
gslpp::complex B0p_MZ2_mHp2_0_mHl2(const double MZ2, const double mHp2, const double mHl2) const
gslpp::matrix< double > ATLAS8_pp_Hpm_tb_e
double ie13c[2][CacheSize]
gslpp::matrix< double > CMS13_VBFVH_H_ZZ_llll
double ip_ex_pp_Hpm_taunu_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a singly charged scalar resonance decaying...
gslpp::matrix< double > temp17e
double rh_VV
Squared relative coupling of to two massive vector bosons.
gslpp::matrix< double > CMS8_gg_H_hh_bbtautau_e
double ip_ex_VBF_H_ZZ_ATLAS8_cache[2][CacheSize]
gslpp::matrix< double > ATLAS8_gg_phi_tautau_e
double ggF_H_WW_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
double ie1ec[2][CacheSize]
double ip_csr_ggH_t_8(double mass)
Interpolating function for the gluon-gluon fusion H cross section ratio of the top-loop and the total...
gslpp::matrix< double > ATLAS8_gg_phi_tautau
gslpp::matrix< double > temp35
double ip_ex_pp_phi_Zga_qqga_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to a boson an...
gslpp::complex B0p_MZ2_mHp2_mHp2_mHl2_cache[4][CacheSize]
double ip_ex_gg_A_Zh_Zbb_ATLAS13_cache[2][CacheSize]
gslpp::matrix< double > ATLAS8_gg_H_WW
double ie21e(double mass)
gslpp::complex B0p_MZ2_0_0_mHl2(const double MZ2, const double mHl2) const
double ip_ex_pp_H_ZZ_llqq_CMS13_cache[2][CacheSize]
double computeThValue()
The quartic coupling .
gslpp::matrix< double > log_cs_ZH_8
gslpp::matrix< double > temp2e
double R_ggF_H_hh_bbbb_CMS13
gslpp::matrix< double > ATLAS8_VBF_H_ZZ_e
double THoEX_ttF_H_tt_ATLAS13
double ip_ex_ggF_phi_Zga_CMS13_cache_e[2][CacheSize]
double R_pp_H_hh_bbgaga_CMS13
double ip_cs_pptobbH_8_cache[2][CacheSize]
double ie40c[2][CacheSize]
gslpp::matrix< double > temp17
double ip_ex_gg_H_WW_lnuqq_ATLAS13_cache[2][CacheSize]
double bbF_A_hZ_bbZ_TH13
Cross section times branching ratio for the process at the LHC with 13 TeV.
double THDM_BR_h_tautau
branching ratio to two leptons in the THDM.
double ip_ex_gg_H_ZZ_qqllnunu_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
gslpp::matrix< double > temp7e
gslpp::complex B00_MZ2_0_MZ2_mHl2(const double MZ2, const double mHl2) const
double ip_ex_gg_H_ZZ_llqq_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double VBF_H_WW_TH8
Cross section times branching ratio for the process at the LHC with 8 TeV.
gslpp::complex B0_MZ2_0_MZ2_mHl2_cache[3][CacheSize]
gslpp::complex B0p_MZ2_mHh2_0_mA2_cache[4][CacheSize]
double ip_ex_pp_phi_Zga_qqga_CMS13_cache_e[2][CacheSize]
double ip_ex_bb_A_Zh_Zbb_ATLAS13_cache_e[2][CacheSize]
gslpp::matrix< double > CMS8_pp_H_hh_gagabb
double ip_ex_pp_H_hh_gagaWW_ATLAS13_cache_e[2][CacheSize]
double THoEX_pp_H_hh_bbVV_CMS13
double ip_ex_VBF_H_ZZ_ATLAS8_e(double mass)
Interpolating function for the expected ATLAS upper limit on a vector boson fusion produced scalar re...
double ip_ex_pp_H_ZZ_llnunu_CMS13_e(double mass)
Interpolating function for the expected CMS upper limit on a scalar resonance decaying to two bosons...
double ip_ex_pp_Hp_tb_ATLAS13_2_e(double mass)
Interpolating function for the expected ATLAS upper limit on a singly charged scalar resonance decayi...
gslpp::matrix< double > ATLAS13_pp_Hpm_taunu
ATLAS observed upper cross section limits at 13 TeV, depending on the charged Higgs mass...
double ip_ex_gg_H_ZZ_llnunu_CMS13_cache[2][CacheSize]
double ip_cs_WtoWH_8(double mass)
Interpolating function for the W associated H production cross section at 8 TeV.
double ip_ex_pp_phi_gaga_ATLAS8_cache_e[2][CacheSize]
gslpp::matrix< double > CMS8_gg_phi_gaga_e
gslpp::matrix< double > temp20
double ip_ex_VBF_H_WW_lnuqq_ATLAS13_cache_e[2][CacheSize]
double interpolate(gslpp::matrix< double > arrayTab, double x)
Linearly interpolates a table with one parameter dimension.
double ip_ex_bb_phi_tautau_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a bottom quark produced scalar resonance...
double ip_cs_pptottA_13(double mass)
Interpolating function for the top associated A production cross section at 13 TeV.
gslpp::complex B0_MZ2_0_mA2_mHh2_cache[4][CacheSize]
gslpp::complex I_HH_D(const double mHh2, const double Ms, const double Mb) const
double ip_ex_gg_H_hh_ATLAS8_cache[2][CacheSize]
double ip_ex_gg_H_WW_ATLAS8(double mass)
Interpolating function for the observed ATLAS upper limit on a gluon-gluon produced scalar resonance ...
double THoEX_bbF_H_bb_CMS8
gslpp::matrix< double > ATLAS13_gg_phi_Zga_llga_e
double ip_ex_VBF_VH_H_ZZ_llll_CMS13_cache[2][CacheSize]
gslpp::complex I_h_U_cache[5][CacheSize]
gslpp::matrix< double > CMS8_pp_Hp_tb_e
double ip_csr_ggH_t_8_cache[2][CacheSize]
double THoEX_ggF_A_tautau_ATLAS13
gslpp::matrix< double > temp35e
double ie24c[2][CacheSize]
gslpp::complex B0_MZ2_mHh2_mHl2_mHl2_cache[4][CacheSize]
gslpp::matrix< double > log_cs_ggA_8
CP-odd Higgs production cross section tables at 8 TeV obtained with HIGLU 4.34, depending on the Higg...
double THoEX_bbF_H_tautau_ATLAS13
gslpp::matrix< double > CMS13_pp_H_ZZ_llnunu_e
complex sqrt(const complex &z)
double ip_ex_bb_A_Zh_Zbb_ATLAS13(double mass)
Interpolating function for the observed ATLAS upper limit on a bb associated pseudoscalar resonance d...